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Ce travail a ete consacre a I'etude de I'effet Hall dans deux systemes differents de 
fermions fortement correles: un systeme constitue de liquides de Luttinger faiblement 
couples et dans un reseau triangulaire bidimensionnel. Dans le but d' avoir les outils 
necessaires pour attaquer ces problemes, on a entrepris une revision des proprietes 
principales des systemes avec des interactions fortes, en nous concentrant dans les 
systemes a deux dimension et a basse dimensionnalite (ID et quasi ID). Etant donne 
que I'effet Hall est avant tout un phenoneme de transport, on a etudie les formalismes 
existants dans la litterature pour traiter les proprietes de transport dans les systemes 
fortement correles, et nous avons consacre tout un chapitre au traitement specifique de 
I'effet Hall. 

Done en un premier temps nous avons etudie le modele des liquides de Luttinger 
faiblement couples. L' etude de I'effet Hall dans un systeme quasi ID a ete partielle- 
ment motivee par les resultats experimentaux obtenus sur des conducteurs organiques 
quasi ID. Ce travail etait base sur 1' application du formalisme de matrice memoire 
pour obtenir la constante de Hall (Ru) dans des chaines a moitie remplis faiblement 
couplees, incluant la diffusion umklapp. La geometric du modele a ete choisie pour 
avoir le courant circulant le long des chaines ID. On a calcule la dependance en tem- 
perature et frequence de Ru en tenant compte des interactions particule-particule et 
particule-reseau (processus umklapp). On a obtenu un coefficient de Hall Ru con- 
stitue d'abbord par un terme d' electrons libres (valeur de bande R^) plus un terme 
de corection avec une dependance en temperature (frequence) donne par une loi de 
puissance, due a la presence de la diffusion umklapp. Ces lois des puissance sont des 
signatures du comportement des liquides Luttinger, oil les exposants dependent des 
parametres d'interactions. Le coefficient de Hall a aussi ete calcule dans le systeme 
sans interaction particule-particule, condisuant a une dependance logarithmique en T 
(oil 00), en accord avec la limite d'interaction nuUe de la loi de puissance. Etant donne 
que les conducteurs organiques quasi ID sont en meme temps des systemes a moitie 
remplis et au quart remplis, nos resultats theoriques ne sont pas directement appUca- 
bles aux messures d'effect Hall existantes, mais ils nous ont permis d'arriver a des 
conclusions interesantes. En premier lieu, la fagon appropiee d' analyzer les messures 
d'effet Hall realisees dans ces sytemes quasi ID (dans la meme geometric du modele), 
est de fitter les deviations de la valeur de bande 
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CHAPTER 1 



Introduction 



It is likely that Edwing Herbert Hall was not aware of the enormous impact that his 
experimental results would have when he measured the Hall voltage for the first time 
in 1879. Not only because the Hall effect has allowed to determine the sign and density 
of the charge carriers of many materials, but also because it lead to a whole new field 
of research on the nature of the particles responsible for the electrical current in any 
system. 

Nowadays, the Hall resistivity has become a common tool when studying the trans- 
port properties of any new material. However real experiments have shown that this 
phenomenon cannot always be interpreted in terms of the density of charge carriers 
and thus, has to be interpreted in a different way. Compounds can have strong interac- 
tions affecting the transport properties. Such strongly interacting systems are now at 
the heart of research in condensed matter physics and the understanding of the effect 
of these interactions is crucial for the description of these materials. 

In this work we undertake the challenging task of studying the Hall effect in two 
different systems taking into account interactions. This requires the study of the exist- 
ing methods for treating strongly interacting systems and their application to the Hall 
effect. We have chosen two model systems which are similar to real compounds where 
the measured Hall effect cannot be readily explained. The work presented here is en- 
tirely theoretical but some applications to real experiments are discussed. The plan of 
the work is presented in what follows. 

We begin this manuscript with a description of the standard theoretical models used 
to describe two-dimensional and low-dimensional fermions systems with strong cor- 
relations. In particular, we review the main properties of the two-dimensional square 
lattice and the 2D triangular lattice, and then those of one-dimensional and quasi one- 
dimensional systems. 

The second chapter is devoted to the transport phenomena in systems with strong 
interactions. For this we explain the linear response theory and the Kubo formulas. 
Such formulas are necessary for the understanding of the memory function formalism 
developed afterwards. We end this chapter with an application of the previously men- 
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tioned theories in the calculation of transport properties in low-dimensional systems. 

The third chapter is completely dedicated to the Hall effect. We begin with the 
explanation of the classical Hall effect. Next, we obtain an expression for the Hall 
resistivity at infinite frequency and finally, we use the memory formalism to include 
interactions in the expression of the Hall resistivity. 

Our first theoretical contribution is presented in the fourth chapter. The Hall effect 
is investigated in a quasi one-dimensional system made of weakly coupled ID chains. 
For this we begin with a review on quasi ID organic conductors, specifically in Bech- 
gaard and Fabre salts, where various Hall measurements have been accomplished. We 
explain the model and methods used to calculate the Hall constant. We then present 
the results and discuss their range of validity and application to real experiments. 

In the fifth chapter we present our second theoretical work. In this case we con- 
sider the two-dimensional triangular lattice as our model system. We start with the 
description of the sodium cobalt oxide, a compound with a triangular lattice structure, 
where the Hall effect have been investigated recently. As before, we introduce the 
model and methods used to calculate the Hall constant. We present the results ob- 
tained and discuss their range of validity and application to the previously mentioned 
experiments. 

We close this manuscript with some general conclusions from the two theoretical 
works presented and perspectives for future works. Some appendices are presented at 
the end. 
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CHAPTER 2 



Strongly correlated systems 



In condensed matter physics, understanding the effects of interactions in real materi- 
als, has occupied theoretical and experimental physicists for more than fifty years now. 
Depending on the strength of interactions the convenient approach to be implemented 
can be completely different. For example, if interactions are very weak, like in metals, 
the free-electron model works well to describe most of the physical properties. The 
dimensionality also plays a crucial role in this choice. In general, high-dimensional 
interacting systems are well described by the Fermi Liquid theory, where electrons 
are dressed by the density fluctuations around them, forming individual objects called 
quasiparticles , as formulated by Landau and others HIEIISI. The fraction of the elec- 
trons that remain in this quasiparticle state is called the residue Z, and gives the ampli- 
tude of the discontinuity in the occupation factor nfc. The big result of the Fermi liquid 
theory is that the properties of the interacting system remain essentially the same of 
the free fermionic particles. 

In one dimension, on the contrary, the Fermi Liquid theory breaks down and is the 
Luttinger Liquid theory that describes the physics, even in systems with very strong 



interactions, as explained in Sec. 2.4. 1 Furthermore, Fermi Liquid theory is not al- 
ways applicable in high-dimensions when one is working with strong interactions. 
The common characteristic of strongly correlated systems, where one finds materials 
as different as heavy fermions, high Tc superconductors, cobalt oxides and organic 
conductors, is a narrow band where fermions interact strongly on a scale of the order 
of electronvolts at short distances. In two-dimensional systems, for example, none of 
the previously mentioned theories is applicable when strong interactions are present. 
Then, basic models such as the Hubbard, t-J, Anderson and Kondo lattice models had 
to be implemented in an effort to understand the properties of these strongly correlated 
layered materials. 

As we are interested in the study of the Hall effect in a quasi one-dimensional 
system and in a two-dimensional triangular lattice, we will devote this chapter to the 
description of the theoretical models used to describe these systems throughout this 
work. First, we will study the main properties of the two-dimensional Hubbard model. 
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Figure 2.1: The Hubbard model: electrons hops from site to site with a hopping ampli- 
tude t. Full (empty) circles denote sites with one (zero) electron. When a particle hops 
into an occupied site, it has to pay a penalty U (on-site Coulomb interaction). This is 
shown in the image on the right. 



focusing on two geometries: the rectangular and the triangular lattice. Then, for the 
case of low-dimensional systems, we will make a short review on the Luttinger liquid 
theory and its extension to describe quasi one-dimensional systems. 



2.1 The Hubbard model 

The Hubbard model is one of the simplest models existing in literature Q to describe 
the physics of interacting fermions on a lattice. Hubbard (1963) proposed an extension 
of the tight-binding model [2], where electrons hop from site to site with a matrix 
element t, adding a term that provides a penalty U for any atomic site occupied by 



more than one electron. This is depicted in Fig. 2.1 In terms of fermionic operators. 



the Hubbard Hamiltonian takes the following form O 

W = X] ^^i^^j^ + U^ni^nii- rii^, (2. 1) 

{ij)a i icr 

where Co (cq,) is the creation (annihilation) fermion operator, = CaCa is the fermionic 



number operator and (ij) are nearest-neighboring sites. The first term in Eq. (2. li 
represents the kinetic energy of the electrons and the second term correspond to the 
interaction between them. In addition, we include a chemical potential /x in order to 
adjust the number of electrons in the system. 

There are extensions of the Hubbard model in which electron hopping beyond 
nearest-neighbors and electron-electron interaction at larger distances, are taken into 
account. We will not treat these models here because they are beyond the scope of this 
work, and we refer the reader to the literature for a complete description of these cases 
El . The Hubbard model was solved exactly in one dimension, by Lieb and Wu (1968) 
m , as we will discuss it in the next section. For dimensions greater than one, many 
attempts have been undertaken to derive a phase diagram at zero temperature using 
random phase approximation (RPA) 121, mean field theory [SJ, numerical methods 
as the dynamical mean field theory (DMFT) (6], etc; but there are still many open 
questions and discrepancies between different approaches. 



4 



When working with the Hubbard model, there are two important energy scales to 
take into account: the on-site Coulomb interaction U and the bandwidth W. They 
determine the behavior of the system by the value of the ratio U /W. If this ratio is 
small, the free-electron model can be used to describe most of the physical properties 
of the system. Simple metals like aluminium or sodium have U <^ W. On the other 
hand if the ratio is large, then U ^ W and electron correlations are dominant. This 
case is not fully understood but some of its properties are currently well characterized. 
We list some of them below. 

The electron density n = N/Ns (N being the number of electrons and Ng the 
number of lattice sites) varies from to 2, due to Pauli principle. When n = 2 there 
are two electrons in each site, they cannot move, and the system is an insulator. At 
half filling, when there is exactly one electron per site (n = 1), the ground state is 
antifenomagnetic (AF), i.e., neighboring spins pointing in opposite directions. Fur- 
thermore, if U ^ W and there is one particle per site, the system is in an insulating 
state known as a Mott insulator ||71|8l. Possible phases of the Hubbard model are: the 
paramagnetic phase, the ferromagnetic phase, the charge-density wave phase (CDW) 
and the spin density wave phase (SPW) 0. Each of these phases has particular prop- 
erties which have been studied extensively in the last decades. Nevertheless a unique 
phase diagram for the Hubbard model for dimensions d > 1, is still lacking. 

One important property of the Hubbard model at half-filling is the particle-hole 
symmetry, which means that the Hamiltonian in Eq. ( 2. 1 1 is invariant under the electron- 
hole transformation: Cjcr —>■ {—l)^dj^, where da (da) creates (annihilates) a hole. This 
property is very peculiar to half-filled bi-partite lattices (lattices that can be decom- 
posed into a disjoint union of two sublattices) and is not fulfilled, for example, on the 
triangular lattice (see Sec. 2.3 1. 

An special limit of the Hubbard model that will be studied in Chapter. |6] is called 
the atomic limit and corresponds to a system where the hopping amplitude is set to 
t = 0. In this case the system consist on isolated sites where particles cannot move 
from site to site. An exact solution can be found for the atomic limit of the Hubbard 
model (each site is considered individually), for any dimension (see Appendix. B.3 1. 
Another exactly solvable limit is the U = limit (tight-binding model), where the 
Hubbard Hamiltonian can be also completely diagonalized. The properties of the Hub- 
bard model strongly depend on the type of the lattice and on the dimensionality. That is 
why we will focus here only on two geometries of the 2D Hubbard model: the square 
lattice and the triangular lattice. But, before entering into these two particular geome- 
tries, we will mention another variant of Hamiltonian ( |2.1[ ) known as the t-J model 
0. 

The t-J model describes a system of interacting fermions on a lattice where the 
on-site interaction U is considered infinite, in order to prohibit double occupancy, i.e., 
sites with nj^njj^ = 1. This constraint is implemented using the Gutzwiller projection 
operator Pq = nj(l — nj^njj). Applying this operator to Hamiltonian (2.1 1 one 
obtains Ht-j = -Pg +"exchange term", where the exchange term is given by 
J '^i ' "^i' with Si the spin operator at site i. The exchange term represents an 
effective antiferromagnetic interaction between neighboring spins, which results from 
the virtual hopping of electrons into and out of doubly-occupied states, and has an 
amplitude J = At^ /U. The t-J model has no exact solution either, but it is widely 
used in the study of fermionic systems with very strong interactions. We will not 
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study the properties of the t-J model here because we are interested in the effect of 
interactions as a function of their strength, and for this the Hubbard model is the more 
accurate model to describe our two-dimensional systems. But it remains a good model 
to make comparisons at very strong interactions, as will be done in Chapter[6] We refer 
the reader to the literature for a complete review on the t-J model lITOlfTTI . 



2.2 The 2D square lattice 

In this section we will study some general properties of the square lattice and defi- 
nitions that will be useful in the rest of this work. Then, we will give the properties 
of the Hubbard model in this particular geometry. The two-dimensional square lattice 



has the structure shown in Fig. 2.1 which is invariant under translation (it has trans- 
lational symmetry). The lattice is defined by two primitive vectors ai = a(l, 0) and 
0-2 = a(0, 1), where a is the intersite distance, called the lattice parameter. The real 
space lattice, often referred as direct lattice, can also be represented in momentum 
space. The latter is known as the reciprocal lattice and is formed by the set of vectors 
K satisfying e*^'^ = 1, with R a vector of the direct lattice Li2 |. The two primitive 
reciprocal vectors are bi = ^(1, 0) and 62 = ^(Oj !)• We will see, in a minute, how 
any vector on k-space is defined by bi and 62- Due to the periodicity of the lattice, 
it properties can be studied in a reduced zone or primitive cell of the reciprocal space 
defined by the Wigner-Seitz cell of the origin, which associates with each lattice point 
all of space which is closer to it than any other lattice point |5]. This is the Brillouin 



zone (see Fig. 2.2 1. 



One important quantity that changes from one geometry to the other, is the disper- 
sion relation or momentum energy of the particles. It also depends on the model used. 
Let us take the previously studied Hubbard model to describe the 2D square lattice. 
The dispersion relation is obtained by transforming the Hamiltonian to momentum 
space, using the Fourier transform of fermionic operators 

where Q is the volume of the system. The inverse Fourier transform is thus given by 
^fccr ~ {^/V^^ Yli ^''"^'c|cr- We will also need the Fourier transform of the number 
operator, which is 

k 



Then, Fourier transforming Hamiltonian (|2.1|) we obtain Q 



W = ^ efccLcfc, + UY^ nk]n_ki - l^Yl ^fc-Cfc,.- (2-4) 



fccr fc fccr 



The coefficient of the first term in Eq. (2.4i (the kinetic energy term) gives the disper 



sion relation of the 2D square lattice. Thus, the dispersion relation is 

fffc = —2t [cos{kxa) + cos(A;ya)] , (2.5) 
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Figure 2.2: (Left) Fermi surface of the 2D square lattice. The red square represent the 
first Brillouin zone for the square lattice and the Fermi surface is the shaded area inside 
this zone. (Right) Density of states of the square lattice. The DOS exhibits a van Hove 
singularity at the center of the band. The bandwidth of the system is equal to 8|t|. 



where a is the lattice parameter (|rj — rj\ = a, with i,j nearest-neighboring sites). In 
Eq. (2.5 1 we supposed an isotropic square lattice. The components and ky of any 
vector in the reciprocal lattice are obtained from the primitive vectors bi and &2 trough 
the relation: kx = nbi/lir and ky = mb2/2vr, where m and n are integers. The 
dispersion relation determines the bandwidth {W) of the system, which is given by the 
range of energies between the maximum and minimum of e^. As can be seen from 
Eq. (2.5 1, the bandwidth of the 2D square lattice is equal to = 8|t|. The dispersion 
relation is often written with reference to the chemical potential, combining the 1st 
and 3rd terms on the right hand side of Eq. (2.4 1: ^fc = —2t [cos(A;a;a) + cos{kya)] — 



/i. The ground state of the non-interacting system ([/ = in Hamiltonian \2A) ) is 
constructed by occupying all the energy-levels with two electrons with opposite spins, 
until reaching the total number of particles N. The highest energy level is called the 
Fermi level and has an energy e^p, with /cp the Fermi momentum. This energy level 
allows to define the "Fermi surface", which is an important quantity when studying real 
materials (nowadays measured experimentally with different methods, such as angle 
resolved photoemission spectroscopy). The Fermi surface is defined by the collection 
of points in the reciprocal space with energy e^p = fi, i.e, = 0. Fig. 2.2 shows 
the Fermi surface and first Brillouin zone of the 2D square lattice described with the 
Hubbard model. Following the Luttinger Theorem, interactions do not change the 
volume of the Fermi surface, only its shape. Another important quantity is the Density 
of states (DOS) of the system, which describes the number of states at each energy 
level that are available to be occupied. It is important to known the shape of the 
DOS when one is working with the Hall effect, due to the relation between the Hall 
resistivity and the electronic density, as we will see in Chapter |4] Fig. 2.2 shows the 
DOS of the non-interacting square lattice, which exhibits a van Hove singularity at the 
center of the band. 

The particle-hole symmetry mentioned in the previous section is fulfilled in the 
square lattice, because it is a bipartite lattice. The particle-hole symmetry at half- 
filling is evident from the shape of the DOS, which is symmetric with respect to the 
center of the band. Although we will not work on the square lattice for the Hubbard 
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a) 



b) 





Figure 2.3: a) Two-dimensional triangular lattice structure with hopping amplitudes t 
and t'. b) Geometrical frustration on a triangle. Should the empty site be spin-up or 
spin-down?. From Ref. lfT3l 



model, it is important to know this geometry because it is the simpler model, it has 
been extensively studied and thus, it always works as a point of reference. Let us now 
move to the triangular lattice, which is the geometry that will be used with the Hubbard 
model in Chapter. [6] 



2.3 The 2D triangular lattice 



The structure of the two-dimensional triangular lattice is shown in Fig. \2.3\ i. In the tri- 
angular lattice, hopping t occurs along the longitudinal direction and hopping t' along 
a direction defined by an angle of 60 degrees from the longitudinal direction, as can be 



seen in Fig. 2.3 1. The lattice is called isotropic when the hopping amplitudes satisfy 
t = t'. The primitive vectors defining the lattice in real space are ai = a(l,0) and 
a2 = a(l, \/3)/2, with a the lattice parameter. And for the reciprocal space vectors 
we have bi — '^'^ ' 



^(1 



1 

V3' 



and &2 = ^(0, ■^)- As before, we obtain the dispersion 



relation from the Fourier transformed Hamiltonian ( 2.4 1 for the 2D triangular lattice. 
Thus, efc is this geometry is given by 



tcos{kxa) + 2t' cos(/c^a/2) cos{kyaV3/2] 



(2.6) 



Thus, the bandwidth is W = 9\t\. With this dispersion relation we obtain the Fermi 
surface shown in Fig. |2.4[ as well as the Brillouin zone which has an hexagonal char- 
acter. The corresponding density of states (DOS) exhibits two van Hove singularities 
which are degenerate when t = t' (see Fig. 2.4 1. Unlike in the square lattice, the 
DOS has no particle-hole symmetry, irrespective of the value of t and t'. Thus, the 
particle-hole symmetry is not fulfilled because the lattice is no a bipartite one. 

Two important properties arise from the topology of the triangular lattice. The 
first one is that it has the smallest number of steps necessary to make a loop, namely 
three. This will be reflected in the study of the Hall effect on the triangular lattice in 
Chapter [6] The second one is called the electronic frustration and it appears because 
the geometrical arrangement in the triangular lattice frustrates the alternation of the up 



and down spins that minimizes the Coulomb repulsion. As shown in Fig. 2.3 two of 
the three electrons in each triangle share the same spin orientation 1.13,1 . The effects 
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Figure 2.4: .(Left) Fermi surface of the 2D triangular lattice. The green hexagon represent 
the first Brillouin zone for the triangular lattice and the Fermi surface is the shaded area 
inside this zone. (Right) Density of states of the triangular lattice. The bandwidth of the 
system is W — %\t\. The density of states (DOS) exhibits two van Hove singularities 
which are degenerate when t — t'. 



of frustration in transport properties are not yet fully understood. We will not treat 
frustration here because it is out of the scope of this work. We refer the reader to the 
literature |[T3l for an explanation on this phenomena. The ground state of the triangular 
lattice is expected to be antiferromagnetic with the spins oriented forming 120 degrees 
between them. 

At this point we have seen numerous properties of the two-dimensional Hubbard 
model in two different geometries. In particular, the properties of the triangular lattice 
will be used in Chapter. [6] In the following, we will study other types of systems called 
low dimensional systems where the physics vary substantially from the one studied 
previously. 

2.4 Low dimensional systems 

Until now we have seen the theoretical models, that will be used in the next chapters, to 
describe strongly correlated systems in two dimensions. The next step, is to investigate 
the rest of the systems that interest us for this work, the one- and quasi one-dimensional 
cases, in order to study their most important properties. In particular those related to 
transport phenomena. We will begin this section by studying the main features of 
electrons moving in one dimension, and the Luttinger liquid theory that allows us to 
describe the physics of these systems. Then, we will focus on a special ID system 
called the one-dimensional Mott insulator. It is essential to understand the physics of 
the latter if one is interested in electrons moving on a lattice, as will be the case in 
Chapter |5] Finally, we will investigate the main properties of quasi one-dimensional 
systems, and the way to extend the Luttinger liquid theory with the aim of describing 
their physics. 

2.4.1 One-dimensional case 

In order to understand the physics and, in particular, the transport properties of strongly 
correlated systems in one dimension, we must first realize that fermions in ID have 
quite different behaviors than fermions in higher dimensions. In one dimension, an 
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Figure 2.5: (a) Single-particle dispersion in a one-dimensional system. The Fermi sur- 
face consists of two points, (b) Linear spectrum of the Luttinger model. The particles are 
now separated in right (R+) and left (L— ) going fermions. The spectrum is extended to 
±00, giving rise to an infinite number of negative states ('unphysical'states). This requires 
the introduction of a cutoff on the momentum to make the model well defined. 



individual electron cannot move without pushing its neighbors, thus any individual 
excitation has to become a collective one. The main consequence of these differences 
is that, instead of the Fermi liquid theory, the proper theoretical description for one- 
dimensional systems is the Luttinger Liquid theory. In the following we will focus on 
the basic ideas necessary to describe transport properties in one-dimensional systems, 
and we refer the reader to the many existing reviews and books in the literature lfT4l 
[15] [161 for a detailed description of all the properties of ID systems. 



Fig. 2.5 a) shows the single -particle dispersion for a one-dimensional system, where 
the Fermi surface consists of only two discrete points. Because the states contributing 
to transport are the low-energy ones, i.e. excitations close to the Fermi level, we may 
replace the original model by one with a purely linear spectrum (see Fig. |2.5fc )) ex- 
tended to ±00. This linearization was first proposed by Luttinger [ 1 8 1 and requires the 
introduction of two species of fermions, right (R) and left (L) going fermions with a 
dispersion relation given by 

e±{k) = VF{±k-kF), (2.7) 

where is the Fermi velocity, kp is the Fermi momentum and + (— ) refers to right 
(left) going fermions. 

Another property of one-dimensional systems is that low-energy particle-hole ex- 
citations, where an electron with momentum k is removed from below the Fermi sur- 
face and created above with momentum k + q, do not exist for < g < kp. This 
arises because the Fermi surface consist of only two points and thus the only places 
where the particle-hole energy can reach zero are q = and q = 2kp. This is de- 



picted in Fig. 2.6 Thus, these low-energy excitations have a well defined momentum 
q, a well-defined energy E{q) q (independent of k) and an energy dispersion SE{q) 
which goes to zero much faster than the average energy. These properties make them 
well defined 'particles 'with a lifetime that increases when one goes closer to the Fermi 
level {E = 0) |[T6l . In the linearized model, these excitations have an energy (for right 
going fermions) 

ER,k{q) = vv{k + q) - vpk = Vpq, (2.8) 
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Figure 2.6: Particle-hole spectrum for high-dimensional systems (a) and one- 
dimensional ones (b). In ID systems, particle-hole excitations have well defined energy 
and momentum for small values of q (low energy modes at q ~ and q ~ 2k p), which is 
not the case in high-dimensional systems (From Ref. I,17j ). 



which is indeed independent of k. They have a well defined momentum q and energy 



E{q) = v^q. The particle-hole spectrum for the Luttinger model is shown in Fig. 2.7 
These excitations are bosonic in nature (product of two fermion operators) and because 
they are well defined, they can be used as a basis to represent the Hamiltonian and 
all fermionic operators. Indeed, the density fluctuations which are a superposition of 
particle-hole excitations, are used as such basis. 

To further understand the previous arguments, let us take a system of fermions 
moving in ID and write down the linearized kinetic energy term in momentum space 

Wkin= vv{±k - kv)c{f^^c^.^^, (2.9) 

k,a,r=R,L 

where the + (— ) sign corresponds to r = R{r = L), right (left) going fermions. To 
transform this term to direct space, the single particle fermion operator tl^a{x) must be 
rewritten in terms of right and left going fermions, taking into account only the parts 
acting close to the Fermi surface 

^'^(^^ = ^Xl^'''^^'^'^ - ^ ^'''"'^'"'^^ e^'^^cfc^ = V'<7ii(2;) + V'aL(x), 

k k^+kp kr^—kp 

(2.10) 

where i7 is the volume of the system. Thus, in terms of ijjR and ijjL, the kinetic energy 
term in direct space is given by 



Wkin = j dx 'Y ^VF'lpl{x)T3{-idx)lp^{x) 



(2.11) 



/,t „/,t 



where V'^ = (V*/? ^I) is the two-component vector composed of right- and left-moving 



fermions and T3 is the Pauli matrix (g i'^)- Eq. (2.11 1 is easily Fourier transformed 
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Figure 2.7: Particle-hole spectrum in the Luttinger model. The states on the left of the 
black line are 'unphysical'states coming from the extension to ±00 of the single-particle 
spectrum (see Fig. |2.5| l. That is why a cutoff must be introduced in momentum to make 
the model well defined 



into Eq. ( |2.9| l using representations (2.10 1. It can be proved, by simple commutation 
relations fT6l, that the kinetic energy term < \2.9} is well represented in a basis generated 
by boson operators bj, and bp'. Tikin — X^p^o ^^f|p|&^^ , showing how this term can be 
quadratic in such a basis. Similarly, in a system of fermions with spins, the interaction 
term has typically the following form 

^int = J dx dx'V{x — x')pa{x)p(ji{x'), (2.12) 



where Pa{x) is the density operator defined as Paix) = 'il)l-{x)il^a{x). This operator is 
made of a product of two fermions operators and thus it can also be represented in the 
boson basis p{q) bg or p{q) ~ bg (or some linear combination). The interaction term 
thus remains quadratic in the boson basis and can be easily diagonalized, as shown 
below. In terms of right and left going fermions the density operator becomes 



(2.13) 

It must be kept in mind that the relevant processes in the interaction are those close 
to the Fermi surface (g ~ and q ~ 2A;f). Because the Fermi surface in ID consist 
only on two points, the interaction term can be decomposed in three different types. 
The first process called couples fermions on the same side of the Fermi surface. 
The second type is called g2 and couples fermions from one side with fermions on the 
other side of the Fermi surface. But each species stays on the same side after the inter- 
action (forward scattering). And finally we have the gi process which corresponds to a 
scattering of 2fcF> where fermions exchange sides (backscattering) |[T6]| . The processes 
between fermions with parallel spins are denoted g\\ and with opposite spins, g±. We 
suppose our system is made of fermions with spin rotational symmetry (for simplic- 
ity), moving in one dimension. Adding the interactions terms to the Hamiltonian we 
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obtain the following expression 



WlD = j dx ^ \vFijl{x)T3{-ida:)'llj^{x) +52]! PaR{x)PaL{x) (2.14) 



+g2± PaR{x)p-aL{x) + ^ Par{x)par{x) + ^ Par{x) P-ar{x) 



r=R,L 



The first term is the kinetic energy ( |2.1 1| ) and the other terms in Eq. (2.14i refer to 



forward scattering. In Hid we have omitted the backscattering terms (gi processes) 
which are, for spin rotationally invariant systems, marginally irrelevant lfT6l . We there- 
fore take gi± = g^ = 0. If only backscattering terms are considered, it must be kept 
in mind that logarithmic corrections can be introduced by the forward scattering terms. 
This will be the case in Chapter [5j in our study of the Hall effect in a quasi ID chain. 

The representation of fermionic operators in term of bosonic fields, is an important 
tool of the Luttinger liquid theory to calculate correlations functions, as explained at 
the end of this section. This technique is known as bosonization |[T5l[T6l . We will not 
present the bosonization method in details here, but only the important results needed 
for the understanding of the transport properties in ID. We refer the reader to the 
literature for a detailed description of this technique |[T5l[T61 . Hamiltonian ( 2.14| ) will 



be used in Chapter [5j both in its fermionic and bosonic representations (see Eq. ( 2.20| )) 
Let us see first how the fermionic fields ipfj = ipa,R + V'ct.l represented in 
terms of bosonic fields denoted 9y and (fields written in the boson basis h and 6^^), 
where v = p{a) denotes the charge (spin) degrees of freedom. The charge and spins 



degrees of freedom must be separated in order to diagonalize Hamiltonian ( 2. 14 1. The 
fermionic field yields 

^ _ ^ ^-^{r(l>p{^)-9p{x)+a[r<l)„(x)-9„(x)]} 

with r = 1) for right (left) moving fermions, and a a cutoff which precise value 
is irrelevant for the low-energy, long-wave length properties of the model. The limit 
a — > should be taken in principle. All the formulas given here are found after taking 
the thermodynamic limit L —>■ oo, with L the size of the system |[T5l[T61 . The fields 
{(pp^Op) and {(t)tj,9^) obey the following commutation relations 



.(x),-V0^(x') 
vr 



i5{x-x'). (2.16) 



The bosonic fields cj) and can be represented in terms of the density operators p^ and 
PL llT6il . giving rise to the following relations 

V</>|(x) = -TT [pR^{x) + PL^{x)\ 

Ve^ (x) = -vr [pR^ (x) - pL^ (x)] , (2.17) 



The Nr terms disappear in the thermodynamic limit L ^ oo 1161 . 
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and the same is valid for the other spin (|). The boson fields for the total charge and 
spin degrees of freedom are related to the formers as 



1 

1 

71 



[,/.t(x) + (^l(x)] 
[01 (x) -(plix)] . 



(2.18) 



and the 6 fields obey the same relations. From this expressions we can already see that 
interactions of the form (2.12), introduce term such as {V(j)u{x))'^ and {V6u{x))'^. 



Using representation (2.15 1, the Hamiltonian Hid can be bosonized obtaining a 
quadratic Hamiltonian of the form 



'Hp + Ti-cr 



dx 
2^ 



u,K,[Ve,{x)f + -^[VMx)f 
A,, 



(2.19) 
(2.20) 



where v = p{a) denotes the charge (spin) degrees of freedom, Uy is a velocity, 
Ky a dimensionless parameter depending on the interactions, and Oy and are the 
"new"bosonic fields. The parameters u and K are given by the interaction parameters 



Uy = VF[{l + yiy/2f -{yy/2f] 



211/2 



Ky 

2/v 



1 + + yy/2 

1 + yAy/2 - yy/2_ 

9i\\ - 5211 =F 92± 

54|| ± m 

gy/{TTVv), 



1/2 



(2.21) 



where the upper sign refers to p and the lower one to a. In the non-interacting case 
= vp and Kp = K^j = 1. For systems with spin rotation symmetry we have 



1. K, 



1 in the absence of interactions and Kp < 1 for 



5111 = ail. and thus 
repulsive interactions. 

The representation of a fermionic Hamiltonian in a quadratic form is one of the 
big advantages of the Luttinger liquid theory. It takes into account all the effects of 
fermionic interactions (momentum conserving interactions) and put them in a simple 



quadratic Hamiltonian of the form of Eq. (2.20 1. Although we will not treat the case 
of gi_\_ scattering here, it is important to know that this term cannot be written in the 



same quadratic form of Hamiltonian (2.191 because it gives a cosine term (called sine- 
Gordon Hamiltonian) and must be treated with renormalization group equations [19 J . 



The splitting of Hamiltonian (2. 19 1 into a charge part Tip and a spin part Tia reveals 



the complete separation of charge and spin degrees of freedom in ID systems, which 
is an important property of Luttinger liquids. It forbids single particle excitations (free 
electron) carrying spin and charge together. 

There are various properties of Luttinger liquids that we will not treat here since 
it would require more than an entire chapter, but there is one special property which 
is very important to know when one is working with one-dimensional systems: the 
power-law decay of correlation functions. In addition to the particle-hole excitations 
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mentioned at the beginning of this section there are also particle-particle excitations. 
These are collective excitations that represent charge-density and spin-density fluc- 



tuations respectively, with g ~ and q ~ 2kp components (see Fig. 2.6 1. When 
calculating charge-density or spin-density correlations functions in one dimension, the 
g ~ part gives a free fermion decay of the correlation and the q ^ 2kF (or multi- 
ples q ~ 2nk^) behaves as a non-universal power law, with an exponent depending on 
interactions. This arises because the charge and spin densities fluctuate in space and 
time. That is why ordered states do not exist in one-dimensional systems llT6l . Let 
us take for example the density-density correlation function of a system described by 
Hamiltonian ( |2.20| ), 

{p{x, r)p(O)) = + — ^ ^^^^ + PqA2 cos{2npox) (^-j 

/ a \ 

+P0A4, cos(47r/5ox) - + • • • (2.22) 



Where ya = ut + aSign(r) and r is an imaginary time related to t by the Wick 
rotation |2 |: t = it + eSign(t). po is the average density of particles. In Chapter [5] 
we will give the technical details necessary to compute a correlation function of this 



type. At this point we are interested in the physical meaning of Eq. (2.221. The g ~ 
part of the correlation decays as This is a Fermi liquid like decay and only 

the amplitude is renormalized by the interactions. The 2kj! and Akji parts, however, 
decay as a non-universal power law, with exponents depending on interactions. In 
this specific correlation, the 4A;f component does not depend on the spin part, thus if 
Kfj = 1 and Kp < 1/3, it becomes the dominant component. A system with this type 
of correlations is referred as a Luttinger liquid. Another important quantity showing 
a power-law decay is the retarded single-particle Green's function, which will be used 
repeated times in the next chapters. It is defined as |j2J 

G',%{x,t) = -ie{t){[i^,,{x,t),4,{0,0)]+), (2-23) 

where []+ is the anticommutator, 6{t) is the Heaviside-function, r = +1 (—1) for 
right (left) movers, and a denotes the spin. The retarded Green's function is usually 
obtained form the Green's function in imaginary time ||2l. The relation between these 



correlation functions will be explained in Sec. 3.1 In a Luttinger liquid, the retarded 
Green's function is given by 



= -i^e-'^Him — (x (2.24) 

ZTT e^o [ e + i{vFt — rx) 

j 1, + i{u,t - rx) V (a + iuut)^ + x^ ) ^ \t ^ -t 
with a the momentum cutoff. The exponent is 

j, = {K, + K~^ -2)/8>0 (2.25) 
For a spin rotation invariant system, = 1 and 7^ = 0. It is not our intention to 



explain how to obtain Eq. (2.24i because is out of the scope of this work. We refer 
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the reader to Ref. fT6l for a detail explanation on the calculation of the ID Green's 



function. We are interested in the asymptotic behavior of Eq (2.24 1, i.e., at x — > oo 



G-^(x^oo)~ (")'''\ (2.26) 



The behavior of correlation functions at large distances will be used constantly in the 
next chapters, to detect power-law behaviors in the calculations using scaling analy- 
sis. Because most of the functions in one dimension are power laws, most result can 
be obtain by simple scaling analysis. The single-particle Green's function in imagi- 



nary time Graix^T) follows the same asymptotic power-law behavior of Eq. (2.26 1. 
From Gra{x, r) we can obtain the occupation factor n{k), which is simply the Fourier 
transform of the equal time Green's function O 

nra{k) = j dxe-'^''Gra{x,T = 0-). (2.27) 

The occupation factor is thus the Fourier transform of a power-law and, for a spin 
rotational invariant system = 1, is given by lfT6l 

n(/fe)~|/fe-A:F|^/'[^''+^'^'l-^/'. (2.28) 

Instead of the discontinuity at kp characteristic of a Fermi liquid, in one-dimension 
one finds a power-law singularity. This means that the quasi-particle residue is Z = 0, 
which is an evidence that single-particle excitation cannot exist in ID. 

At this point we can summarize the features of Luttinger liquids that we have 
discussed until now. Firstly, we saw that in one dimension the dispersion relation can 
be linearized keeping the same low-energy properties of the system. This requires the 
definition of right- and left-going fermions. We learned that particle-hole excitations 
are well defined in ID and thus can be used as a basis to represent fermionic operators, 
thought a transformation known as bosonization. And we saw that a fundamental 
property of Luttinger liquid is the power-law decay of correlation functions. Now we 
will discussed the case of fermions moving on a lattice. For this we dedicate the next 
section to the one-dimensional Mott insulator and all the relevant properties for the 
study of transport physics in it. 

2.4.2 ID Mott insulator 

Let us begin by considering a system made of fermions moving on a one-dimensional 



lattice. This system is described with the Hubbard model studied in Sec. 2.1 We 
will apply the bosonization transformation to the Hubbard Hamiltonian of Eq. ( 2. 1 1 
in one dimension, where the sum over j disappears. The result is again a quadratic 
Hamiltonian of the form given in Eq. ( 2.19| ), but now parameters u and K are: 



UpKp 


= UaKa 


Up/Kp 


= vf(^1 


Ua/K„ 


= vf(i 



VF 

u 

TTVf 

u 

TTVf 



(2.29) 
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(a) 



Figure 2.8: (a) One dimensional Mott insulator: at half filling (one particle per site) 
due to the repulsion between particles, the system prefers to localize one particle on each 
lattice site. In this case, antiferromagnetism order is favored, (b) If the system is doped 
with holes or electrons, the particles (electrons or holes) can move and we thus have a 
metallic state (From Ref. 1 171 ). 



functions of the Coulomb interaction U. 

If there is one particle per site (half-filled system) and the repulsion U is larger 
than the kinetic energy t, the particles are localized on the lattice sites to minimize 
repulsion and the system is called a Mott insulator ||7]|. This is depicted in Fig. |2.8[ In 



the Mott insulating state a gap A appears in the charge excitation spectrum, but the spin 
properties are totally unaffected. To attain a metallic state, the system must be weakly 
doped away from half-filled case in order to have particles propagating. For other 
commensurate fillings, like quarter-filling (one particle every two sites), it is necessary 
to have a nearest neighbor repulsion (V) to form a Mott insulator lfT6l . The transition 
between the insulator and metallic state is known as Mott transition fS^ and is depicted 



in Fig. 2.9 It can occur changing the value of the repulsive interaction, in this case 



is called a Mott-U transition; or changing the doping 6 and is called Mott-(5 transition 



(see Fig. 2.9 1. Here we will not treat the case of doping and thus we refer to the Mott- 
U transition as just the Mott transition. The properties of Luttinger liquid discussed in 
the previous section remain essentially the same for the ID Hubbard model, knowing 



that interaction Ki, and velocity parameters, depend on U as shown in Eq. (2.29 1. 

To study transport properties in commensurate systems (as will be done in Chapter 
[3]l, we must take into account the effect of the underlaying lattice. Due to the presence 
of the lattice, the wavevector is defined modulo a vector of the reciprocal lattice, i.e., 
271 /a in one dimension with a the lattice spacing. As a consequence, new interac- 
tion processes will appear in which the total momentum is not conserved, such that 
^1 + ^2 + ^3 + ^4 = Q, where Q is a vector of the reciprocal lattice (in momentum 



conserving processes, such as (2.121, ki + k2 + k^, + = 0). In other words, the 
particles exchange momentum with the lattice. These processes are known as umklapp 
processes and they are the only ones responsible for dissipation in ID systems |[T6l . 
The umklapp processes appear only at commensurate fillings. For example in a half- 
filled system we have Akp = 2-K/a and it corresponds to a process where two electrons 
jump from one side to the other of the Fermi surface, transferring a momentum Akp to 
the lattice [il6l fTTIl . This is illustrated in Fig. 2.10 The umklapp term modifies only 



the charge part of Hamiltonian ( 2.19) , by adding a term of the form 



Wi/2n = gi/2n / c^x cos(n V80p(x) ) , (2.30) 
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Figure 2.9: Phase diagram of the one-dimensional Mott insulator, n denotes the order 
of the commensurability: n = 1 for 1/2-filled and rt = 2 for 1/4-filled systems. Int 
denotes the repulsive interaction between particles, /i is the chemical potential, S is the 
doping and A is the Mott gap. Ml and LL correspond to the Mott insulator and Luttinger 
liquid (metallic) phases, respectively. The critical exponent Kc and velocity Uc depend 
on whether it is a Mott-U or Mott-5 transition (From Ref. lfT6l ). 



(a) 





Figure 2.10: (a) Umklapp processes in high-dimensional systems appear regardless of 
the filling (provided \kp\ is large enough), (b) In one dimension an umklapp process ap- 
pears when two particles are scattered from one side to the other on the Fermi surface 
ikp = In I a, that is, for half filling. There are also umklapp processes for other com- 
mensurate fillings (From Ref. [ 17j). 
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where n is the order of the commensurabihty (n = 1 for half-filling and n = 2 for 
quarter-filling), gi/2n '^^e coupling constant, and the commensurability n is related 
to the reciprocal lattice vectors by 2pkF = 2TTq/a, where p = 2n and q are integers. 
In the case of one particle per site (n = 1, half-filled band), we rewrite Eq. (2.30 1 as 

m 



n 



2g3 



dx cos{V8<f>p{x)). 



(2.31) 



(2™)2 

We have replaced the amplitude g^^^g by the factor 2g^/{2TTa)'^, where the coefficient 
gg is of the order of the Coulomb interaction U and a is the lattice spacing. The 



umklapp process of Eq. (2.31 1 will be the only term, producing dissipation, considered 
in our study of the Hall effect in a quasi one-dimensional system. 

After reviewing the necessary tools to understand the transport properties of fermions 
moving on a one-dimensional lattice, we will describe an extension of the pure ID case 
which is the lowest dimensional situation where the Hall effect has a sense: quasi one- 
dimensional systems. This will be the type of strongly correlated low-dimensional 
system treated on Chapter |5] 



2.4.3 Quasi one-dimensional case 

In the group of low-dimensional systems we find, beside the purely one-dimensional 
ones, systems with a dimensionahty between one and two or one and three. These sys- 
tems are known as quasi one-dimensional and they consist of one-dimensional chains 
coupled by an interchain coupling term of the form 



-Hi. 



dx t_L 



(2.32) 



with (/i, v) a pair of chains and t^^^^y, the hopping integral between these two chains. 
The simplest example of a quasi ID system is composed of two one-dimensional 
chains coupled to form a fermionic ladder [16|. Here we are interested in systems 
made of a large or infinite number of coupled chains. In Chapter [5j we will use 
Eq. (2.32 1 together with Hamiltonian ( 2. 14| ) to describe our system of weakly coupled 
one-dimensional chains. 



In addition to the coupling term (2.32 1, there exist hopping processes of second 
order in t^, where two particles jump between chains. These terms have a well de- 



fined classical limit (unlike the coupling term in Eq. (2.32 1) and for a sufficiently large 
number of chains they can be treated using a mean field approximation ||20| . Their 
most important effect is to drive the system into an ordered state 1 17 1. The ratio tx/iy, 
where t\\ is the intrachain hopping, determines the effective dimensionality of the sys- 
tem. When the system is nearly isotropic, t±_ ~ ty, we deal with a high-dimensional 
situation where even small interchain interactions must be taken into account because 
they lead to a Fermi liquid state or another correlated state. In the opposite limit 
t_L <C the system is highly anisotropic and the chains are in a well defined Lut- 
tinger liquid regime. Thus, the processes of second order in t±_ are less likely to occurs 
and can be neglected. This is the limit that interest us for the discussion of trans- 
port in quasi one-dimensional systems. We will therefore focus only on the interchain 



hopping term of Eq. (2.32 1, and we refer the reader to the literature for an extensive 
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Figure 2.11: Open Fermi surface of a quasi one-dimensional system, (a) If the tempera- 
ture (or any other energy scale) is larger than the warping of the Fermi surface produced by 
the interchain hopping, the system cannot feel it and thus behaves as a one-dimensional 
one. (b) For temperature/energy smaller than the warping, the system feel the two- or 
three-dimensional nature of the dispersion and thus behaves as a high-dimensional sys- 
tem (From Ref. ifTTl). 



discussion of the effects of interchain interaction terms |[T6l [TTl . Therefore, in the rest 
of this manuscript the word "interaction"wilI only refer to intrachain interaction. 

One important effect of the coupling term ( 2.32| ) is to induce a dimensional crossover, 



as a function of decreasing temperature (or another energy scale), from a one-dimensional 
regime to a high-dimensional one (2D or 3D). To understand this dimensional crossover 
we consider the dispersion relation of the system in the absence of interactions: 

e(A;||, k±) = — 2t|| cos(fc||a) — 2t± cos{k±b) (2.33) 

where b is the distance between chains. In the limit t± <^ the Fermi surface given 
by Eq. ( |2.33 1 is an open surface as the one shown in Fig. 2.1 1 If the temperature (or 



energy) is larger than the warping of the Fermi surface due to interchain hopping, the 



system cannot be sensitive to it (see Fig. 2.11 1 and thus feels a flat Fermi surface cor- 



responding to a one-dimensional regime. In contrast, when the temperature (energy) 
is much smaller than the warping, the system behaves as a two- or three-dimensional 
one. This dimensional crossover occurs, in the non-interacting case, at an energy scale 
of the order of the interchain hopping t±. Interactions renormalize this energy scale to 
(unities are taken in order to have /c^ = 1) 

E* = T,ir^W i^j , (2.34) 

where ( = J2u ^'^'^ = (-^f + ~ 2)/8 > For a spin rotation invariant 
system = 1 and 7^ = 0. Furthermore, in the non-interacting case C = and we 
recover T^i ~ t±. For interacting systems we have C > and thus the energy scale at 
which the dimensional crossover takes place is reduced, making the system effectively 
more one-dimensional. This can be understood in the following manner: the interchain 
coupling involves single-particle hopping processes which are unstable excitations in 
a Luttinger liquid, then in order to have an electron jumping from one chain to the 
other a collective excitation must brake and then recombine in the new chain. This 
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Figure 2.12: Deconfinement transition: quantum phase transition that takes place at 
r = as a function of the interchain hopping t^. At finite temperatures, one can have 
a crossover between a Luttinger liquid (LL) and a Mott insulator (Ml) if t± < t*j_ by 
lowering the temperature (arrow 1), or between a LL and a high metallic phase (HDM) if 
tj_ > t*j_ (arrow 2). A is the Mott gap (From Ref. fTTl). 



makes the single-particle hopping very difficult. As mentioned previously, there exist 
hopping processes of second order in t± which are more favorable, where two particles 
jump between chains. These processes lead to an ordered state at a temperature Tx2 
which depends on the precise coupling (spin-spin, Josephson term or density-density) 
|[T6.I . In the limit t± <C the dimensional crossover always occurs first (Txi > Tx2) 
because two-particle processes are of the order <^ 1 and thus are less likely to 
occur. 

Another property of quasi ID systems that we shall mention in this section is the 
deconfinement transition. It consists of a quantum phase transition where the system 
goes from a one-dimensional Mott insulator to a high dimensional metal with increas- 
ing interchain hopping 1 16l. It appears in commensurate systems which are generally, 
as we saw in the previous section, Mott insulators with a gap A in the charge excita- 
tions. A qualitative picture of the deconfinement transition is given in Fig. 2.12 This 
transition is more complex than the dimensional crossover, because when it occurs 
electrons are able at the same time, to leave the chains and to conduct. The critical 
value of t_i_ at which the deconfinement occurs is t\. One should solve the fulled cou- 
pled problem to obtain a critical value for t\. It is quite difficult to extract physical 
properties in the deconfined phase. One known quantity is the transverse conductivity 
that will be discussed in the next chapter. 

At this stage, we have seen the main properties of low-dimensional systems related 
to transport physics. We know how to describe a strongly correlated quasi ID system, 
made of weakly coupled Luttinger liquids. In addition, we can introduce the effects 
of the lattice through the umklapp scattering review in this chapter. All these will be 
implemented in Chapter |5] 
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CHAPTER 3 



Transport in strongly correlated systems 



In order to study the Hall effect in strongly correlated systems, which is above all a 
transport measurement, we must first investigate the various existing methods to treat 
transport in systems with strong interactions. In this chapter we make a review of two 
related methods and then use them to obtain the conductivity in a one-dimensional 
system, with the purpose of clarifying their domain of applicability. 

A discussion about transport in correlated systems cannot start without mentioning 
the most simple theory used to study transport in nearly-free electron systems. This is 
the Drude theory of metals |[T2l (Drude constructed his theory on electrical and thermal 
conduction in 1900). This theory is the precursor of all modern transport theories and 
is still widely used. The Drude model assumes electrons are classical particles that 
can move freely and experience collisions with other electrons with a probability per 
unit time 1/r. The time r is known as the relaxation time. In accordance with Ohm's 
law there is a relation between the potential drop V and the current / flowing along a 
wire, which is F = IR, where R is the resistance of the wire. This relation can be 
rewritten in terms of the applied electric field E and the current density j as E = pj, 
where p is the resistivity of the wire. If n electrons of charge — |e| per unit volume 
move with an average velocity v, then the current density is proportional to v and is 
given by j = —\e\nv. After each collision, the momentum of the electrons is changed 
by Ap = FAt, where F = —\e\E and the average time At between collisions is r. 
Then, the average velocity of the electrons is given by v^^g = —\e\ET/m, where m is 
the electron mass. Thus the current density can be rewritten as 

j = ^^E. (3.1) 

m 

This formula is usually given in terms of the conductivity a = 1/p. 

I |2 

n e T 

j = aE; a = ^^. (3.2) 

m 

The formula on the right is the known as the Drude formula and gives an estimate of 
the dc electrical conductivity in terms of known quantities, except for the relaxation 
time r which is generally determined experimentally [12,1 . 
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In the case of an applied time-dependent Electric field, E{t) = Re(£J(a;)e~*'^*), 
the current induced in a metal j{Lo) is given by: j{u}) = a{io)E{uj), where (T(a;) is 
known as the frequency dependent (or AC) conductivity. In the Drude theory of Metal 
(t{lo) is equal to |[T2t 

= ■■ — > o'o = ■ (3.3) 

1 — ILOT m 



This result reduces to the DC result of Eq. ( |3.2| ) at w = 0. In the next sections we 
describe two quantum theories to compute the conductivity a. First, we discuss the 
linear response theory used to compute ac conductivities {a{uj)), and then we explain 
the memory function formalism that allows to compute cf{u}) including the effect of 
interactions. Finally we applied all these methods to the study of transport in low- 
dimensional systems. 



3.1 Linear response theory and Kubo formulas 

When an external field is applied to a physical system, its response will depend on the 
magnitude of the perturbation. Many experiments in condensed matter physics mea- 
sure responses to applied fields. For example, a sample is placed in an electric field, a 
magnetic field or a temperature gradient to measure its electric current, magnetization 
or thermal current, respectively. In linear response theory, the applied field is assumed 
sufficiently small, so that the system response increases linearly with the intensity of 
the field, and non-linear terms can be neglected. 

Although linear response theory applies to all possible fields Q (Kubo first derived 
his equations for the electrical conductivity in solids in 1957-1959), we will focus here 
on the electrical conductivity tensor cj(r, r', t — t') which relates the current response 
ja{f, t) to an electric field Ep{r' ,t') trough 



jair,t) 



[ dr' f dt'aap{r-r';t-t')Efs{r',t'), (3.4) 



where Ep{r', t') is the total electric field, i.e the sum of the applied external field E 



i(ext) 



and the fields created by the charged displaced in the solid. In Eq. (3.4i we write the 
conductivity as a function of the difference (r — r') because we assume translational 
invariance in the solid. 

The main goal is to obtain an expression for the ac conductivity (Jq, ,g(q, c<j), that 
is the Fourier transform of a a/sir — r' ,t — t'). Let us consider a system described 
by a Hamiltonian Ti, where an electric field is applied. For this formalism to be well 
defined, the system is considered unperturbed at t = —oo and the perturbation, in 
this case the electric field, is turned on slowly to be totally present at t 0. Then, 
at large times, t = +oo, the perturbation is turned off and the system returns to its 
unperturbed state. A(r, t) is the time dependent vector potential related to the electric 
field by E{r, t) = — 5A(r, t)/dt. The current is defined as the functional derivative 
of the Hamiltonian with respect to the vector potential (unities are chosen in order to 
have c = 1, h = 1 and ks = 1) 

Ja[r,t) -- 



6Aa{r,t) 



(3.5) 

A'='=0 
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and it can be itself proportional to A. We also define the operator of the total current 
as Jq = / drja{r). The vector potential also enters in the Hamiltonian modifying 
the momentum 11 of the particles by 11 — \e\A, with |e| the elementary charge. Thus, 
to capture all the linear response to the vector potential, one should make a first order 
expansion in both the Hamiltonian and the current. In the case of the current, the 
expansion of the average value is given by 



{ja{r,t)) = Y, / dr'dt 



P 



SAf3{r>,t> 



(3.6) 



with (• • • ) the thermodynamic average taken with respect to the Hamiltonian H. Here 
we suppose that TL can be diagonalized in some basis for the average (• • • ) to be 
computed analytically. In the same manner, the total Hamiltonian is expanded to first 
order in A giving 



n = n[A = 0]- I drJ2j>,t)Aa{r,t) 



(3.7) 



where j°(r, t) is the part of the current independent of A. Using the above result in 
Eq. (|3.6[) we obtain the following expression for the average current 



0„(r,t)) = ^ / dr'dt' 



/3 



=0] 



Ap{r',t'). 



(3.8) 

{j^{r,t); j^{r' ,t')) ret stands for the retarded current-cuiTcnt correlation function ID 
defined as 

{ja{r,ty,j^{r' ,t'))ret = -ie{t - t') {[ja{r,t),jf,{r' ,t')]) . (3.9) 
In Eq. (|3.9|) we dropped the superscript to lighten notation. Using the definition of 



the current given in Eq. ( |3.5| ), we rewrite the first term in Eq. ( |3.8| ) and we obtain for 
the average current 



{ja{r,t)) 



J dr'dt' 



5^n 



6A^{r,t)6Af^{r',t') / ^ 



-{ja{r,t);jf3{r',t')U] Afs{r',t'). 



(3.10) 



As pointed out before, the Hamiltonian is a function of 11— | e| A and thus the functional 
derivative with respect to A can be written in terms of the momentum operator 11, 



A=0 



6{r - r')6it - t')5^p. 



5A^{r,t)5Ap{r',t') 
The average current is finally given by the following expression 



(3.11) 



Oa(r-,t)) 



J dr'dt' -eH^p(^-^y{r-r')5{t-t') 



-{3a{r,t)-jp{r',t'))r,;\ Ap{r',t'). 



(3.12) 



25 



Making the Fourier transform 

Oa(Q,'^)) = jdrj (ite-^''-V-*(j„(r,t)), (3.13) 

and using the relation between the vector potential and the electric field in Fourier 
space: Aa{q,io) = Eci{q,io) /iuj, we finally encounter an expression for the con- 
ductivity matrix (Tap describing the linear-response current ja{<l, induced by an ac 
electric field Ea{q, uj) 

ja{q, = ^ CTapiq, ^)E(3{q, w), (3.14) 

where cia^ is 

0-Q/3(q, = — [XQ/3(0)(5a/3 - Xa/3(q, ^)] ■ (3.15) 



Eq. ( |3.15 1 is known as the Kubo formula. The first term is called the diamagnetic term 



and is purely imaginary 

2 / \ _ / '^^W 



(3.16) 

A=0 



The second term is the retarded current-current correlation function in Fourier space 

poo 

Xa/3{q,io) = -i dte'^'-'-''^e{t-t'){[Uq,t)Jfs{q,t')]). (3.17) 



The real part of the conductivity is totally given by this term. The dc conductivity is 
found by taking the limit q —>^ first and then the limit tj ^ 0. If the order of these 
limits is interchanged, one would obtain the thermodynamic response of the system, 
where the limit lj ^ is taken first in order to have a static perturbation with g 7^ 0. 

The retarded current-current correlation function is usually computed in imaginary 
time T using the Matsubara formalism [2J, because it is the standard way to perform 
the calculations at nonzero temperatures. For this we first define the current-current 
correlation function in imaginary time 

Xap{q,r) = -{Trjl{q,T)jp{q,0)), (3.18) 

where < r < /? and /3 is the inverse temperature (/? = 1/T). Tr is the r-ordering 
operator, which arranges operators with earliest r to the right |^. Then, we perform 



the Fourier transform of Eq. ( 3. 18 1: 



Xaf3iq,iuJn) = dre'^^-^Xapiq^T). (3.19) 

The frequencies iuj^ are called imaginary Matsubara frequencies and have the fol- 
lowing values: for bosons (j„ = Inir/ (i and for fermions ujn = (2n + 1)tt/ (i. To 
recover the retarded correlation function from the Matsubara function we have to 
change the imaginary frequencies to real frequencies making the analytical continu- 
ation: iiOn + 16, 

Xaf3{q,i<^n) Xafs{q,Uj): (3.20) 
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with 5 = 0^. Finally by inserting the result in the Kubo formula (3.15 1 we obtain the 
conductivity. 

Until now, we have seen how the conductivity can be expressed in terms of a 
current-current correlation function plus a diamagnetic term, via the Kubo formula 
(3.15 1. This formalism is suitable for systems with a Hamiltonian that can be diago- 
nalized, i.e., in which the averages (• • • ) can be computed analytically. However, if 
we are interested in strongly correlated systems, as is the case in this work, we must 
consider interacting terms in the Hamiltonian that cannot be diagonalized. Thus we 
must search for other approaches where the conductivity can be calculated with in- 
teractions included. The Kubo formula is, however, the most general expression to 
compute ac conductivities and serves as a basis for the other approaches. The next 
section is devoted to one of this approaches known as the memory matrix formalism. 



3.2 Memory function formalism 

To understand the essence of this formalism, we first need to study some important 
properties of correlation functions. The current-current correlation function defined 



in Eqs. (3.17l-(3.19l, belongs to a larger group of functions known as susceptibilities 
denoted as xab{^) = B) at q = 0. They are holomorphic functions for all 
complex frequencies lu. Integrating by parts Eq. ( 3.17| ) we can extract the behavior at 
large frequencies (this is done in detail in Sec. |4.2|) 



{[A,B]) , {[[A,n],B]) 
Xab{uj) = \ y for Lo ^ oo, (3.21) 

CO 

with [A, B] the commutator between operators A and B at the equal time, and H the 



total Hamiltonian of the system. If A and B commute, it is the second term in (3.21 1 
which gives the behavior at high frequency. In the imaginary-time representation of 
the susceptibilities xab{t) = —{TrA{T)B{0)), the operators obey the Heisenberg 
time evolution A{t) = e^'^Ae"^'^, which implies the following equation of motion 
9r^(T) = [H,A{t)]. With this, one can demonstrate that susceptibilities obey the 
following equation of motion 

uj{A; B) = {[A, B]) - {[H, A]-B) = {[A, B]) + {A; [H, B]). (3.22) 

In order to prove this, we have taken the time derivative of xab{t), applied the time 
homogeneity property {TrA{T)B{{))) = (Tt-A(0)-B(— r)) and then Fourier trans- 
formed the whole equation. We recall that the retarded correlation function x{'^) is 
obtained from the Matsubara function xi^^n) by the analytical continuation iun 

uj + i6. 

In general, susceptibilities can be represented as a spectral integral (also called 
Lehmann representation (1954)) of the form 

Xi^^n) = / 7, (3.23) 

where the spectral function S{uj) is given by the imaginary part of the retarded corre- 
lation function defined in Eq. ( |3.17| ) 

S{uj) = -2\m[x{u:)] and ± i<5) = Re[xH] ± ilm[x(c^)]. (3.24) 
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One can also demonstrate that x{^) verifies tlie following symmetry properties 



X{^^) = X{-^) and x*(w) = x(^^*)- 



(3.25) 



Thus, Re[x(w)] andlm[x(w)] are real and satisfy Re =Re[x(— w)] andlm[x(tij)] 
— Im[x(— w)]. All these susceptibility properties will help us define the memory func- 
tion properly. 

Now let us come back to the conductivity a. We suppose a system described by a 
Hamiltonian H = Hq + Hint, where T^int is an interaction term that we do not know 
how to treat exactly. We will start with the case where the conductivity goes only 
along one direction and thus the Kubo formula (3.15 1 becomes aapiio) = (7{uj) = 
[x(0) — x{^)] /^i^- As pointed out before, the static conductivity, in a normal conduc- 
tor, is given by the diamagnetic term x(0) (w = term in Eq. (3.15 1) and for all uj we 
have: 

XH/X(0). (3.26) 

The memory function formalism consist on representing the conductivity by a relax- 
ation or memory function M{u)) 1211 . as mentioned before. Let us consider the func- 
tion 

x(o) - x(w) ■ 



iM{uj) 



(3.27) 



Due to inequality (3.26 1 the memory function is also an holomorphic function for all 



frequencies uj. The behavior of x{^) at large frequencies shown in Eq. (3.21 1 gives 
an asymptotically decrease of M{uj) as l/w^ when — > oo (because the commutator 
[J, J] is zero). The fact that M{uj) vanishes at infinite frequency will be determinant in 
the study of the high-frequency Hall effect in Chapter |4] The susceptibility properties 
shown in Eq. ( |3.25| ) imply the following symmetry relations for the memory function: 
AI*{io) = M{u!*) and M{uj) = —M{—u>). Furthermore, it can also be represented by 
an spectral integral 

, fduj'S{io') 



2tt uj' 



UJ 



(3.28) 



where again 
S{uj') = -2Im[M(u;)] and 



M{uj ± i6) = Re[M(a;)] ± ilm[M{uj)]. (3.29) 

Thus, Re[M(a;)] andIm[M(u;)] are real functions satisfying Re[M(u;)] = -Re[M(-cj)] 
and Im[M(a;)] = Im[M(— w)]. From Eq. (3.27) we can rewrite the susceptibility in 
terms of M {uj), 

UJ + iM[uJ) 

Now that we have represented the current-current correlation function in terms of the 
memory function, it can replaced in the Kubo formula to derived an expression for the 
conductivity in terms of M{uj), 



a{uj) 



. x(o) 

UJ + iM{uj) 



(3.31) 



This representation of the conductivity provides a correct way to make perturbative 
expansions in small parameters (like a coupling constant or density), that are not feasi- 
ble on susceptibilities due to their singular character at small frequencies. Comparing 
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Eq. (3.31 1 to the semiclassical fomiula (3.3 1, derived from the Drude model, it is ob- 
vious that the memory function plays the role of the relaxation time r in the Drude 
conductivity. The advantage of the memory function formalism is that iM{u) gives 
a practical way to obtain the conductivity which is expected to be non-singular when 
a; — > 0. Eq. (3.31 1 will be frequently used in this work to compute conductivity. 

The next step is to see how the memory function is computed by treating the inter- 
action term "Hint perturbatively. From the equation of motion ( 3.22[ ) one can write 

io{J;J) = -{K;J) (3.32) 

where K = [H, J] is an operator known as residual force. It is given by the part of 
the Hamiltonian that does not commute with the current, i.e., the interaction term 
(K = [Hint, J])- Furthermore, using the same equations of motion for operators K 
and J, we have 

uj{K;J) = {[K,J]) + {K;K). (3.33) 
And from Eqs. (3.32i and ( 3.33[ ) at w = we obtain 

{[K,J]) = -{K;K)^=o. (3.34) 

Using Eqs. (3.32 i-(3.34l, and the current-current correlator defined in (3.17 1 at g = 
(where operator j becomes a total current operator J) we find 

{K; K) - {K; K)^=o 



UJ 



(3.35) 



Until now we have worked with exact relations. However, we do not know how to com- 
pute the correlator {K; K) with interaction terms present in the Hamiltonian. Thus, 
we have to make some approximations. First, we expand iM{uj) in Eq. (3.27 1 at high 
enough frequencies where |x(i^)/x(0)| is small 

x(0)^M(c^) -c^xH. (3.36) 

With this approximation we arrive to the following expression for the memory function 

1 {K-K)^-{K-K)^=o 



iM{uj) ~ 



x(o) 



(3.37) 



Because the operator K is proportional to the interacting term Hint, the thermodynam- 
ical average in correlation {K; K) can be computed with the Hamiltonian Hq (denoted 
(• • • )°) to get the lowest order in the interaction parameters, 

1 {K-K)l-{K-K)l^, 



iM{u) 



x(o) 



(3.38) 



Finally, calculating iM{ijj) and replacing it on Eq. (3.31 1, we obtain a result for the 



ac conductivity at second order in the interaction parameters (each K contributes with 
one parameter), which is the lowest non-zero order for the memory function. The 
application of these results will be much more clear in Sec. |3.3[ where it will be used 
to obtain the conductivity in a one-dimensional system. 

In this entire discussion we have supposed a system with longitudinal conductivity, 
where the memory matrix reduces to a scalar function. However, if one is interested 
in the conductivity in a plane (a more general case), the previous results must be reob- 
tained using relations involving matrices. In the next section we will study the memory 
function formalism using the conductivity tensor in order to obtain a matrix expression 
for M(lj). 
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3.2.1 Matrix representation of the Memory function 



In the study of the Hall effect, one is confronted with a problem of conduction along 
different directions in the presence of an applied magnetic field, as will be explained 
in Chapter |4] In this case, the previous derivation for the memory function must be 
remade starting from the conductivity tensor, in order to get a matrix representation 
for M{lo). The conductivity tensor is given by 

o- = f """" '^"^^ . (3.39) 

It is important to observe, for the following derivations, that the non-diagonal terms of 
a appear when a magnetic field is applied perpendicular to the x-y plane. Thus, for 
the rest of this section we will suppose that a magnetic field is indeed applied. Let us 
rewrite the conductivity tensor in terms of the memory matrix M{u}) as 

icr(tj) = x(0)[tjl + n + iM(a;)]"^ (3.40) 



Eq. ( 3.40| ) as well as the following definitions, are enclosed in a general theory called 



the Mori theory, and details can be found in Ref. Il22ll . As mentioned previously, the 
advantage provided by the memory matrix formalism is the possibility to make finite- 
order perturbation expansions which are singular in the conductivities due to their 



resonance structure 121]. x(0) in Eq. (3.40 1 is a diagonal matrix composed of the 



diamagnetic susceptibilities in each direction. 



The matrix in (3.40 1 is called frequency matrix and is defined in terms of the 



equal-time current-current correlator as ll22l 

^,u = ^^{[J,,JA)- (3.42) 

The frequency matrix gives the behavior of cr{oj) at high frequencies because the mem- 
ory matrix vanishes as l/o;^ when uj — > oo, as mentioned in the previous section. From 



Eq. (3.401 one can directly express the memory matrix M{ijj) in terms of the conduc- 
tivity tensor, 

iM{u) = -ia-^{u:)x{0)-u:l-n. (3.43) 



The diagonal terms Mxxi'-j-^) and Myy{uj), are given by definition (3.27i with x{^) 
replaced by Xxx{^) and Xyyi'^)' respectively (idem for x(0)). For the off-diagonal 
terms we have Myx{io) = —Mxy{uj) (due to Gxy = —(^yx), and we take Mxy{uj) 
because it will be the term necessary in the description of the Hall effect in Chapter]?] 
Mxy{u>) written in terms of the conductivities gives 

^Mxy{u.) = , YTt!:'^ f ^ - ^3.44) 
axx{uJ)cryy{uJ) + criyito) 

Expressing the conductivities in terms of current susceptibilities by means of the Kubo 
formula a^,^ = ^[Xfi.v ~ ^fivXiii^)]^ the above expression leads to 



^Xy{^)Xxy{i^) 
[Xx(0) - Xxx{^)] [Xy{^) - Xyy{^)] 



'Mxyiu;) = -A.v-^ A..v^; _ ^^^^ ^3_45^ 
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This representation, together with Eq. ( |3.27 1 for the longitudinal terms, completely de 



fines the Memory matrix. The parity property in Eq. ( 3.25| ) is valid for the longitudinal 

Xxx{^) and Xyy{^)^ giving for the diagonal term M^xi^j) = —Mxx{—<jj), and the 
same for Myy{uj). In the other hand, the transverse susceptibility satisfies Xxy{^) = 
~Xxy{—^) and from Eq. (3.45 1 it is easy to prove that Mxy{uj) = Mxy{—u)). 



Following the same procedure applied in the longitudinal case, we rewrite Eq. ( 3.45 1 
at high enough frequencies, such that |x^/^(w)/Xm(0)I small. In this expansion we 
use the equation of motion for the susceptibilities, as well as the relation \Hq, J^] = 
—^vfj,Jv, with = x,y and summation over repeated indices is implied. This rela- 
tion is the precondition necessary to obtain a regular expansion of the memory function 
for all frequencies to leading order in the interaction term [22] . The obtention of the 



following result will be presented in detail in Sec. 4.3 The off-diagonal term of the 
memory matrix thus reads 

iMxy[u:)^ 1 (K^^^ ^^^^^ 

where are the residual forces operators defined in the previous section. In this 
case, they are given by the part of the Hamiltonian which in the absence of magnetic 
field does not commute with the currents, i.e. = ["Hint, J^i]- When a magnetic 
field B is applied, the non-interacting Hamiltonian does not commute anymore with 
the currents, that is why it must be emphasized that K operators are computed with 
S = 0. The quantity {K^] Ky) stands for the retarded correlation function of the 
operators K^. The obtention of this correlator will be the key point in the study of the 
Hall effect in a quasi ID system in Chapter [5] 

Expression ( |3.46 1 does not contain the term at = present in Eq. (3.37 1, be- 



cause the off-diagonal terms of the memory matrix are even in frequency (therefore 
ujiMxy{co) is odd). Only the diagonal terms ai^e odd in cu and thus have the residual 



forces correlator evaluated at zero frequency. The terms omitted in Eq. (3.46 1 are ei- 
ther of second order in |x/x/t(i^)/XAt(0)l' or of second order in the magnetic field. As 
will be seen in the next chapter, Eq. p.46| ) is the memory matrix element necessary to 
compute the Hall resistivity. 

Now that we have a technique to compute longitudinal and transverse conductivi- 
ties in strongly correlated systems, to leading order in the interaction term, we will see 
in the next sections one example of their application in low-dimensional systems. 



3.3 Transport in low-dimensional systems 

After reviewing the main properties of low-dimensional systems in Chapter [2] we will 
devote this section to the study of "transport" properties in these systems. For this 
we make use of the two formalisms discussed previously. We will begin studying 
transport in purely one-dimensional systems and then in quasi one-dimensional ones. 
We concentrate on systems with commensurate fillings because in Chapter [5] we will 
investigate the Hall effect on quasi ID commensurate materials. 
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3.3.1 ID systems without umklapp scattering 



Transport properties are commonly used as a probe for Luttinger liquid (LL) behavior 
in low-dimensional systems. The measured conductivity or resistivity helps determin- 



ing the energy scales at which the sample is in a Luttinger liquid regime. In Sec. 2.4. 1 



we saw that the umklapp scattering is the only interacting term producing dissipation 
in a Luttinger liquid. Thus, if one removes this term for a moment and thus momen- 
tum is conserved in all interaction processes, the system should behave as a perfect 
conductor. In order to prove this we will compute the conductivity in a LL without 



umklapp scattering, making use of the Kubo formula presented in Sec. 3.1 First, we 



calculate the current using the charge part of the ID Hamiltonian (2.19 1 (because the 
transport properties do not affect the spin degrees of freedom) in the form (again we 
put h = 1 and c = 1) 



n. 



^ ^n,K,[vrn,(x)]^ + -g[V<Ap(x)]^ 



(3.47) 



where n(x, t) = {1/'k)V9{x, t). The vector potential enters only in the 11 part of the 
Hamiltonian via the substitution Ilp{x, t) —>■ Ilp{x, t) — eA{x, t) /tt. Then, the current 



operator defined in Eq. ( 3.5 1 is given by 

j{x,t) = e{V2upKp)Ilp{x,t), 



(3.48) 



where the factor \/2 comes from the sum over spins (Hp = (11^ + Ilj)/ lfT6l ). In 
the same way we calculate the diamagnetic term x(0)» defined in (3. 16 1, which for 
fermions with spins gives 



X(0) 



2e'^UpKp 

TT 



(3.49) 



As we saw in Sec. 3.1 the conductivity is equal to the diamagnetic term plus the 
current-current correlation term. Using ( 3.48| ) the current-current correlation function 



in imaginary time (defined in Eq. ( 3. 18 1) thus result 

X{x, x'; t-t') = - {euKf (T,n(a;, r)n(x', r')). 



(3.50) 



We remove the subscript p in order to lighten the notation, remembering that all these 
transport calculations affect only the charge part of the ID Hamiltonian. At this point 



we can write down the Kubo formula (3.15 1 in bosonization language [1231: 



a{uj) 



e^2uK 



TT 



(3.51) 



where x{^) is the Fourier transform of Eq. (3.50 1. Note that 8e uK plays the role of 



the plasma frequency in the usual formulas for the conductivity fTl\ . 

As we need to evaluate time-ordered correlation functions x{^) order to obtain 
the conductivity, we give below the necessary tools to compute correlation functions 
using functional integral techniques. Some important results of functional integration 
are presented here but we refer the reader to the literature for a complete review L24,l . 
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The partition function Z = Tr (e '^^) of Hamiltonian (3.47 1 represented via a func- 
tional integral is ll24l 



Z = j V(j){x,T)VU{x,T) 



(3.52) 



V</'(x,t) -H(0(x,r),n(x,r))]. 



S is the action in imaginary time associated with Hamiltonian (3.47i. Time-ordered 
correlation functions of two operators A and B, functions of the operators (/> and 11, 
are defined in the functional integration formalism by the formula 

1 



(r,A(x,T)s(o,o)) 



z 



r)Pn(x, t)A{<P, Ii)^^,B{<p, n)o,o e-^ , (3.53) 



where Z and S are given in Eq. (3.52 1. From now on we simply denote (Tr...) as 
(...). An important advantage of this technique, is that A and B on the right hand 
side of Eq. (3.53 1 are the value of operators A and B, respectively, and thus have 



the properties of scalar fields, that are much easier to deal with than operators 11241 . 
Another useful formula from functional integration (which we will not prove here) is 
the one corresponding to the Fourier transform of the correlation function (n(r)n(r')) 

/ Vu[q]u*{qi)u{q2)e-'2^^ A{g)u' (g)u{g) 



{u*{qi)u{q2)) 



1 



Vulqle 



4E„A(gK(g)«(<?) 



Aiqi)6, 



(3.54) 



where A{q) is a diagonal matrix and, for real fields n(r), we have u*{q) = u{—q). 
Now we are ready to compute the correlation function in Eq. ( 3.50| ). 

(n(x, T)n(x', r')) = ^ y ^'/'(^' r)Pn(x, T)n(x, r)n(x', r ) e"^. (3.55) 

To express the Fourier transforms we use the following notation: r = (x,nr), q = 
{k,ujn/u), and e**?'^ = e*^'^^-'^"^). Then n(r) = Egn(q)e^''-^, with /3 = l/T 
and i7 the volume of the system. Thus, the above expression becomes 



(n(n)n(r2)) = ^ ^ (n*(qi)n(q2)>e-*'^^^^e 



qii<j2 



(n*(qi)n(q2)) 



V<Piq)VU{q)U*{q,)U{q2)e- 



(3.56) 
(3.57) 



The action given in Eq. ( |3.52| ) must be also written in Fourier space, 

-y 



s 



u:n<t>{qm-q) " ^n(q)n(-g) 



2ttK 



k'<p{q)H-q) 



(3.58) 

Next we complete the squares on the 11 part of the action to get ride of the linear term. 
The action thus result 

S = -^"V^^'H<l)<H-q)+i^y'^e4>{q)4,{-q) (3-59) 



q 

n(q) 



" (3n^ 2TrK 

Q 



(5VL 



UJ-n 



uKtt 



[n(-Q) 



+ 



OJr, 
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Making the change of variable n(q) = n(g) — :;^^4'{fl), we obtain an action com- 
pletely separable in Sff, and S-jj part, 



1 /I 



E 



{uKTr)fl{q)U{-q) + 
/ uKir 



+ 



uKtt ttK 







_u^,u^j (nq, 0q) 
uKtt Ktt ^ 



(3.60) 



The matrix representation defines the diagonal matrix ^(q) to be inserted in formula 
(3.54 1. Using the same procedure, the partition function can also be separated in Z = 
Zfj,Z]j- We rewrite the correlator (3.57 1 in the new variables and simplify it using 
(1 /Z^) f Vcf) e-^'t- = 1 (the same is valid for ft). 



(n*(9i)n(g2)) 



1 



z 



n 



Pn(q)n*(qi)n(q2)e-^ft 



{uKtt)' 



(3.61) 



At this point, we can make use of formula ( |3.54 i. It is evident that the cr ossed corre- 
lations (n</>*) and (0n*) vanish because A is a diagonal matrix (see Eq. (3.60l). The 
(n*n) correlation function thus gives 

,2 



(n*(qi)n(q2)) 



1 



2 



(3.62) 



Moving to real space and multiplying the result (3.62i by (euK) we obtain for the 
current-current correlation function of Eq. p.50 1 

ieuKf{U{x, t)U{x', t')) = ~ {dr(p{x, r)a,0(x', r')) + ^—6{x - x')5{t - r') . 

vr^ vr 

(3.63) 

After summing over spins the above expression (a factor \/2 appears for each 11 and (f)), 
the second term on the right cancels exactly the diamagnetic term given in (3.49 1, when 
using Eq. (3.51 1 to obtain (t{lo). Now we Make the Fourier transform of Eq. (3.63 1, 



dxe 



ikx 



dre*'^"^(eni^)2(n(x,T)n(0,0)). 



(3.64) 



Thus the ac conductivity (at A; = 0) of a one-dimensional system without umklapp 
scattering gives 

2e2 



(t{uj) 



2 

2' 



-i{ijO + i6){(f>*{k = 0,U!n)(t>{k = 0, Wn)) 



. , tt2uK 
;i[uj + id) 



IT- 



iuJn — >uj~\-i6' 



(3.65) 



(3.66) 



Here we have used again the functional integration formula ( 3.54| ) to calculate the 
{(j)*4') correlation using the (p part of the action given in (3.60 1 (and Z^) at A; = to get 
the long wavelength ac conductivity. Finally, after analytical continuation we obtain 



, , 2iuK 9, 
a{uj) = = e^{2uK) 



+ iP- 



TTCJ 



(3.67) 
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where P denotes the Cauchy principal value. As expected, the one-dimensional system 
is a perfect conductor with an infinite static conductivity given by a delta function peak 
at a; = (Drude peak). The weight of this peak is 2e'^uK. 

If we now include umklapp scattering in the system we would no longer have a 
perfect conductor. Such scattering processes are responsible for the appearance of a 
finite resistivity because they produce momentum relaxation. The next section is ded- 
icated to the obtention of the conductivity in a ID system in the presence of umklapp 
scattering. 



3.3.2 ID systems with umklapp scattering 



As we saw in Sec. 2.4.2[ the umklapp processes arise due to the presence of the lattice 



in systems with commensurate fillings. Here we consider a 1/2-filled system, but the 
subsequent derivations can be done for any commensurate filling [,16.| . In the following 
we make use of the memory matrix formalism presented in Sec. |3.2[ because the cal- 
culation of the exact correlation function is not feasible when umklapp scattering 
is present in the system. Let us first write down the expression for the conductivity in 
the memory function formalism, Eq. ( |3.31 1, in bosonization language 



a(uj) = — . , ^, , ■ (3.68) 

^ ^ vr uj + iM{oj) 

In a one-dimensional system M{io) is just a scalar function, as the one discussed in 
Sec. 3.2 given by[^ 



-x(o) 

where the force (the subscript u is there to avoid confusion with the LL parameter 
Kp) is Ku = \H.u^ J\ and Hu is the umklapp scattering term in a 1/2-filled system, 

Uu = 77^^ / a!xcos(^/80p(x)). (3.70) 
(2vra)^ J 

From the above expressions it is clear that if the current commutes with the Hamilto- 
nian, as in a ID system without umklapp scattering, the memory function is zero and 
one recovers a perfect conductor. 

The operators are easily computed calculating the commutator between the 
umklapp operator (3.70 1 and the current ( 3.48| ), using the commutation relations (2.161 



between bosonic fields 4> and Q (remember that 11 = V^/vr), 

= \Hu,A = (^^fj2(^P-^p)^ y c^a;sin(\/8(/)p(x,r)). (3.71) 
The next step is to calculate the correlation {Ku^Ky) and then insert the result in 



Eq. (3.69 1 to finally obtain the conductivity. We will not compute here the full expres- 



sion for the correlation {Ku, Ku) which is obtained in detail in Refs. ll23llT6ll . Instead, 



'This expression for M {u) differs from tire expression in reference 1 16| by a factor i. We make tliis 
cfioice of notation to be in accordance with the 2D case where the memory matrix is defined as iMai3{uj)- 
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Figure 3.1: ac conductivity in a one-dimensional Mott insulator at commensurate fill- 
ings. At frequencies smaller than the optical Mott gap 2A, the conductivity is zero and 
at frequencies larger than 2A, it decays as a power-law with an interaction-dependent 
exponent fj, — 3 — An^Kp, with n the commensurability order (From Ref. [17j ). 



we use scaling arguments to determine the exponents of the power-law frequency (or 
temperature) dependence of the conductivity. Each force depends linearly on the 
parame ter g3 , thus the memory function is of order g^. The function sm{\/84>p{x, r)) 



in Eq. (3.71 1 behaves at large distances as 



16 ] (where the cutoff a is of the 



order of the lattice parameter), giving for the correlator 



/3 



dx / dre*""^(T,i^„(x,r)E:„(0,0)) 



■*Lj+i5 



(3.72) 

The factor —2 in the exponent comes from the integrals in r and x. Finally, we obtain 
the following frequency dependence for the memory function 

1 



iM{uj) 



g\—uj 



(3.73) 



At high frequency (u » T) the expression for the conductivity ( 3.68| ) must be ex- 



panded in Lu^^: a-{Lu) ~ x(0) [l/i^ + iM{u)/u>'^ + •••]. Therefore, (t{uj) is given 
by a power-law dependence with an non-universal exponent depending on interactions 
(typical behavior of a LL): 

aioj) ~ ^uj^^p-^. (3.74) 



Result ( 3.74 1 is valid as long as a perturbative expansion in 513 is reasonable, i.e. , when 
the umklapp operator is irrelevant (or marginal) |[T6l . In this case the system is a per- 
fect conductor with a regular part given by the power-law dependence of Eq. ( |3.74 ): 
Re[cr(i^)] = DS{lo) + (Treg(w). There is still a peak at w = (Drude peak). In the op- 
posite case, if the umklapp operator is relevant and leads to the opening of a gap in the 
excitation spectrum, the above procedure will only be valid at energy scales larger than 



the gap. For instance, in a Mott insulator Eq. (3.74 1 is only valid for frequencies larger 
than the optical Mott gap {lo » 2A). The conductivity is thus strongly affected by the 
Mott transition. As mentioned before, we consider here only the case of the Mott-U 
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T 



Figure 3.2: Temperature dependence of the dc resistivity in a one-dimensional Mott 
insulator with commensurability n. At temperatures below the Mott gap A, the number of 
carriers is exponentially small giving an exponential increase in resistivity. At T > A the 
dc resistivity shows a power-law behavior typical of a Luttinger Liquid (From Ref. lil7J ). 



transition where the ID Mott insulator becomes a Luttinger Liquid (metallic phase) at 



a critical value of the interactions Kc = Xjr? (see Fig. 2.9 1 [16|. Fig. 3.1 shows the 
zero temperature ac conductivity in a Mott insulator. The ac conductivity is zero for 
frequencies smaller than the optical Mott gap 2A. At lij > 2A, the conductivity decays 
with the power-law given in Eq. p.74| ). 

The temperature dependence (at a; ^ T) or dc conductivity is obtained by the 



same scaling method, giving cr(T) ~ T'^~'^^'' j g\. Fig. 3.2 shows the temperature 
dependence of the dc resistivity, which is related to the conductivity as p = l/cj, in 
a ID Mott insulator. At temperatures smaller than the Mott gap A, the resistivity 
increases exponentially, and for temperatures larger than A there is again a power-law 
behavior 

p[T) ^ (.3 

At the Mott transition, Kp = Kc = the ac conductivity and dc resistivity lead 

back to universal exponents: ct^lu) ln(u;)^) and p{T) ~ r/ln(l/T)^, respec- 

tively fT6l. 

Now that we have studied the transport properties of one-dimensional systems, 
especially those of ID Mott insulators, we acquired the necessary tools to investigate 
the conductivity of quasi ID systems. These transport properties are essential for the 
study of the Hall effect in weakly coupled Luttinger liquids in Chapter|5] 
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3.3.3 Transport in quasi one-dimensional systems 



In quasi one-dimensional systems, transport properties change from the pure ID case. 
In addition to the conductivity along the chains, which remains essentially the same as 
studied before, there is also a transverse conductivity a± . This transverse conductivity 
is obtained at high temperatures or frequencies by making a perturbative expansion in 
the interchain coupling term, or by using a mean field approach Il25ll26]| . In this section 
we will use again a scaling analysis to obtain the frequency (temperature) dependence 
of the transverse conductivity a± and we refer the reader to the literature L25t i26il for 
the complete expressions. 



In order to obtain a± we use again the Kubo formula (3.15 1. Let us consider a 



system composed of weakly coupled ID chains described by Hamiltonian of the form: 
n = Tiw + W_L (see Eqs. pll) and ^232) ), 



n 



dx 



C. 



(3.76) 



where j is the chain index. We suppose a system with spin rotation symmetry gi_i_ = 
= (Ka = 1). We will keep the fermionic representation (to maintain the same 
representation of the coupling term Tl±), knowing that the terms corresponding to 



the ID chains (first and second terms in Eq. (3.76l) are easily bosonized in the form 



( 2. 19 1, as explained in Sec. 2.4. 1 With this Hamiltonian we can easily calculate the 
diamagnetic term in the transverse direction using definition ( 3.16| ), 



-2e'^t±al 



dx 



(^J^(x)V'iT(a;)) +h.( 



(3.77) 



with Qy the lattice parameter in the transverse direction. The thermodynamical average 



in Eq. (3.77 1 must be obtained to first order in the coupling term (first order mt± « 
1), using standard perturbation theory [2]. For this we must expand the action in the 
thermodynamical average to first order in t± 



{Wo^{x)4^i^{x)) 



-t± 



/3 



',t 



(3.78) 

The thermodynamical averages on the right hand side of Eq. ( |3.78 1 are taken with re- 
spect to the ID part of the Hamiltonian. The quantity (V'j(x)'i/'t(a;)) is just the equal- 
time Green's function of the ID chains, who's Fourier transform gives the occupation 
factor n{k) ~ max[(5fc, r]''+\ with r]= l{Kp + K'^) - \, (see Eq. |2.28| ). Using sim- 
ple scaling arguments (we must add a factor to account for the x and ri integrals), 
we find the temperature dependence of the transverse diamagnetic term 



X±(0) 



^2 rp-2rp2ri+2 



^2 j,-l+l{Kp+K-^) _ 



(3.79) 



Note that the scaling analysis does not give the correct dimensions and non-interacting 
Umit for x±(0)- The detailed calculation gives ~ t\J^ [25.1 , but we are interested 
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here in the power-law dependence and for this the scaling analysis is sufficient. Re- 
sult (3.79) will be used in Chapter |5] to study the temperature dependence of the Hall 
coefficient. 

To obtain the conductivity a± we need the current-current correlation function 
X±{^) defined in Sec 3.1 The current along the transverse direction is obtained using 
definition ( 3.19 1 and Hamiltonian ( 3.76 ), 

J± = -iet±ay / ['>p]A^)^j+iai^) -h.c.^ . (3.80) 

Thus, the current-current correlation function in frequency is given by (see definition 



Km) 



(3.81) 



dre'^--{Trj[{r)J^m 

- iijji — >uj+i5 

Because each current operator contributes with a t_\_, this term is already of second 
order in this parameter {t<^ << 1) and thus the averages can be taken considering only 
the ID part of the Hamiltonian. Taking r > 0, we have 



(T.Jl(r)J^(O)) 



-ehlal j dx((^]+iJx,T)V,+iJO,0))(Vj,(0,0)^^.Jx,r)) 
-ehWy [ dx Graix,T) [1 - G,,(x,r)] , (3.82) 

D T 



r=R,L 

with G{x, T)ra ~ max[aj, T]''+^, with r/ = \{Kp + K~^) — ^, the Green's function 
of a one-dimensional system, with spin rotational symmetry (Ka = 1). Now we can 
apply scaling arguments to get the frequency (temperature) dependence of the current- 
current correlation function 



(3.83) 



Using the Kubo formula for the transverse conductivity a± = [x±(0) — x±('^)]/^'^' 
the result is, as for the longitudinal conductivity, a power law 



(3.84) 



E* is the energy scale of the dimensional crossover: at T, u; < E* the system is in 
a high-dimensional metallic state, whereas at T, u; > E* the system is in a Luttinger 
liquid state. A full expression of a± can be found in Ref. ESll . At this stage we have 
a general knowledge of the transport properties in low-dimensional system and, in the 
next chapter, we will take the challenge of studying the Hall effect in these systems. 
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CHAPTER 4 



The Hall effect 



4.1 The classical Hall effect 

The Hall effect was found by E. H. Hall in 1879 using an experimental setup similar to 
the one shown in Fig. |4.1| An electric field is applied along a wire extending along 
the X-axis and a magnetic field B is applied along the z-axis. The electric field drives 
a current density jx along the x-direction and the magnetic field deflects the electrons 
in the negative y-direction with a Lorentz force equal to 

F = -i^t-xS, (4.1) 

c 

where |e| is the electronic charge and v is the average velocity of the charges. Because 
current cannot flow along the y-axis, the charges accumulate on the sides of the wire. 
As they accumulate, an electric field builds up in the y-direction that opposes the mo- 
tion of the charges and their further accumulation. At equilibrium, this transverse field 
Ey balances the Lorentz force. The Hall resistivity pyx is defined as the ratio between 
the transverse electric field and the current along the wire ||2l 

Pyx = ^. (4.2) 



J, 



X 



The Hall coefficient Ru (or Hall constant) is just the Hall resistivity divided by the 
magnetic field. 

Assuming that the motion of the electrons is classical, their equation of motion reads 

vxB 



mv 



E + 



— (4.4) 
r 



where r is the relaxation time for scattering, that was defined in the description of the 
Drude model, at the beginning of the previous chapter. In the steady state the time 
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Figure 4.1: Hall's experiment: an electric field is applied in a metallic sample along 
the X-axis and a magnetic field B is applied along the z-axis. The electric field drives a 
current density along the x-direction and the magnetic field deflects the electrons in 
the negative y-direction with a Lorentz force equal to —{\e\/c)vxB (v is the charge's 
velocity). The electrons accumulate on the sides creating an electric field Ey in the y- 
direction. At equilibrium, this transverse electric field balances the Lorentz force. 



derivatives ai^e zero and we have 

= 
= 



VyB' 


mvx 


c 


T 


VxB 


mvy 


c 


T 



(4.5) 
(4.6) 

The current density along the x-axis is given by jx = —\e\nvx. Furthermore, the 
current density along the y-direction is Vy = 0. Then, multiplying Eqs. (4.5 1 and (4.6 1 
by n|e|r/m we obtain 

ixB 

jx, Ey = -^ (4.7) 



n\e\c 



where a = n\e\'^T /m is the Drude model dc conductivity in the absence of magnetic 



field (see Eq. (3.2i). Then, the Hall coefficient results 

-Rh = r-;-- 



n e c 



(4.8) 



As can be seen in Eq. (4.8 1 the Hall constant does not depend on the parameters of the 
metal or the dimensions of the sample, except for the density of carrier n. As a result, 
i?H gives a direct measurement of the number and sign of the carriers in the metal. 

There are additional phenomena, related to the Hall effect, known as the quantum 
Hall effect and the fractional quantum Hall effect that occur in highly correlated elec- 
tron systems where electrons are allowed to move only in a plane and a strong magnetic 
field is applied In this work we study the regime where the applied magnetic field 
B is supposed weak enough for the linear response theory to be valid and the quantum 
effects do not appear. 

The Hall constant result obtained in the classical approximation assumes a system 
where electrons move as free particles. Therefore, it can only be used in metals where 
the free-electron model works well. In order to understand the Hall effect in correlated 
systems, we must take into account interactions, which is generally a complicated task. 
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The next sections are devoted to the investigation of two different approaches for the 
problem of the Hall effect in the presence of interactions. The first approach consist on 
a high-frequency expansion where the Hall coefficient is obtain in the limit of infinite- 
frequency. And the second approach is the obtention of Ru using the memory function 
formalism studied in Chapter [3] 



4.2 The infinite-frequency Hall constant 

The conductivity was expressed, in the previous chapter, by the tensor cr^jy 



'^yx ^yy 



(4.9) 



The Hall resistivity pyx = cf^, written in terms of the conductivity tensor, is equal to 



Pyx 



a 



xy 



^xx^yy ~1~ ^xy^yx 



(4.10) 



To obtain a result for pyx by means of Eq. (4. 10 1 when interactions are present in 
the system is generally not feasible, thus some approximations must be made. In 
the following, we develop a high-frequency expansion of the conductivity, in order to 
obtain a simpler expression for the ac Hall constant, as was proposed by Shastry et al. 
in Ref. |27l. Let us write the conductivities using the Kubo formula Eq. (3.15 1. The 
longitudinal conductivity along the x-direction is 



(^xx^^) 



1 



lUJ 



[Xx(0) - Xxx{uj)] , 



(4.11) 



where XxxC"^) is given by Eq. (3.17 1 at g = 0. Without lost of generality we can take 
t' = in Eq. (3.17 1. Integrating by parts, we have 



Xxx{uj) 



ijj 



([Jx(t),Jx(0)]) 



+ i 



dt- 



([[H,Jx(t)], Jx(0)]). (4.12) 



We have used the Heisenberg representation for the current operator J(t) = e J(0)e 
with 7i the total Hamiltonian of the system. The first term in ( 4.12[ ) vanishes because 
the currents commute at t = and, as we said at the beginning of Sec. 3.15[ the per- 
turbation (here the electric field) is assumed to vanish when t —>■ itcxD, implying the 
vanishing of the current. Integrating by parts again, the second term on the right side 
of Eq. ( |4.12 1 becomes 



Aujt 



Xxx{uj) 



([[Jx(t),W],Jx(0)]) 



+ i 



dt — 



([[W,[Jx(t),W]],J..(0)]). 

(4.13) 

As before, the current vanishes at t = oo and therefore only the i = contribution 
remains in the first term. Now, we can write the high-frequency expansion of (Txx(i^) 
as 

{[[Jx,n]^Jx\) 



1 



Xx(0) + 



(4.14) 
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The longitudinal conductivity in the y-direction, ayy{uj), has the same high-frequency 
expansion with Xa;(0) and replaced by Xj/(0) and Jy, respectively. From the Kubo 
formula we have for the transverse conductivity 



(4.15) 



Integrating by parts the definition of Xxy{^), as before, the following high-frequency 
expansion results 



C^xyiyj) 



I 

UJ 



UJ 



(4.16) 



We can now insert the frequency expansion of the conductivities in Eq. (4.101. Since 



the conductivity Gxy is proportional to the magnetic field, the factor Oxy<yyx in the 



denominator of Eq. (4.10) is of second order in B and thus can be neglected. Keeping 



only the linear term in B and zero-order term in uj, we obtain for the Hall constant 



i {[Jx,Jy\) 

'SXx(0)Xj/(0)' 



(4.17) 



This expression for is known as the infinite-frequency Hall constant ||27|. We 
will use it in Chapter [5] to calculate the high-frequency Hall coefficient on a triangular 
lattice. For this high-frequency expansion to be valid in an experimental measurement, 
it requires a probe frequency larger than any other energy scale in the system, and it 
can be measured using the Faraday rotation experiment. 



Although Eq. (4. 17 1 seems to be a much simpler expression for the Hall coefficient 



when comparing with (4.10), we must remember that when interactions are present in 
the system, they must be taking into account when calculating the thermodynamical 
average (•••). This is usually not a simple task. If interactions are small, they can be 
treated using perturbation theory; if not, other formulas to compute iin must be used. 
In the next section, we will use the memory function formalism, studied in Sec.|3.2[ to 



obtain an expression for the Hall constant at leading order in the interaction term. 



4.3 The Hall constant in the memory function approach 



As in the previous section, our main goal is to calculate the Hall resistivity py^, and 
consequently Ry^, in terms of the conductivity tensor a^y using Eq. (4. 10 1. However, 



in this section we ought to obtain an expression for the Hall constant where interac- 
tions can be included directly in the calculations, going beyond the infinite-frequency 
limit. For this we use of the memory matrix formalism. With the definitions given 



in Eqs. ( 3.40 )-( 3.43 ) we arrived at an expression for the off-diagonal memory matrix 



element in terms of the conductivities 



iMxy{uj) 



iXy{^)<^xy{^) 



^xx 



^xy = iXy(^)Pxy{yj) - ^xy 



(4.18) 



From this expression it is straightforward to rewrite the Hall coefficient, Ru = p^y/B, 
as 

Qxy + iMxyiuj) 



Ruiuj) 



1 



lim 



B 



(4.19) 
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Since the memory matrix vanishes as uj~'^ at high frequency, we see from Eq. (3.42i 
that the infinite-frequency Hall constant, Rui^^ oo), is indeed given by Eq. (4.17 1. 
The result ( 4. 19| ) will be used to compute the Hall constant in a strongly correlated 
quasi one-dimensional system in Chapter [5] 

Let us now return to the memory matrix element necessary to calculate Ru, given 
in Eq. ( [3l3] ), 

[Xx{0) - Xxx{uj)] [x?/(0) - Xyyi^)] 

This expression will be expanded at high frequencies to obtain the result of iMxy{uj) 
given in Eq. (3.461. At high enough frequencies, where the quantity |x/x/x('^)/XAt(0)l 
can be considered small, the off-diagonal memory matrix term can be written as 



iMxy{uj) 



xy 



(4.20) 



iMxy{uj) ~ 



^Xxy{^) 

Xx(0) 



1 + 



X.(0) Xs/(0) 



17 



xy 



(4.21) 



Owing to the equations of motion given in Eq. (3.22i, for the susceptibility Xxy{'-^), 
the above expression becomes 



{[Jxj Jy]) 

Xx{0) 



{[7i, Jx];Jy) 

XxiO) 



1 + 



Xx 



Xx{0) Xy{0) 



{[Jxj Jy]) 



Xx{0) 

(4.22) 

The two terms involving [J^ , Jy] cancel. In the following we will write the Hamilto- 
nian asH = Ho + Hint, 



Xa 



iMxy{u;)^nxy ^^^^^ . ^^(0) 

{[Ho, Jx]; Jy) , {[Hint, Jx]', Ji 



(4.23) 



Xx{0) 



+ 



Xx{0) 



1 + 



Xi 



X^) _^ Xyyi^) 



Xx{0) Xy{0) 



Because there is always a magnetic field present in the Hall experiment, the non- 
interacting part of the Hamiltonian Hq does not commute with the currents. How- 
ever, following Mori's formalism II22I . the commutator [Ho, J^] can be expressed as 
1281 ED 



[Ho, JfM 



-j^n' 



(4.24) 



with fi,v = x,y and summation over repeated indices is implied. The superscript 
means that the average (• • • ) must be computed using the non-interacting Hamiltonian 
Hq. With this expression as well as the symmetry property Xx{0)^xy = —Xy{^)^yx, 
Eq. (4.23 1 can be rewritten in the form 

Xxxiu; 



XxiO) 

{Jy, Jy) (['^int; Jx]^ Jy 



Xyyi^) 
Xy{0) 



(4.25) 



1 + 



Xx 



"(^) _^ Xyy{i^) 



Xx{0) Xy{0) 



Xy{0) Xx{0) 

The second term in the right hand side is cancelled by the first term in the square 
brackets. Now let us keep only the leading terms in |x^^(^)/Xm(0)I (the correlator 
([Hint, Jx]; Jy) is of the order of x^j^t), 

{[Hint,Jx];Jy) , Xxx{^^) 



iMxy{uj) 



Xx(0) 



+ 9. 



xy 



Xx(o) ■ 



(4.26) 
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As before, we use the equation of motion for the correlator ([Hint, Jx] ; Jy) and for the 
susceptibility Xxx{^)- Introducing the residual force operator = [Hm\.,Jx], we 
have 

~ '''''' u;Xx{0)Xx{0)~''''^J^ 

The second term and last term in the right hand side cancel. Each frequency matrix 
is of first order in the magnetic field because the commutator [J^, Jy] is proportional 
to B. Thus, the fourth term above is of second order in the magnetic field. Then, at 
high frequency, we end up with two terms for the off-diagonal memory matrix element 
at first order in B 

iMxyiuj) ~ - (4.27) 



The second term in Eq. (4.27 1 is of first order in the interactions parameters. How- 



ever, it can be proved that the leading terms of the memory matrix are of second 



order in the interaction parameters, as obtained in Sec. 3.3 for a ID system. In or- 
der to see this, let us take a system where electrons interact via a Hubbard interaction 
[7^- ct^Cj|. The correlator {[K^, Jy]) correspond to a first-order expression of the 
frequency matrix (22 ]. The frequency matrix is easily traced back to the number oper- 
ator njcr = (c|o-Cj^), whose first-order contribution in U vanishes. In consequence, the 



only remaining term in Eq. ( |4.27 1 is of second-order in the interaction parameters, 

^Mxy{u) ^ (4.28) 
This is the expression that will be used in the following chapter to compute the mem- 



ory matrix contribution to the Hall constant Ru, by means of Eq. (4. 19 1. Because 
Eq. ( |4.28| ) is already given at second order in the scattering parameters, the thermo- 
dynamical average of K operators can be computed with the free Hamiltonian, i.e., 
putting interaction terms to zero. This facilitates enormously the calculations, because 
these averages are usually straightforward to compute without interaction. Even more 
in one-dimensional systems, as explained in the preceding chapters. 

At this point we have reviewed all the theoretical tools necessary for our study of 
the Hall effect in two different models of strongly correlated systems. In the next chap- 
ters, we will see how all the previously studied formulas can be applied to an specific 
theoretical model or to Hall measurements made in real experimental compounds. 
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CHAPTER 5 



Hall effect in strongly correlated quasi 1 D systems 



In the last decades, various experimental realizations of low-dimensional systems have 
been achieved. Realizations as the organic conductors [291, carbon nanotubes [30], ul- 
tra cold atomic gases |3r|, quantum wires [321, quantum dots [[33, 34] and others, have 
largely stimulated the research in low-dimensional physics. Among these, the organic 
conductors have become the model systems for the study of quasi one-dimensional 
physics due to their highly anisotropic molecular structures. They have been exten- 
sively studied for more than twenty years now, since the discovery of a superconduct- 
ing state in their phase diagram ||29]| . We will start this section with a short review on 
the physical properties of these organic conductors and some experimental facts that 
motivated the theoretical study of the Hall effect in these strongly correlated quasi ID 
systems. 



5.1 Quasi one-dimensional organic conductors 

Within the various families of organic conductors [35J we will focus here in the prop- 
erties of Bechgaard (TMTSF-X) and Fabre (TMTTF-X) salts, which are very suitable 
compounds for the study of quasi- ID physics, as will see in the following. TMTTF-X 
stands for tetramethyltetrathiafulvalene and TMTSF-X for tetramethyltetraselenaful- 
valene. In these formulas, the X denotes an inorganic anion with various possible 
symmetries: spherical (PFg, AsFg, SbFg, TaFg), tetrahedral (BF4, CIO4, Re04) or 
triangular (NO3). TMTTF-X and TMTSF-X salts belong to a same family forming 



the generic (TM)2X phase diagram shown in Fig. 5.2 These salts have three main 
conducting directions: the stacking direction of the molecules called a, the transverse 
direction called b and the direction c, perpendicular to the ab plane. The hopping in- 
tegrals are of the order of ta = 3000K, tb = 300K and tc = 20K. The molecular 



structure of (TM)2X conductors is depicted in Fig 5.1 With such anisotropy in their 
molecular structure, where the overlap between electron clouds along a is 10 times 
larger than the overlap between the stacks in the transverse b direction and 150 times 
larger than the one along c, the electronic structure can be seen as one-dimensional 
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Se Se 




Figure 5.1: a) TMTSF molecule; b) view along the stacks, a direction; c) Along the c 
direction the stacks are separated by the X anions (PF^T in the picture). 



with a slightly warped Fermi surface (see Fig. 2.11 1. This is what make these com- 
pounds quasi ID materials. 

There exist a variety of ordered states in the phase diagram of (TM)2X compounds 
(see Fig. 5.2 1: spin-Peierls (SP), spin density wave (SDW), charge localized (loc), 
charge ordered (CO), antiferromagnet (AFM) and superconducting state (SC). In ad- 
dition, there is a metallic state which description changes from a one-dimensional 
Luttinger liquid to a two- or three-dimensional Fermi liquid (dimensional crossover 
studied in Sec. 2.4.3 1. Going from left to right in the phase diagram, the materials get 
less one-dimensional due to the increasing interaction in the second and third direc- 
tions. In order to study the Hall effect in these systems, we will concentrate on the 
properties of the "normal phase"(metallic state) at high temperature, where the one- 
dimensional behavior appears. We refer the reader to the literature [29] for a review 



on the different ordered states appearing in the phase diagram of Fig. 5.2 



The insulating properties of the TMTTF compounds evolve towards those of TMTSF, 
which are good conductors, trough an insulator-metal transition when increasing pres- 
sure (or changing the anions X), as can be seen in Fig. |5.2| The insulating behavior of 
the TMTTF family, as in (TMTTF)2PF6 compound, is expected for a one-dimensional 
Mott insulator. Such Mott insulator behavior could come from the 1/4-filled nature 
of the band or from the 1 /2-filled nature, due to a small dimerization existing in these 
molecules, of the order of ~ lOOK |29|. This indicates that interactions have a 
large effect in the properties of TMTTF family. In the case of the TMTSF family, due 
to their metallic behavior at ambient pressure, the role of interactions is more compli- 
cated to understand. 

Studies of the longitudinal transport (a direction) have revealed signatures of LL 
properties [3711381 129] . Fig. |5.3| shows the optical conductivity along the chain axis (a- 
axis) in the TMTSF family, where a power-law decay is found for cr(w), in agreement 
with the predicted conductivity in a Luttinger liquid (see Eq. ( 3.74| )). The exponent 
of the power-law behavior allowed an experimental determination of the Luttinger pa- 
rameter Kp ~ 0.23 1,37,1 . consistent with an interpretation of the insulating state as a 
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Figure 5.2: Phase diagram for the (TM)2X compounds. For the different compounds 
the ambient-pressure position is indicated. Here loc stands for charge locaHzation, CO 
for charge ordering, SP for spin-Peierls, AFM for antiferromagnet, SDW for spin density 
wave, and SC for superconductor The clear phase transitions are indicated with solid 
lines and the dashed lines indicate crossovers. From Ref. If36l . 



quarter-filled Mott insulator in the TMTSF family. Photoemission data has revealed 
a similar value for Kp, but in a very large range of energies ll39i . making these re- 
sults more difficult to interpret. Transport transverse to the chains has given access 
to the dimensional crossover between a pure ID behavior and a more conventional 
high-dimensional one f40l l4n[T7ll29l . Optical conductivity measurements give a di- 
rect evidence of a deconfinement transition between a one-dimensional insulator and 
a high-dimensional metallic regime, when the observed gap is of the order of the in- 
terchain hopping ll39l . Finally, measures of resistivity along the chains and spin sus- 
ceptibility, have shown evidence of the spin-charge separation characteristic of a LL 

m. 

To probe further the consequences of correlations in these compounds, several 
groups have undertaken the challenging measurement of the Hall coefficient Ru{T) 
||42l|43l|44l|45l. In particular, two measurements of Ru{T) were made in 2000 by 
Moser et al. Bll and Mihaly et al. [43], in the organic conductor (TMTSF)2PF6. 



The results are shown in Figs. 5.4 and |5.5[ respectively. In the former, the current 
is applied along the a axis and the magnetic field along the c axis. Thus, the Hall 
voltage develops along the b axis. They investigated the temperature dependence of i?H 
between and 300 K. In the normal state of (TMTSF)2PF6, at temperatures T > 130 
K (a dimensional crossover to a 2D metal is expected at T ~ 130 K), the Hall constant 
was found to be temperature-dependent (see Fig. 5.4 1, increasing with T. The sign 
of i?H was found to be positive (holelike). In the second Hall experiment, shown in 
Fig. 5.5 the current was applied along the c direction and the magnetic field parallel to 
the most conducting direction (a axis). The Hall voltage was measured along b. They 
obtained a Hall constant independent of temperature for T > lOOK and a positive sign 
for Ryi. These results, different depending on the direction of the applied magnetic 
field, proved difficult to interpret due to a lack of theoretical understanding of this 
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Temperature (K) 

Figure 5.4: Temperature dependence of the Hall constant in (TMTSF)2PF6 measured 
by Moser et al. in Ref. [42 1. The current is applied along a and the magnetic field along 
c. The Hall voltage is measured along h. The dashed fine is a power law fit with 
a = 0.73. 
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Figure 5.5: Temperature dependence of i?H in the normal phase of (TMTSF)2PF6 mea- 
sured by Mihaly et al. in Ref. B3l . The current is applied along c and the magnetic field 
along a. The Hall voltage is measured along b. 



problem. This prompted for a detailed theoretical analysis of the Hall effect in quasi- 
ID systems. A first move in this direction was reported in Ref. ll25l where the Hall 
coefficient of dissipationless weakly-coupled ID interacting chains was computed and 
found to be T-independent and equal to the band value iZ^ = 1/ nec. This surprising 
result shows that in this case Ru, unlike other transport properties, is insensitive to 
interactions. However the assumption of dissipationless chains is clearly too crude 
to be compared with realistic systems for which a finite resistivity is induced by the 
umklapp interactions ll23l . 

This chapter is the object of publications f46l and (4T\. In the following we ex- 
amine the effect of umklapp scattering (see Sec. 2.4. 1| ) on the temperature dependence 
of the Hall coefficient in quasi-lD conductors, and we discuss the applications to the 
Hall experiments mentioned above. 



5.2 The Hall effect in weakly coupled Luttinger liquids 

We consider a model composed of weakly coupled 1/2 -filled ID chains. We take 
a 1/2-filled band because the umklapp scattering at 1/4-filling is much more compli- 
cated to treat in the calculations. But understanding the effect of the 1/2-filled umklapp 
scattering already gives an idea of the 1 / 4-filled case. With this model we compute 
Rh{T) to leading order in the umklapp scattering using the memory function approach 



explained in Sec. 3.2 112111471 . We find that umklapp processes induce a T-dependent 



correction to the free-fermion value R^. This correction decreases with increasing 



temperature as a power-law with an exponent depending on interactions (Fig. 5.7 1. At 
the end, we discuss the implications for quasi- ID compounds. 
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Figure 5.6: Schematics of the model. The chains and the current / go along the x-axis, 
the magnetic field B is applied along the z-axis, and the Hall voltage is measured along 
the y-axis. 



5.2.1 Model and methods 



Our model is sketched in Fig. 5.6 We consider ID chains coupled by a hopping 
amplitude t± supposedly small compared to the in-chain kinetic energy. As explained 
in Sec. |2.4.l"| the usual Luttinger liquid model of the ID chains assumes that electrons 
have a linear dispersion with a velocity vp. For a strictly linear band, however, the Hall 
coefficient vanishes identically owing to particle-hole symmetry. A band curvature 
close to the Fermi momenta ^kp is thus necessary to get a finite Ru. We therefore 



take for the ID chains of Fig. 5.6 the dispersion 



(5.1) 



The upper (lower) sign corresponds to right (left) moving electrons. Eq. (5.1 1 can be 
regarded as the minimal model which gives rise to a Hall effect, while retaining most 
of the formal simplicity of the original LL theory, and its wide domain of validity. 
In particular, this model is clearly sufficient at low temperatures (compared to the 
electron bandwidth) since then only electrons close to the Fermi points contribute to 
the conductivities. 

Our purpose is to treat the umklapp term perturbatively. We express the Hamilto- 
nian as Ho + Ti-u where Hu is the umklapp scattering term and Hq reads 



-ieA 



+ h.c. 



. (5.2) 



P.' A ■ dl. We 

J 1 



Bxay with ay the interchain 



In Eq. (5.2 1 j is the chain index, T3 is a Pauli matrix, and Ajj 
choose the Landau gauge Ay = Bx, such that Aj jj^i 

spacing, if)^ = (■i/'Jj V'l) ^ two-component vector composed of right- and left- 
moving electrons. The second term in Eq. (|5.2[) is the band curvature, the third term 



is the forward scattering and the last term corresponds to the coupling between the 
chains (Eq. ( |2.32 ) with a Peierls phase due to the presence of the magnetic field). As 
mentioned before, we omit the backscattering terms {gi processes) which are, for spin 
rotationally invariant systems, marginally irrelevant f\6\. We therefore take gi_\_ = 
= 0. At 1/2 filling the umklapp term reads (the bosonized version is given in 
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Eq. (2.311) 



nu = f [dxY^ {i^ln4-aRi^JaL^j,-.L + h-c) • (5.3) 

We will compute the ac Hall constant using the memory matrix formalism, where the 



Hall coefficient Ru = Pyx/B is given by (see Eq. (4.19 1) 

1 



iXy{0) 



lim ^""^ ~^ iMxyiiv) 



B 



(5.4) 



From Hamiltonian (5.2 1 and the definition for the diamagnetic term, given in 



Sec. 3.1 we obtain the longitudinal and transverse diamagnetic terms 



Xx{0) 
Xy{0) 



2eH±al / dx(V'5t(x)^it(x)e-'^-^"«^+h.c.). 



(5.5a) 
(5.5b) 



For the longitudinal diamagnetic term (5.5ai we have used the relation between the 
electron density n and the Fermi momentum: nuy = kp/n (since nay is the density 
per one chain). It is also easily obtained from Eq. (3.491 with UpKp = vp. For the 



evaluation of the frequency matrix in Eq. ( |5.4| ), we write down the current operators, 
making the functional derivatives of the Hamiltonian with respect to the vector poten- 
tial, as explained in Sec. |3.1[ 



Jx = e j dxY2'^]cTi^)[^F^3 + '^a{-idx)ri]4jj^{ 
Jy = -iet±ay J dxY^ (vj^V^j+^^.e 



^•^■+1 - h.c. 



(5.6a) 
(5.6b) 



To obtain the frequency matrix, we must calculate the commutator between these cur- 
rent operators. Thus, using the standard commutation relation for fermionic operators. 



the expression resulting from Eq. ( 3.42 1 for i^xy is then 

"y"" I /„/,t l^\^-ieBayX 



xy 



2TTeat±ayB 

>j 



dx {'^pQ^{x)^l^^{x)e 



+ h.c.). 



(5.7) 



At this stage we can already evaluate the high-frequency limit of Ru, because the 
memory matrix vanishes as if lij — > oo. Thus, the effects of the umklapp dis- 
appear at high frequency, and in this limit one recovers from Eqs ( 5.4f 5.7| ) the result 
obtained for dissipationless chains in Ref. Il25l . namely that the Hall coefficient equals 
the band value Ry^: 



Rn{oo) 



R'i 



vraa,. 



(5.8) 
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5.2.2 Rn in the presence of umklapp and forward scattering 



As explained in Sec. 4.3 the memory matrix element Mxy{uj), necessary to obtain 
i?H(w), can be reduced at high frequencies and linear order in B to the calculation of 



the following average 



iMxy{uj) 



1 {K^;Ky 



(5.9) 



with K^j^ the residual forces operators, which in this case are K^^ = \Hu^ </^], and 
{Kx^ Ky) stands for the retarded correlation function of the operators K^. Using the 
umklapp term (5.3 1 and the currents ( |5.6[ )p]we find 



= 2euF53 Ma; ^ (^V'jaij^j -crijV'i -aLV'jaL - h.c.^ (5.10a) 

■J I, r D 

(5.10b) 



jcr h=L,R 



Note that each of the K\ is of first order in 513, hence M^y is of order g\. The quan- 
tity {Kx\Ky) entering Eq. (5.9 1 is the real-frequency, long- wavelength limit of the 



correlator, which we evaluate as 



{K^\Ky) 



dT e'''^TrK,{T)Ky{Q)) 



(5.11) 



where denotes the Matsubara frequency. The correlator {K^] Ky), at first order in 



t^, vanishes for i3 = or a = as shown in Appendix A.l Thus, retaining only 



leading-order terms in t_\_ and a, the first nonvanishing contribution in Eq. ( |5.11| ) is of 
order at\glB, and involves three spatial and three time integrations, which we were 
not able to perform in full. Based on a scaling analysis, we can nevertheless extract 
the temperature (or frequency) dependence of this contribution. 

We evaluate the correlator {K^; Ky) to first order in a and t^. Let's denote by Tia 
the curvature [second term in Eq. (|5.2|)], by T-L±_ the inter-chain hopping [fourth term in 



Eq. (5.2 1], and by TIq the remaining parts of the Hamiltonian, Tiio = 'Hq — Tia — 'H±. 



Standard perturbation theory yields t2J 



{K,;Ky 



dr e 



ifir 



dndT2 {TrK,{T)Ky{0)n±{n)nair2)) (5.12) 



where the average is taken with respect to Hq. The latter corresponds to a ID chain 



and can be easily bosonized, as shown in Sec. 2.4. 1 With the help of representation 



'Here we use the following properties for commutators: [AB, CD] — A[B, CD] + [A, CD]B, and 
[A, BC] = [A, B] + C - B[A, C] + , where [A, B]+ denotes the anticommutator AB + BA. 
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( [2.15| l we bosonize each operator in Eq. ( |5.12[ ): 

Mevpgs 

I ax 

JO" 



K 



J dx sin (^V^(f)p{x)^ 



(5.13) 



^ (27ra)2 ^ ^ J " \ 



(5.14) 



where j and / are neighboring chains, b = 1) for right(left) moving fermions, 

and ejj' = (— 1)-' . For the coupling term we have 



t± 

27ra 



^-ieBayX^^{bM^)-9pix)+alb(f>„{x)-e„{x)]}. 



Y^jdx 

jab 

^^{h,j,,{x)-ep{x)+a[bcl,^{x)-e^{x)]}^^^ + h c ) 



and for the band curvature term we take B8l 



a 
2TTa 



dx 



(5.15) 



(5.16) 



Next we will use a very helpful identity for the calculation of correlators between func- 
tions of the fields cf) and 9. The following identity can be proved using the functional 
integral technique presented in Sec. [33] and its fully demonstrated in Ref. lfT6l 



exp { - ^ - [AnA^K + BnB.mK-^)Fi{ 



fn Tm ) 



+ {AnBm + BnAm)F2{r. 



n ' mj ( 1 



(5.17) 



where r = (x, ut), the notation J2' means that the sum is restricted to those terms for 
which J2n -^n = X^n = 0, and Fi 2 are universal functions which for (x, ut) ^> a 
are, 



Fi(r) 
F2{r) 



log 



-zArg 



sinh^(— — ) + sin^( 
pu 



7r2a2 V 



,'ITyas . . .TTX, 

tan(— — ) + ?tanh( — ) 
pu pu 



vrr , 

J' 



(5.18) 
(5.19) 



where ya = ut + aSign(r), and a is a momentum cutoff. The resulting expression for 
the correlator in Eq. ( 5.12| ) is (see appendix A.2 1 



{K,;Ky) ~ B j d\d\id\2e-^'^^^'^'\ 



(5.20) 



1^2 1 



The factor \r\ results from the linearization in the magnetic field B. In Eq. (5.201 we 
have discarded all factors involving the F2 function, since they correspond to angular 
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max(T*, A) 



Figure 5.7: Correction of the high-temperature/high-frequency Hall coefficient Rh by 
the umklapp scattering in weakly-coupled Luttinger liquids. R*^ is the value of the Hall 
coefficient in the absence of umklapp scattering, Eq. (5.8 1, and W is the electron band- 
width. Our approach breaks down below some crossover scale (dashed line, see text). In 



this figure we have assumed that A in Eq. (5.24 1 is negative 



integrals of the r variables and therefore do not contribute to the scaling dimension. At 
distances much larger than the cutoff a we have e~^^^^'^^ ~ {a/\r\)^, and therefore 
we find the high temperature, high frequency behavior as 



{Kx;Ky) ~ Bmax{uj,T) 



-?,+iKp-UKp~K-^) 



(5.21) 



As done in Sec. 3.3.3 we follow the same procedure for the diamagnetic term Xy(0) — 
however at zeroth order in a and B — and find 

Xy(0) ~max(w,T)-i+5(-^''+^^"'). (5.22) 

Combining these expressions and collecting the relevant prefactors we deduce 

1 {Kx;Ky 



(5.23) 



where Kp is the LL parameter in the charge sector (see Sec. 2.4. 1 1. In the absence 



of interactions we have Kp = 1, while Kp < 1 (Kp > 1) for repulsive (attractive) 
interactions. If the interactions are strong and repulsive (Kp <^ 1) the exponent in 
Eq. (5.23 1 changes due to the contraction |[T6l of the operators in K^ and Ky, which 



gives the relevant power-law in this case. Together with Eqs (5.9 1 and (5.4 1, Eq. (5.23 1 
leads to our final expression for the Hall coefficient: 



R-H 



3Xo 



(5.24) 



with W the electron bandwidth. 



Eq. ( 5.24 1 shows that in 1/2-filled quasi-lD systems the umklapp scattering changes 
the absolute value of the Hall coefficient with respect to the band value, which is only 
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recovered at high temperature or frequency. Note that Eq. ( |5.24 ) also describes the 
frequency dependence of Ru provided T is replaced by uj. The backscattering term gi 
(neglected here) could possibly give rise to multiplicative logarithmic corrections to 



the power law in Eq. (5.24i IIT6ll . The sign of the dimensionless prefactor A can only 
be determined through a complete evaluation of {Kx] Ky) in Eq. ( |5.11| , and is for the 
time being unknown. The available experimental data are consistent with Eq. (|5.24[) if 



one assumes that A is negative (see below), as illustrated in Fig. 5.7 



Eq. (5.241 would imply that in the non-interacting limit Kp ^ 1 (52 ^ 0) the cor- 
rection to the Hall coefficient behaves as \og{T /W). In order to check this prediction 
we have evaluated the correlator in Eq. ( |5.11| ) for §2 = 0. This is done in the following 
section. 



5.2.3 i?H in the presence of ukmlapp without forward scattering 



Here we provide the calculation of i?H to leading order in but in the absence of for- 
ward scattering, g2 = 0. Using Eqs ( |5.4| ) and ( |5.8[ ) we can express the zero-frequency 
Hall coefficient in terms of R^^ and Re[Ma;y(iO"'")]. We then perform a Kramers- 
Kronig transform, insert the free-fermion values of the diamagnetic susceptibilities, 
Xa;(0) = — 2e^fF/(vrajy) and Xy(0) = — 4e^t5_^y/(7rfF), and use Eq. ( 3.46| to arrive 
at 



iiH(O) 



1 + 



1 



^Im ( i{Kx] Ky/o 



(5.25) 



where {K^; Ky)^ is to be evaluated to first order in B. From Eq. (5.10 1 one sees that 
{Kx] Ky)Q involves 8 fermion fields and can be represented by diagrams like the one 
displayed in Fig. |5.8| There are 32 different diagrams, but all of them can be expressed 
in terms of only one function A{i^l, B), whose expression is given by the diagram in 
Fig. 5.8 This is done in Appendix |A.3| We thus obtain, 

i?H(0) =R'kU- ^ /" ^ Im U'(a; + iO+) - A{-uj - iO+)l 1 (5.26) 
[ ea J J 

where A'{iQ) = dA{iQ,, B) /dB\B=o and we have pulled all prefactors fromEq. (5. 10 1, 
as well as a factor t± from the diagram, out of the definition of function A. The explicit 
expression of A' is (see Appendix A.3 1 



A' {in) 



(27r) 



dkidk2dq 



d^+{h) 1 



dki 



1 



1 



g3 z — ^ III,-. 

1 



W2 - i+{^2) iv'i - i-{k2 - q) iy\ + 1^2 - i^s + i^, - ^-{ki + q) ' 



(5.27) 



The frequency summations in Eq. (5.27) are elementary, and the various momen- 
tum integrals can also be evaluated analytically to first order in a, yielding (see Ap- 
pendix A.3 1 

Rh{0) 



r9j 



1 

16 



93 y 



doj (/3a;/4)2 -sinh2(/3a;/4) 
uj tanh(/3a;/4)sinh2(/3u;/4) 



(5.28) 
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Figure 5.8: Example of a diagram appearing in Eq. ( |5.1l| at first order in t± and for 
52 = 0. The full (dashed) lines correspond to free right (left) moving fermions, j is the 
chain index, and the arrows represent up and down spins. The magnetic field increases 
the momentum of the electron by dk — eBay. 



The remaining energy integral is divergent and must be regularized. Cutting the inte- 
gral at the bandwidth W and assuming T <^ W we obtain the asymptotic behavior 
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1 + 



1 



93 Y 



log 



wj 



(5.29) 



consistent with Eq. ( |5.24 ). For non-interacting electrons, though, we see that the rel- 
ative correction induced by the 1/2-filling umklapp is positive at T < W. Since all 



properties are analytic in Kp, we can also deduce from Eqs (5.24i and (5.29 1 that A 



tends to [24(1 - Kp)]~^ in the limit Kp ^ 1. Note that Eq. ( |5.29| would also apply to 
models in which g2 
when 53 < 52- 



(73, such as the Hubbard model, while Eq. ( |5.24 ) is valid only 



5.2.4 Discussion and perspectives 



The result of Eq. (5.24i shows that in 1/2-filled quasi ID systems the umklapp pro- 



cesses induce a correction to the free-fermion value (band value R^) of the Hall coef- 
ficient Ru, which depends on temperature as a power- law with an exponent depending 
on interactions. At high temperatures or frequencies, Ru approaches the band value 
as shown in Fig. |5.7[ implying that any fitting of experimental data must be done with 
respect to the value of Ru at high temperature or frequency. 

To study the range of validity of our result, one must consider that at low tempera- 
ture the quasi- ID systems generally enter either in an insulating state characterized by 
a Mott gap A, or in a coherent two- or three-dimensional phase below a temperature 



T* controlled by t±, as explained in Sec. 2.4.3 In either case our model of weakly 



coupled LL is no longer valid, as illustrated in Fig. 5.7 The variations of Ru below 
max(T* , A) can be very pronounced, and depend strongly on the details of the materi- 
als. When the ground state is insulating, for instance, R}i{T) is expected to go through 
a minimum and diverge like e^/-^ as T ^ 0, reflecting the exponentially small carrier 
density. Other behaviors, such as a change of sign due to the formation of an ordered 



state or nesting in the FL regime [49], can also occur. The validity of Eq. (5.24) is 
therefore limited to the LL domain: max(T*, A) < max(r, u>) <C W. 
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For the case A > T*, we estimate the change of i?H with respect to at the 
crossover scale A, for a system with <^ U, where U is the Coulomb repulsion. The 
umklapp-induced Mott gap in 1/2-filled systems is given by A/1^ ~ [53/(^^f)]^ with 
X = [2(1 — Kp)]^^ 1 16 1. We thus find that the largest correction is ~ [53/('?rt'F)] 2 and 
has a universal exponent. On the other hand, Ru approaches the asymptotic value Ry^ 
quite slowly, and according to Eq. (5.24i a correction of ~ [93 / {t^vp)]"^ still exists at 
temperatures comparable to the bandwidth. 

The available Hall data in the TM family and in the geometry of the present anal- 
ysis Il42l l44ll show a weak correction to the free fermion value which depends on 
temperature. Some attempts to fit this behavior to a power law have been reported 
(see Fig. 5.4 1 p^. However the analysis was performed by fitting Ru{T) to a power 
law starting at zero temperature. As explained above, the proper way to analyze the 
Hall effect in such quasi- ID systems is to fit the deviations from the band value start- 
ing from the high temperature limit. It would be interesting to check whether a new 
analysis of the data would provide good agreement with our results. However in these 
compounds both 1 / 4-filling and 1 / 2-filling umklapp processes are present. For the 
longitudinal transport, the 1 /4-filling contribution dominates lITTl . For the Hall effect, 
the analysis in the presence of 1/ 4-filling umklapp is considerably more involved, but 
a crude evaluation of the scaling properties of the corresponding memory matrix gives 
also a weak power-law correction with an exponent 2 — 16Kp + {Kp + Kp^)/2, and 
thus similar effects, regardless of the dominant umklapp. The observed data is thus 
consistent with the expected corrections coming from LL behavior. However more 
work, both experimental and theoretical, is needed for the TM family because of this 
additional complication. 



Our result Eq. (5.24i is however directly relevant for 1/2-filled organic conduc- 
tors such as (TTM-TTP)l3 and (DMTSA)BF4 1351. Hall measurements for these com- 
pounds still remain to be performed. Comparison of the Hall effect in these compounds 
with the one in 1/ 4-filled non-dimerized systems ll50l l29l for which only 1 / 4-filling 
umklapp is present, could also help in understanding the dominant processes for the 
TM family. 



The other type of Hall measurements, shown in Fig. 5.5 were done in a different 
geometry from the one used in the present theoretical work. In this case the current 
flows along the least conducting direction c and the magnetic field is applied parallel 
to a, as explained before. This implies a cunent flowing along a direction which is 
incoherent at high temperature (T > T*) and thus, the results for R^ are expected 
to differ from Eq ( 5.24| ). To describe the system in this geometry, the Hamiltonian 
must have two coupling terms corresponding to the hopping along h and c directions 
and one must choose a gauge for the magnetic field that determines the Peierl's phase 
appearing in the Hamiltonian, as done for the model studied here. The longitudinal and 
transverse currents will have the form of Eq. (5.6b 1 (with a different Peierl's phase), 
each one with the respective hopping amplitude ti^ or tc- The diamagnetic terms will 
be both similar to Eq. (5.5b I, with replaced by th or tc and ay replaced by the 
respective lattice parameter. Since in this geometry the magnetic field goes along the 
one-dimensional chains, the question is whether or not signatures of Luttinger liquid 
behavior will appear on i?H- With the information given until now we can compute 
the first term of Eq. (|5.4|) (the high-frequency i?H of Eq. (|4. 17|)). For this we must 
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calculate the commutator between the currents, obtaining 



Rh[u; ^oo) = h-, (5.30) 

where a and /3 are chain indexes, and ab and Oc are the lattice parameters along b and 



c directions, respectively. It is evident that the average in the numerator of Eq. (5.30 1 



must be expanded to first order in 4 and tc in order to obtain a nonzero result. Im- 



mediately, many questions on the derivation of expression (5.30 1 appear: can we treat 
tc perturbatively even if the current goes along c?, is the ID Hamiltonian term nec- 
essary or just a Hubbard Hamiltonian will be sufficient?, should one consider the ID 
properties of the chains to obtain the averages, even if they do not play a crucial role 
in this geometry?, is the above expression temperature independent?, etc. These and 
many other questions must be answered in order to understand the Hall effect in the 
geometry consider by Mihaly et al. in Ref . H3l . In this work we did not solve this 
problem, but we hope to do it in the near future. 



5.2.5 Conclusions 

In this chapter we have accomplished a theoretical study of the Hall effect in a system 
made of weakly coupled 1 /2-filled chains, in the presence of umklapp scattering. We 
obtained a Hall coefficient i?H given by the free-fermion value (band value plus 
a correction term with a power-law dependence on temperature (or frequency), due 
to the presence of umklapp scattering. This power-law is a characteristic behavior of 
Luttinger liquids, where the exponent depends on the interaction parameters. The Hall 
constant was also computed for the system without forward scattering, resulting in a 
logarithmic dependence on T (or uj), in agreement with the zero interaction limit of 
the power-law. The results are not directly applicable to the Hall data in the quasi 
ID organic conductors, reviewed at the beginning of this chapter, because these are 
1/2- and 1/4-filled compounds, but they allowed us to reach the following conclusion: 
signatures of LL behavior (power-law dependence) are expected to appear, at high 
temperatures, in Hall measurements made in the geometry considered here. 
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CHAPTER 6 



Hall effect on the strongly correlated 2D triangular lattice 



As we have seen in the preceding chapters, the interpretation of the Hall effect in 
strongly correlated systems can result in a very complicated task. Interactions can 
have a large effect in the Hall resistivity, and this effect seams to increase as the di- 
mensionality of the system decreases. Furthermore, the understanding of these effects 
is crucial for the investigation of transport properties in strongly correlated systems, in 
particular for the Hall effect. 

We now know that there exist a variety of strongly correlated systems with differ- 
ent geometries. Among these, the triangular lattice exhibits a unique property: it has 
the smallest possible closed loop with an odd number of steps (namely 3), as explained 
in Sec. |2.3| Anderson proposed that the model could have a spin-liquid ground state 
at commensurate fillings such as one electron per site [51]. These peculiarities make 
the triangular lattice a very interesting system for the investigation of the Hall effect. 
In particular, important differences between the Hall effect in the square and triangu- 
lar lattices were pointed out by Shastry et al. in Ref. ETl . Since we have already 
studied the Hall effect in a quasi ID system, which can be also considered as a highly 
anisotropic square lattice, the investigation of the same phenomena in the triangular 
lattice can be a helpful way to understand the relation between the Hall resistivity and 
the geometry of the underlying lattice in two-dimensional systems. 

There exist a variety of compounds with structures resembling the triangular lat- 
tice. Among these the C0O2 layered compounds where the recent discovery of super- 
conductivity in the hydrated Na^Co02 [52] have motivated a large number of works. 
These materials, also known as cobaltates, are good realizations of an isotropic 2D tri- 
angular lattice and have been extensively investigated, both experimentally 153115411551 
and theoretically fSSj |57l ESl , in the past years. There are also organic conductors of 
the BEDT (bis(ethylenedithio)) family [,591 where one finds various structures resem- 
bling the anisotropic triangular system. This anisotropy, together with the dimerization 
of the molecules present in some compounds of the BEDT family, make these materi- 
als much more complicated to describe from the theoretical point of view. That is why 
we will concentrate here exclusively in the Naa;Co02 compound. 
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Figure 6.1: Molecular structure of 7-Na^Co02 phase with a space group symmetry of 
PQ^/mmc and lattice constants a — 2.84 A and c = 10.81 A. The C0O2 layers are inter- 
calated with insulating layers of Na+ ions. Nai and Na2 differ from their crystallographic 
positions: Nai is situated on the vertical crossing Co atoms above and below, and Na2 is 
slightly displaced with respect to the vertical crossing the center of the triangle formed by 
the Co atoms. From Ref. I16O1I 



In the following we make a review of the most important properties of Na3;Co02, 
which is the material considered in this work for the application of our theoretical 
study of the Hall effect on the two-dimensional triangular lattice. 



6.1 A triangular lattice compound: Naa;Co02 

The crystal structure of Naa;Co02 consist on two-dimensional C0O2 layers of edge- 
sharing tilted octahedra. Each octahedra is composed by a Cobalt ion surrounded by 
six Oxygen atoms at the vertices CoOg- Within each C0O2 layer, the Co ions occupy 
the sites of a two-dimensional triangular lattice. The C0O2 layers are separated by 
insulating layers of Na+ ions. There exist four phases of Na2^Co02, with slightly 
different structures, called a, a' , j3 and 7. They differ by the stacking order of C0O2 
layers and Na-0 environments ||6T]| . In this work we will focus only on the 7 phase 



because it is the one used in the Hall measurements. Fig. 6.1 shows the structure 
of the 7-Na2,.Co02, which has an hexagonal structure with a space group symmetry 
of PQ^/mmc and lattice constants a = 2.84 A and c = 10.81 A. Band-structure 
calculations f56l show that the O 2p orbital states lie far below the Co 3(i states and 
the chemical potential falls within the band formed from t2g states in Co. Hence the 
electrons donated by the Na ions are distributed among the Co ions, a fraction {6) of 
which are in the Co^+ state in which S = 1/2, while the rest (1 — 5) are Co'^^ with 
5 = 0. The elementary charge-transport process is the hopping of a hole from Co^"*" 
to Co^+. A large on-site repulsion excludes double occupancy of a site by the holes. 



Different ordered states appear as a function of the doping x of Na ions. Fig. 6.2 
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Figure 6.2: Phase Diagram of Naa;Co02. The different order states are found changing 
the doping of sodium atoms (x). From Ref. lfT3l 



shows the phase diagram for Naa:Co02. A paramagnetic metal is found forx<^ and 
for X > I the system behaves as a Curie -Weiss metal [621 . These two metallic states 
are separated by a narrow charge-ordered insulating state at x = ^. This compound 
becomes superconducting below 5K, when intercalated with water forming Na3;Co02 • 
yH20, for J < X < |. For values of the doping above |, a Spin Density Wave metallic 
state appears f62', T3l. Na2.Co02 presents also an unusual enhanced thermopower at 
a; ~ 2/3 which has been recently related to the spin entropy carried by the holes in the 
Curie-Weiss phase Il63l . Nao.7Co02 is the host compound from which Na^ is varied 
to achieve superconductivity in the hydrated compound. 

The effective Hubbard interaction U have been estimated in band-structure calcu- 
lations |[56l to be of the order of [/ ~ 5-8 eV and a Fermi surface with hexagonal 
character in agreement with photoemission measurements ll55l 1541 . The dispersion 
behavior seen in photoemission data is consistent with a negative sign of the single- 
particle hopping of the order of t = 10 it 2 meV ||55]| . Therefore, the bandwidth (W) 
is estimated between 70 to 100 meV making this system a real strongly correlated 
one with U » W. Several Hall measurements have been undertaken in Nao.7Co02 
||63 , 641. The anomalous linear increase of the dc Hall coefficient and a recent infrared 
Hall measurement 1641 have motivated recent theoretical works ll57l l58l l65l l66l with 
the aim of investigating the role of correlations in the Hall effect, but many questions 



remain open. Fig. 6.3 shows the Hall constant measured by Choi et. al. Ii64il at an 



infrared frequency and at w = 0. 

This chapter is the object of publication |[67]| . In the following we make a theoret- 
ical study of the Hall effect in a 2D triangular lattice where electrons interact via an 
onsite Coulomb repulsion U. We calculate Ru in the high frequency limit studied in 
chapter [4j and we cover the whole range of interaction values using several approxi- 
mation schemes. 
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Figure 6.3: Temperature dependence of the Hall constant measured in Nao.7Co02 by 
Choi et. al. |64|. (Left) Infrared Hall constant measured at — 1100 cm^^. The upper 
(lower) curve correspond to the real (imaginary) part of R\i{Lo) . (Right) dc Hall constant 
(w = 0). The scale on the left differs from the one shown in Ref [64 1 because it was 
corrected by the authors. 



6.2 Model and methods 

For this theoretical study, we consider an anisotropic triangular lattice with nearest- 
neighbor hopping amplitudes t and t' and an on-site Hubbard interaction U, with the 
structure sketched in Fig. |6.4| The system is described by the Hubbard model on the 



triangular lattice studied in Sec. 2.3 



W = - ^ tijc\^c-^ + UY^ Ui^rii^ (6. 1) 

where (cq) is the creation (annihilation) fermion operator, is the fermionic num- 
ber operator and {ij) are nearest-neighboring sites. The dispersion relation for this 
model (see Sec |2.3|) is 



Efc = —2tcos{kxa) — 4t' cos{kxa/2) cos{kyaV3/2). (6.2) 



We assume that a current I flows along the x axis and a dc magnetic field B is 
applied along z, hence a Hall voltage develops along the y axis (see Fig. 6.4). We 
use the vector potential A = + where for the magnetic part we choose, 

as in the previous chapter, the Landau gauge A'"^^ = Bxy, and A^^ describes the 
electric field. The coupling between the lattice fermions and the electromagnetic field 
induces a Peierls phase in the hopping amplitudes which change according to Uj — > 
tij exp(— ie fj' A • dl). 
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Figure 6.4: Top: Two-dimensional triangular lattice, a is the lattice parameter, t and 
t' are the hopping amplitudes for bonds along the x direction and for ±60° bonds, re- 
spectively. The unit-cell area is 5 = a^\/3/2. The current / flows along the x axis, the 
magnetic field B is applied along the z axis, and the Hall voltage is measured along the 
y axis. Bottom: Non-interacting density of states of the model in the cases t = t' = —1 



and < = 2t' = — 1, presented in Sec. 2.3 



We use Eqs. (3.5 1 and (3.161 for the total current operator, = / (irj^(r), and the 
diamagnetic susceptibilities XAt(O), respectively. Performing the functional derivatives 
we find for the currents 



Jr = ea 



2*5Zi^fc^ sin(/c^a) + 

fccr 



(6.3a) 



J 1] 



-eaVst' ( ^ + f ) {^clc,^,^e^'^^''^ +h.c.) , (6.3b) 



fccr 



where we have defined the vector rj = (ry,0) with rj = \/^eBa/2. The diamagnetic 
susceptibilities are: 



Xy{0) 



4e^ 1 \ - . . ( K.a 



|^2t cos{kxa) + t' cos 



cos ( kyVS 



k^a 

N ^cos^— icos 
fc 



(nfc) 
(6.4a) 
(6.4b) 



As shown in Sec. 4.2 it is possible to rewrite Ru as a high-frequency series where the 
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infinite-frequency limit reads 



oo 



lim 



{[JxiJy]) \ 



BNSxMXy{^)) 



(6.5) 



The remaining contributions are expressed in terms of a memory matrix (see Sec. |4.3| ). 
i^Hl"^ — > oo) is expected to provide tiie dominant contribution at any finite frequency. 
The memory matrix formalism allows in principle to go beyond the infinite frequency 
approximation and compute corrections at finite frequency ll28l l46l . It leads, in par- 
ticular, to corrections due to interactions that vanish identically if ?7 = 0. These 
corrections do not affect the sign of Rn which should be entirely determined by 
R\{{u} — > oo). In the following we shall consider only the infinite-frequency con- 
tribution to Rn, Eq. (6.5 1, and adopt the notation i?H('j-' oo) = i?H. 

Strictly speaking, our results are valid provided the probing frequency is larger 
than any other energy scale in the system, uj > max{f7, t, T}. The last two conditions, 
io > max{t, T}, are easily fulfilled experimentally in known triangular compounds, 
while the condition lo > U is more problematic. However, as we will discuss in Sec. 
6.41 in certain limits our results coincide with those obtained in Ref. ll57l under the 



opposite assumption lu <^ U, showing that this condition is not stringent. 



In order to evaluate Eq. ( 6.5 1, we calculate the commutator [J^ , Jy] from Eqs. ( 6.3 1, 



and we use the diamagnetic susceptibilities of Eq. (6.4i to arrive at 



T.kAk{nk) 



e NUk^kink. 



jfJ2k^k{nk)' 



(6.6) 



with 



Ak 

Bk 
Ck 



cos 



cos{kxa) cos ( ky 



+ -{t' /t) cos{kxa) + cos (^kyVSc 



2 cos{kxa) + {t' /t) cos 



cos ( kyVS 



cos I j cos (^kyVs*^^ . 



(6.7) 



As can be seen from Eq. (6.6 1 the high frequency Hall coefficient depends only on the 
distribution function (uk) as well as some geometrical factors. The interaction term in 
Eq. (6.1 1 therefore only influences Ru through its effect on (nfc). Another implication 



of Eq. (6.6 1 is that at low temperature the behavior of Rn can be interpreted in terms 



of an effective carrier concentration, as in the non-interacting case. 



6.3 Results 

In the following we evaluate Ru in the whole domain of interaction values U with 
respect to the bandwidth W = 9\t\ of the system, by using four different approaches: 
exact calculation at ?7 = 0, a perturbative expansion of the self-energy at U < W, 
a local approximation to the self-energy, treated with dynamical mean field theory 
(DMFT) atU >W, and finally the atomic hmit of the self-energy atU :$>W. 
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6.3.1 Non-interacting case 

In the non-interacting case there are various Umits in which we can obtain analytical 
results for (i?H at [/ = 0): at zero temperature and band fillings near n = and n = 
2, and at high temperature T » W. For intermediate fillings and temperatures, we 
compute numerically by performing the sum in Eq. (6.6 1 on a dense 2048 x 2048 



discrete fc-point mesh. Some details on the numerical sum over momentum are given 



in Appendix B.l 



Zero temperature 

Here we restrict for simplicity to the isotropic case t' = t and we set the lattice param- 



eter a = 1. Close to the band edges we can expand the various integrands of Eq. (6.7 1 
and thus perform the k integrals. 

Near the bottom of the band the Fermi surface is made of two nearly circular elec- 
tron pockets around (^,0) and (x' ^^ch pocket we have ■?fc = ^fc — ~ 
St — — where k is the momentum measured from the pocket center, and there- 
fore = ^{3 — fi/t). The corresponding electron density is n = /cp/vr. Writing 
similar expansions of A^, B^, and Cfc close to the pocket center and performing the 
Brillouin zone integrations, we obtain the non-interacting Hall coefficient at low elec- 
tron density: 



rUt = 0) = - 

ne 



1 — + 0{n 



2^ 



(6.8) 



At sufficiently low density we recover, in the above expression, the classical result 
R^ = 1/ne. 

Near the top of the band the Fermi surface is a nearly circular hole pocket centered 
at = (0,0). Close to this point we have ^fc w —6t + ^tk"^ — /x, and therefore 
A^p = |(6+/i/t). The corresponding density is obtained by subtracting the contribution 
of the hole pocket from the maximum density: rih = 2 — kp/2TT. Similarly, for the 
functions Aj^, Bf,, and we have to subtract the contribution of the hole pocket from 
the contribution of the whole Brillouin zone, which turns out to be zero because 

Y^Ak = Y,Bk = Y.^k = 0- (6.9) 



Thus, for low hole densities Uh = 2 — n we find that the non-interacting Hall coefficient 
is given by 



R^{T = 0) = - — 



■h) 



(6.10) 



and as n/j ^ we have Ry^ = — l/n/^e. 

The complete density dependence of R^ calculated numerically at zero tempera- 



ture from Eq. (|6.6|) is displayed in Fig. 6.5 and compared to the limiting cases Eqs. ( 6.8 1 



and ( 6.10| ). It is clear from this figure that the infinite-frequency i?H follows the well- 



known dependence of the dc Hall coefficient iinCf^ = 0) the carrier charge density. 
This indicates a weak frequency dependence of the non-interacting Hall coefficient at 
zero temperature, since the dc result is recovered from the infinite frequency limit of 
i?H- Furthermore this suggests, as we will discuss in more details below, that the fre- 
quency dependence should not be too crucial, even in the presence of interactions, for 
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Figure 6.5: Non-interacting Hall coefficient i?^ at zero temperature as a function of the 
electron density n, for an isotropic triangular lattice with t' = t = —1. The dashed lines 
indicates the classical behavior at low electron and hole carrier densities. 



most band fillings. At [/ = the sign of the Hall coefficient is entirely given by the 



sign of the carriers, and it can be seen from Fig. 6.5 how the sign changes at quarter 
filling when the Fermi energy crosses the van Hove singularity of the DOS, and the 
Fermi surface shape evolves from electron to hole like. 



High temperature 

If T S> t the distribution function (nfc), which reduces to the Fermi distribution at 
U = 0, can be expanded in power of /3 = 1/T. This expansion must be done at 
constant density n, which requires that remains finite as /? — > 0, in other words 
/i ~ T at high temperature. Taking this into account we can deduce the relation 
between and n, exp(— = 2/n — 1, and write the Fermi distribution as 



(nk) 



n 
2 



n(2 



(6.11) 



Due to Eq. ( 6.9 1 the fc-independent terms in Eq. ( 6. 1 1 1 do not contribute to R^, which 



in this case takes the form: 



i?o (T > = -4T 



S 



1 



e n{2 -n) ^J2k Bk^kj^ Efc Ck^k ' 
Performing the Brillouin zone integrations we obtain 



(6.12) 



rUt » t) 



T/t 



1 



n(2 



n] 



2 2 + {t'/tf 



(6.13) 



This result is plotted in Fig. 6.6 together with the numerically calculated full temper- 



ature and density dependence. The most striking feature of Eq. (6.13 1 is the linear 
increase of R^ with T. The same linear behavior was obtained in Ref. ||68J at a; = 0, 
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Figure 6.6: Temperature and electron density dependence of the non-interacting Hall 
coefficient R^, for the isotropic triangular lattice with t = t' = —1. The dashed line 



shows the asymptotic behavior described by Eq. (6.13 1. Inset: Temperature and density 
dependence of the chemical potential /i, illustrating the relation fi ^ T ?A high tempera- 
ture. 



indicating a weak frequency dependence of at high temperature. Our result shows 
that the T-linear dependence of i?H is not due to interactions but to the peculiar topol- 
ogy of the triangular lattice. The sign of at high T is determined by the sign of t, 
irrespective of the density (see Fig. 6.6 1. We attribute this property to the fact that at 
high enough temperature the full band contributes to the Hall effect; hence the sign of 
Rn reflects the dominant nature, electron or hole-like, of the band. As is clear from 
Fig. 6.4 for t < the band is dominantly hole-like, while for t > it is electron-like. 

The relevance of result ( |6. 13| ) is that even without interactions, the Hall coefficient 
has a linear dependence at high temperature due to the geometry of the system, em- 
phasizing the peculiarity of the triangular lattice. By contrast, on the square lattice the 

same analysis yields a T-independent non-interacting R^^ = | ^ — ^1^2-11) 
temperature. 



6.3.2 Weakly interacting regime 

When interactions are present, the distribution function (rife) can be expressed in terms 
of the one-electron self-energy Ti{k, iojn) as:||2l 

with ujn = (2n + 1)itT the odd Matsubara frequencies. In the weak coupling regime 
U < W, we evaluate the self-energy using conventional perturbation theory in U and 
we keep only the lowest order contributions of order U^. For a local interaction like 



the Hubbard term in Eq. (6.1 1 there is only one diagram which is drawn in Fig. 6.7 
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Figure 6.7: Brillouin zone average of the real and imaginary parts of the self-energy 
Eq. (6.15 1 at low temperature T — 0.1, calculated using a 64 x 64 fe-point mesh (solid 



lines). The small standard deviations cr (shaded curves) illustrate the weak momentum 
dependence of S(fe, iw„). The dashed lines show the local self-energy resulting from the 
DMFT calculation (see Sec. 6.3.3 1. The density was set to n = 1.54, which is the value 



for Nao.7Co02 (see Sec. 6.4 1. Inset: Feynman diagram corresponding to Eq. ( 6. 15 i 



The standard diagrammatic rules yield the following expression for the self-energy: 

E(fc,iu;„) = > ; (6.15) 



fclfc2 



/(^fc2) [/(^fcl) ~ /(^fc+fci-fc2)] ~ /('^fci)/(-Cfc+fci-fc2) 



where /(^fc) is the Fermi distribution function. 



The numerical evaluation of Eq. ( 6. 15 1 is demanding due to the double momentum 
integration. This is particularly time consuming because our calculations are done at 
fixed density, and thus require to calculate Ti{k,iuJn) many times in order to deter- 
mine the chemical potential. However it turns out that the momentum dependence of 
iujn) in Eq. (6.15 I is weak. This is illustrated in Fig. 6.7 where we plot the Bril- 



louin zone average of the self-energy, S(ia;„), as well as its standard deviation. The 
weak momentum dependence allows us to compute iujn) on a coarse (typically 



16 X 16) fc-point mesh, and then to interpolate using splines (see Appendix B.l I onto 
a dense mesh for the evaluation of (n/,) and eventually R^. The Matsubara sum in 
Eq. ( 6.14| ) also requires special attention: the formal regularization of the divergence 
through the exponential factor is not suitable for a numerical evaluation of the sum. 



We therefore rewrite Eq. (6.141 as 



{nk) 



1 1 ^ 



1 



lUJr, 



(6.16) 



The ujn sum is now convergent and can be efficiently calculated via the truncation at 
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some large frequency and the analytical evaluation of the remaining terms using an 
asymptotic expansion of the self-energy. 

The i?H resulting from perturbation theory are valid in the regime U < W <^ uj, 
with W = 9\t\ the bandwidth of the system. As already anticipated the effect of a 
small U on the distribution (rifc) is a subtle broadening, and as a result the dependence 
of i?H on U is very weak at low U. Fig. 6.8 provides an illustration of this weak 
dependence. As a consequence the non-interacting results of Sec. 6.3.1 are expected 
to give a fairly good account of the Hall effect for an interaction strength smaller than 
the bandwidth W. 

An important observation which we can make from our perturbative calculations 
is that the momentum dependence of the self-energy is very small, i.e. the self-energy 
is almost local in real space. This suggests to approach the strong-coupling regime 
U >W by assuming that the self-energy is exactly local. In the following section we 
study such local approximations to the self-energy, and we compare them to the result 
of the perturbation theory. 



6.3.3 Strongly interacting regime 

Assuming that the self-energy is local in first approximation, we investigate here two 
models for and their implications for the Hall coefficient i?H. The first ap- 

proach is based on the dynamical mean field theory (DMFT) 1 6 1 and requires to solve 
a difficult self-consistent quantum impurity problem. Due to numerical difficulties this 
method cannot be pushed to very high interactions and/or very low temperature. Our 
second approach is based on a simple analytical form for which is expected 

to be valid at f/ » W, and allows us to express (nfc) analytically in this limit. 



DMFT 

The DMFT approximation provides the exact solution of the problem under the as- 
sumption that the self-energy is local. In this framework the self-energy is expressed 
as: 

T.{iuJn) = GQ^ilLVn) - G~^{i^^n) (6.17) 

where Go is an effective propagator describing the time evolution of the fermions in 
the absence of interaction, and Q is the full propagator, which takes into account the 
local Hubbard interaction. The calculation of Q from a given Qq amounts to solve the 
problem of a quantum impurity embedded in a bath. We do it by means of the quantum 
Monte Carlo Hirsh-Fye algorithm ll69ll as described in Ref. f^l (see Appendix |B.2|). 



From the requirement that Q coincides with the local Green's function of the lattice, 



I.e. 



one can deduce the self-consistency condition 

QQ^{iuJn) = 1/D[iujn - + S(ia;„), (6.19) 

where D{z) = f d^D{^)/{z - is the Hilbert transform of the DOS D{^) cor- 



responding to the triangular lattice and shown in Fig. 6.4 Once the self-consistent 
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Figure 6.8: Evolution of the high-frequency Hall coefficient with U calculated using 
different approximations at T = |<| and n = 1.54, for an isotropic triangular lattice with 

t^t' = -1. 



is injected in Eq. ( 6. 16 1 to 



QQ{iLUn) is obtained, the corresponding self-energy T,{iLiir, 
compute i^H- 

In Fig. |6.7| we compare the DMFT self-energy with the Brillouin zone average of 
the perturbative expression Eq. ( |6.15 1, both calculated at U = 4. It can be seen that 
the frequency dependence and the order of magnitude of the two quantities are very 
similar, suggesting that the self-energy is dominated by the C/^ term and therefore the 
domain of validity of the perturbation theory is not limited to very small U. On the 
other hand it shows that the DMFT, although it is not a perturbative approach, provides 
a smooth transition from the weak to the strong-coupling regimes. This is further illus- 



trated in Fig. 6.8 where we see that the values of i?H calculated by perturbation theory 
and DMFT coincide up to [7 f» 4|t|. At not too low temperature the DMFT calcula- 
tion is reliable up to interaction strengths comparable to the bandwidth W. We have 
performed DMFT calculations at [7 > W, but since these results could be affected 
by systematic statistical errors in the Monte-Carlo summation, they are not shown in 
Fig. Appendix [R!2|). At f/ > it is expected that the DMFT resuh approaches 
the atomic limit in which accurate calculations can be performed, as discussed in the 
next paragraph. 



Atomic limit 

In the limit of very strong interactions U ^ W we. assume that the self-energy ap- 



proaches its atomic limit given by the expression: (see Appendix B.3 1: 



nU 



+ 



n/2(l -n/2)C/2 



(6.20) 



2 ' ioon + fJ-at - - n/2)U 
with /Uat the chemical potential in the atomic limit, not to be confused with the lattice 



chemical potential /i. Using this expression in Eq. (6.14) it is possible to evaluate 
analytically the sum over Matsubara frequencies and thus to obtain a closed expression 
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Figure 6.9: Temperature and density dependence of the high frequency Hall coefficient 
calculated in the atomic limit at U — 500 The dashed line shows the asymptotic 
and S — I 



behavior, Eq. (6.21 



1|. For densities close to half-filling (n — 1.05), 



RniT) deviates from the asymptotic behavior (see text). 



for (nfc) (Appendix B.3 1. In Fig. 6.8 we sliow tlie Hall coefficient calculated witli tlie 



atomic limit of the self-energy in tiie whole range of interaction values. i?H obviously 
converges to the non-interacting limit at low U since the atomic self-energy vanishes 
at [/ = 0, and provides a good interpolation between the weak and the strong-coupling 
regimes. At intermediate values U ^ W the atomic limit is not reliable, although it 



gives the correct order of magnitude for i?H. Fig. 6.8 also shows that saturates at 
sufficiently large U. 



In Fig. 6.9 we display the temperature and density dependence of i?H at ?7 = 
500|t|, which is a typical value for the cobaltate compounds as discussed in the next 
section. We have selected four densities corresponding to the bottom and top of the 
lower and upper Hubbard bands (see also Fig. 6.11 below). Like for U = we find 
a T-linear increase of Ra at T > W. Due to the Mott gap, however, the density 
dependence of the slope is not the same as for f/ = 0. The slope can be obtained 
explicitly by sending f/ to +oo and performing the high-temperature expansion as in 
Sec. l6.3.1I The result is 



T}U = 00 



(r>t) 



T/t 



1 



e 6(1-6) 4 2 + (tVt)2' 



(6.21) 



1-1 



very similar to Eq. (6.13 1 except that the slope oc [4(^(1 — (5)] ^ replaces [2n(2 — n) 
where 6 = In — 1| measures the departure from half-filling. The U = oo result 



of Eq. ( |6.21| ) is displayed in Fig. |6.9[ and correctly describes our high-temperature 
results at [/ = 500|t|. The differences observed at n = 1.05 in Fig. 6.9 reflect the 



fact that close to half-filling the slope of the high temperature i?H depends strongly on 
the interaction and is not saturated even at f7 = 500|t| (see also Fig. 7). Away from 



half-filling the U dependence of the slope is weaker, and Eq. (6.21 1 is valid for lower 
interaction values. 
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Figure 6.10: Density dependence of the Hall coefficient within the T-linear regime at 
U ^ and in the atomic limit atU^ \t\ (solid lines), compared to the result of the t-J 
model (Ref. Il57l . dashed line). 



6.4 Discussion and perspectives 

The various approximations presented above allow us to calculate the Hall coefficient 
on the triangular lattice for all interactions strengths U and all temperatures T. The 
main limitation of our approach, in view of a comparison with experimental systems, 
is that our results are in principle valid in the limit W,U <^ to, because they are based 
on a high-frequency expansion. The first criterion, W <C w, is not too difficult to 
satisfy for realistic compounds if the measurement of the Hall effect is performed at 
optical frequencies. The second criterion, [/ <C seems more problematic since 
interaction strengths can be as large as several electron volts, at the upper edge of the 



mid-ultraviolet frequency domain. However, we have seen (Fig. 6.5 1 that at [/ = and 
T = the Hall coefficient calculated at w = oo coincides with the w = dc value, 
and at f7 = and T » t, we obtained the a; = results of Ref. [68|. All this suggests 
that the frequency dependence of Ru is weak in the non-interacting case. 

At the other extreme of the parameter space, U = oo and T » W, we can com- 
pare the result of the atomic limit approximation with the result of the t-J model ll57l . 
In the latter model U is considered infinite from the outset, so that the high-frequency 
and high temperature expansion of Ref. |57| is in fact valid at frequencies uj <U . We 



plot in Fig. 6.10 the density dependence of obtained in both models at f/ = cxd and 
T > W. The small quantitative difference between the atomic limit at f/ = cxo and the 
t-J model shows that these two ways of treating the U = oo limit are not equivalent: 
they differ, in particular, in the renormalization of the kinetic energy by the interaction. 
However the two models give away from half filling very similar behaviors. This re- 
inforces the idea that the frequency dependence of Ru is weak. Exact diagonalization 
on small clusters also indicate such a weak frequency dependence |[58l . This strongly 
suggests that our results could also be valid at lo < U, and therefore be relevant to 
interpret experiments performed in this regime. The atomic-limit approach has the 
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Figure 6.11: (a) Density dependence of the Hall coefficient at T = and U = 500\t\ 
calculated in the atomic limit approximation, (b) Density of states for various electron 
densities, showing the lower (LHB) and upper (UHB) Hubbard bands. The shaded regions 
indicate the occupied states and the position of the chemical potential. 



advantage to give access to the full temperature dependence (Fig. 6.9 1 as well as the 
[/-dependence as shown in Fig. 6.10 while the calculation of Ref. lISTIl is valid at 
[/ = oo and T » VF. 

The evolution of Ru with temperature is of particular interest since it is most easily 
probed experimentally. A linear increase of Ru with temperature, without saturation 
at high T, was reported in Ref. [57] for the t-J model. Our results show that the 
Coulomb interaction is not responsible for this effect which is also present at [/ = 
(Fig. 6.6 1 and is therefore a consequence of the peculiar geometry of the triangular 
lattice. However the interaction controls the density dependence of the slope which 
changes smoothly from [2n(2 — n)]~^ at f/ = to [4(5(1 — S)]~^ at J7 = oo. This is 
further corroborated in Fig. |6.10| 

The sign of Ru turns out to be independent of n and U at high temperature, unlike 
in the square lattice where Ru changes sign at n = 1. The situation is different at 
T = 0. In the non-interacting case Ru changes sign at quarter filling and can be simply 
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interpreted in terms of the carrier density (Fig. 6.5 1. We have also investigated the 
T = density dependence of i?H at large U as shown in Fig. 6.11 The interpretation in 
terms of the carrier density remains qualitatively valid, provided one takes into account 
the splitting of the DOS into the lower and upper Hubbard bands. These two bands 



are displayed in Fig. 6.113, where it can also be seen that the DOS keeps qualitatively 
the same shape as for [7 = 0, but the width of each band varies strongly with the 
density n. Due to this band renormalization the sign change of Ru at n < 1 does 



not occur at quarter filling, but a little below. Comparing Fig. 6.11 1 with Fig. 6.10 
one easily understands why the temperature dependence of Ru is more pronounced 
slightly above n = and n = 1, where Ru changes from negative at T = to positive 
at high T, than slightly below n = 1 and n = 2 where its stays positive (see also 



Fig. 6.91 



Let us now discuss the application of our theoretical results to the Hall measure- 
ments performed at finite frequency by Choi et al. f64l on the cobaltate Nao.7Co02 



(see Fig. 6.3 1. As mentioned before, ARPES measurements [55] indicate that the trian- 
gular lattice is isotropic with an estimated hopping amplitude of t = —10 meV and an 
effective Hubbard energy [7 ~ 5 eV. From the radius of the Fermi-surface hole pocket 
observed in ARPES, fcp = 0.65 it 0.1 A^^, we deduce an electron density n = 1.54. 
Choi et al. measured the temperature dependence of both the dc and ac Hall coeffi- 
cients up to room temperature (see Fig. |6.3| ). The ac measurement was performed at 
u> = 1100 cm^^ 12|t|. The experimental conditions thus satisfy T,W < uj <^ U . 



In Fig. 6. 12 we plotted together the experimental and theoretical curves. As can be 
seen in the figure, there is a factor of 10 between the measured Ra at finite and zero 
frequency. The order of magnitude of our theoretical curve, obtained with the param- 
eters measured with ARPES, is in agreement with the dc Hall data and not with the 
ac Hall data, as expected. We do not have, at the moment, any clear explanation for 
this discrepancy. One would have to extend the theoretical approach in order to cover 
the domain of intermediate frequencies and more infrared measurements are needed 
in order to assure that the experimental data is correct. In any case, we compare our 
results with the dc Hall data at high temperatures and the results are discussed below. 

The behavior of the dc i?H above T = 250 K is consistent with the linear increase 
predicted by the various theoretical models. By adjusting these models on the dc 



experimental data at high temperature (dotted line in Fig. 6. 12 1 we obtain independent 



determinations of the hopping amplitude t, namely t = —7.4 meV using the atomic 



limit model Eq. ( 6.21 1 and t = — 5.7 meV using the t-J model. This values are in good 
agreement with the ARPES results. We note, however, that there are discrepancies 
between different sets of experimental data |[64ll63]| . 

The organic conductors of the BEDT family present several compounds with an 
anisotropic triangular structure. Unfortunately we are not aware of any measurements 
of the ac Hall effect which we could compare to our calculations, although measure- 
ments have been done at zero frequency in these materials. f/Ol ITTl 



6.5 Conclusions 

The theoretical high-frequency Hall coefficient in the two-dimensional triangular lat- 
tice exhibits two different characteristic behaviors at low and high temperatures: near 
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Figure 6.12: Comparison of the Hall coefficient in Nao.7Co02 as measured by Choi 
et al. |64| at w = (empty circles) and lo — 1100 cm^^ (full circles) with available 
theoretical models. The blue solid line is the high-frequency result in the atomic limit 
and the dashed line shows the high-temperature result in the t-J model 157 1. Theoretical 
curves are calculated alt — t' — —10 meV, [/ = 5 eV, and n — 1.54. The dotted line 
shows the best linear fit of the experimental data at high temperature (see text). 




T = 0, i?H resembles the classical dc Hall coefficient 1/qn* where q and n* are the 
carrier charge and density, respectively; at temperatures higher than the bandwidth, i?H 
shows a remarkable T- linear behavior with a density- and interaction-dependent slope. 
These conclusions apply provided the probing frequency is larger that the other energy 
scales of the problem, and that the electron self-energy remains essentially local for 
strong interactions. 

Although we argued that the frequency dependence of i?H is probably weak, it is 
clear that for understanding the anomalously large i?H(w) measured experimentally in 
Nao,7Co02 in the mid-infrared range, one would have to extend the approach in order 
to cover the domain of intermediate frequencies. Concerning the self-energy, we do 
not exclude that this quantity could present a non-negligible momentum dependence 
for strong interactions. This could affect the Hall coefficient, especially at low tem- 
perature. Such a momentum dependence is indeed expected for ordered ground states 
at and close to half-filling due to antiferromagnetism. However, the influence of a 
momentum-dependent self-energy would be reduced on the Hall effect, because the 
expression Eq. ( |6.6| ) for the Hall coefficient averages the distribution function over the 
Brillouin zone. In addition we do not expect a strong momentum dependence far from 
half-filling, as in the case of the n = 1.54 cobaltate considered here. 

Another possibility is that the simple one -band model considered in this study 
would not suffice to capture the detailed properties of the materials 11721 . Experiments 
conducted as a function of u), as well as measurements of other materials with a trian- 
gular structure, would be very helpful to elucidate the peculiarities of the Hall effect 
in triangular compounds. 
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CHAPTER 7 



Conclusions 



This work was devoted to the investigation of the Hall effect in two different theoreti- 
cal models of strongly correlated systems: a system made of weakly coupled Luttinger 
liquids and the two-dimensional triangular lattice. In order to have the necessary theo- 
retical tools to tackle these problems, we started with a review on the main properties 
of systems with strong interactions, focusing on two-dimensional and low-dimensional 
systems (ID and quasi ID). Because the Hall effect is first of all a transport phe- 
nomenon, we reviewed the formalisms existing in the literature to treat transport prop- 
erties in strongly interacting systems, and we dedicated a whole chapter to the specific 
treatment of the Hall effect. The conclusions corresponding to each of the models 
studied are presented in what follows. 

The study of the Hall effect in a quasi ID system was in part motivated by various 
Hall measurements made on quasi ID organic conductors. This work consisted on 
the application of the memory matrix formaUsm to obtain the Hall constant (Rh) in 
weakly coupled 1/2-filled chains, in the presence of umklapp scattering. The geome- 
try of the model was chosen in order to have the current flowing along the ID chains. 
We computed the temperature and frequency dependence of Ru taking into account 
particle-particle interactions and particle-lattice interactions (umklapp processes). We 
obtained a Hall coefficient i?H given by the free-fermion value (band value R^) plus 
a correction term with a power-law dependence on temperature (or frequency), due 
to the presence of umklapp scattering. These power-law dependencies are signatures 
of Luttinger liquid behavior, where the exponents depend on the interaction parame- 
ters. The Hall coefficient was also computed for the system without forward scattering 
(particle-particle interaction), resulting in a logarithmic dependence on T (or u), in 
agreement with the zero interaction limit of the power-law. Because the organic con- 
ductors are 1/2-filled and 1/4- filled systems at the same time, our theoretical results 
are not directly applicable to explain the available Hall data, but they allowed us to 
make valuable conclusions. Firstly, the proper way to analyze the Hall data, made 
in such quasi- ID systems (in the same geometry of our model), is to fit the deviations 
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from the band value starting from the high temperature limit. And secondly, signatures 
of LL behavior (power-law dependence) are expected to appear, at high temperatures, 
in Hall measurements made in the geometry considered here. A second geometry, 
where the magnetic field is applied along the ID chains, can be considered for this 
system and the question of the appearance of LL signatures in this case remains open. 

The second model studied in this work was the two-dimensional triangular lat- 
tice with an onsite interaction U. This theoretical study was partially motivated by 
Hall experiments made at finite and zero frequency on the sodium cobalt oxide. We 
computed the infinite-frequency Hall constant, where the probing frequency must be 
larger that the other energy scales of the problem, for all interactions strengths U and 
all temperatures T. We obtained an Rii{oj — > oo) with two different characteristic 
behaviors at low and high temperatures: near T = 0, i?H resembles the classical dc 
Hall coefficient 1/qn* where q and n* are the carrier charge and density, respectively; 
at temperatures higher than the bandwidth, presented a T-Unear behavior with a 
density- and interaction-dependent slope. These results were applied to the Hall data 
measured on cobaltates at finite (infrared) and zero frequency. The finite frequency 
data showed an order of magnitude of difference with our theoretical results. To un- 
derstand this discrepancy, one would have to extend the approach in order to cover the 
domain of intermediate frequencies. The behavior of the dc i?H at high temperature 
was consistent with the linear increase predicted theoretically. By adjusting the model 
on the dc experimental data at high temperature we obtain independent determinations 
of the hopping amplitude in accordance with the experimental values. An extension 
to the work done here can be the calculation of the memory matrix term for the Hall 
constant to leading order in U. This would give the behavior of the Hall constant at 
lower frequencies. 

After summarizing the conclusions corresponding to each one of the models stud- 
ied, we must point out that the study of the Hall effect in strongly correlated system is 
a very comphcated problem from the theoretical point of view. Interactions can have 
a large effect on the transport properties. In particular, for the Hall resistivity these 
effects seem to increase when the dimensionaUty of the system decreases. However, 
with this work we have learned that there are feasible ways to treat this problem and 
that each geometry must be treated in a different way. In any case, much more needs 
to be done in order to understand form a theoretical point of view the large number of 
experiments and the real meaning of the Hall resistivity in strongly correlated systems. 
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APPENDIX A 



Appendix for the study of the Hall effect in quasi 1 D 

systems 



A.l Correlator {Kx] Ky) at zero order in a and B 



The operators , Ky and Ti.^ are given in boson representation in Eqs . ( 5 . 1 3 1 and (5.151. Here 



we replace b = +1(— 1) by r = +1(— 1) for right(left) moving fermions in the expression for 
Ky. To obtain the correlator at zero order in the band curvature term and magnetic field, we 
take Eq. (5.12 1 with Ha = and we put _B = in expressions (5.13 1 and (5.15 1, 



{K,;Ky} = - J dre^'^^ J dTi (TrK,{T)Ky{0)^°^nf\n)) . (A.l) 

We denote by (...) the time-ordered correlation function and the superindice (0) means zero 
order in B. All the terms are multiplied by the factor: (2et;F.93/(27ra)^) (iet^^g-^ay/ (27ra)^) 
(t±/2i:a). We will use the letter Ci to denote it in the following calculations. In bosonization 
language, the average gives 



Pp{ri) _ g-iVS4>p(ri) 



(A.2) 



|^g75[3'-0p('-2)-ep(r2)+(T(-r0„(r2)-e„(r2))],g-;;^[-r0p(r2)-ep(r2)+<T(-r0^(r2)-e,.(r2))]j^i + /l C " 

X ^ g;)5['-'0p('-3)-9p('-3)+<T'(r'0„(r3)-e„(r3))]^.^-^[r>p(r3)-6/p(r3)+c7'(r'0„(r3)-e„(r3))]^.^^ 
r' ,a' 

+ g-75[''Vp('-3)-»p(''3)+'^'('-''?i-('-3)-e<,(r3))]^g;i5[r-'0p(r3)-ep(r3)+a'(r>„(r3)-e„(r3))]^.^^ 



with Ti = {xi,UTi) Due to result (5.17i, the surviving terms are those satisfying ^^Ai 



Bi = 0, where and Bi are defined in Eq. (5.17 1. Separating spin and charge parts, we 



81 



have (the common prefactor Ci is omitted in other to lighten notation) 

X! ( [e*^'^"''^''!^ |^g75[-30p(r2)-ep(r2)+(T(0<,(r2)-e<,(r2))]^^-^[</.p(r2)-ep(r2)+<T(0„(r2)-e„(r-2))]^ + i 

+ g-75[3<^p('-2)-ep('-2)+f^(-0.(r-2)-e„(r2))],g^h0p(r2)-ep(r2)+<T(-0„(r2)-6/„(r2))],^i^ 
_g-iy80p,,(ri) |^g;j5[3</.p(r2)-ep(r2)+cT(-0^(r2)-e„(r2))],g-^[-0p(r2)-9p(r2)+(T(-0„(r2)-0„(r2))]^. + i 
_^ g-;75[-30p(r2)-ep(r2)+<T(0<,(r2)-e„(r2))]jg^[</.p(r2)-ep(r2)+a(0„(r2)-e„(r2))],^iy ^ 

^ g7l[''Vp('-3)-ep(r3)+<T'(r>„(r3)-e„(r3))]^^-^[r'0p(r3)-ep(r3)+<T'(r'0^(r3)-0^(r3))]^.^^ 
r'.cr' 

+ g-75[''''^p('-3)-ep('-3)+'^'('-'0-('-3)-0.(r3))],.g^['-'0p('-3)-ep(r3)+fT'(rV.('-3)-e.(r3))]^.^^ \^ 

Now, making expHcitly all the products in the above expression we obtain four different terms, 

^ ^^giy80P,3('-l) |^g;^[-30p(r2)-ep(r2)+fT(0„(r2)-9„(r2))],g-^[0p(r2)-ep(r2)+<T(0<,(r2)-e<,(r2))], 
^ ^g-;7f K('-3)-ep(r3)+<T'(0„(r3)-9„(r3))] ^g^[0p(r3)-ep(r3)+fT'(</,„(r3)-e„(r3))]^.^^ 
_^ ^^giy80P,,('-l) |^g-;^[3<>p(r2)-ep(r2)+<T(-<^„(r2)-e„(r2))],g^[-0p(r2)-ep(r2)+<T(-</.<,(r2)-e<,(r2))], 
^ ^g75[-'^p('-3)-ep(r3)W(-0„(r3)-e„(r3))]^^-^[-0p(r3)-6/p(r3)+t7'(-0„(r3)-e„(r3))]^^^ 
_ ^^g-iy80P,3('-l) ^g;)j[3<^p(r2)-ep(r2)+£r(-0„(r2)-9„(r2))],g-^h0p(r2)-ep(r2)+<T(-0<,(r2)-e<,(r2))], 
^ ^g-^[-0p('-3)-ep(r3)W(-0„(r3)-e„(r3))]^^^[-0p(r3)-ep(r3)+tT'(-0„(r3)-e„(r3))]^^^ 
_ ^^g-iys^p.jlri) |^g-^[-3<#)p(r2)-ep(r2)+a(</,<,(r2)-e„(r2))]^g^[0p(r2)-ep(r2)+<T(</.„(r2)-e„(r2))]j 
^ ^g7lK('^3)-ep('-3)+'T'(</.„(r3)-9„(r3))]^.^-^[</.p(r3)-ep(r3)+o-'(0„(r3)-e„(r3))]^^^ 



Summing over the spins a and a', and using the result shown in Eq. (5.17 1, we finally obtain 
that at zero order in the band curvature (a) and magnetic field (B), the correlator {K^; Ky) 
vanishes, 

= +2e"3Kp'P^i(''i-'-2)+2^P-Pi('-i-'~3)-(Kp--ffp-')Fi(r2-r3)] ^ 

g-i[-2F2(ri-r2)+2F2(ri-r3)]g-[(iK<,+ i/<'-')i=^l(r2-r3)+F2(r2-r3)] 

+ 2e"5[6^p-P'i('-i-''2)+2^p'f'i(''i-'-3)-(^p--^p"')^i(''2-'~3)] X 

g-i[2F2(ri-r2)-2F2(ri-r3)]g-[(iA:„ + iK-l)Fi(r2-r3)-F2(r2-r3)] 
_ g-i[6KpFi(ri-r2)+2XpFi(ri-r3)-(Kp-Kp-l)Fi(r2-r3)] ^ 

g-i[2F2(ri-r2)-2F2(ri-r3)]g-[(iA:„ + iK-i)Fi(r2-r3)-F2(r2-r3)] 
_ g-5[6-ffp-Fi('-i-'-2)+2/<'pFi(ri-r3)-(Kp--R-p-i)Fi(r2-r3)] ^ 

g-i[-2F2(ri-r2)+2F2(ri-r3)]g-[(|-ft:<,+ iif-i)Fi(r2-r3)+F2(r2-r3)] ^ q 
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This can also be verified applying spatial inversion and particle-hole symmetry in correlator 



A.2 Full expression for the correlator {K^] Ky) at first order 
in a, B and 



Here we will repeat the same type of calculations done in Appendix A.l but in this case 
we consider the band curvature term in the correlator {K^; Ky) and we take to first order in 
B the terms depending on the magnetic field {Ky and H.^). Thus the correlator is given in 



Eq. (5.12 1. Again we take the bosonized form for each operator (see Sec. 5.2.2 1 and we put all 



the prefactors together in a constant named C2. The resulting expression is 



gi\/80p(ri) _ g-i\/8(/)p(ri) 



E 



) V2 



(r2)-ep(r2)+(T(-r0„(r2)-e„(r2))],g-^hr0p(r2)-ep(r2)+<T(-r0^(r2)-e<,(r2))]j^i _ ^ ^ 
[r'0p(r3)-8p(r3)+<T'(r>„(r3)~e„(r3))]_-^[r'0p(r3)-ep(r3)+<T'(r'0„(r3)-e<,(r3))]^.^^ 



+ e ^ 



;75['->p('-3)-ep(r3)+<T'(r>,(r3)-e„(r3))]_^[r'0p(r3)-ep(r3)+<T'(r>„(r3)-e„(r3))] 



x(V0p(r4))3) 

There are again four surviving terms, all multiplied by {ieB{x2 + x^)ay/c) C2 



(A.3) 



= ^ ^e*^'^"-"*''!^ ^g^h30p(r2)-9p(r2)+<T(0„(r2)-e„(r2))]„g-^[0p(r2)-0p(r2)+<7(0^(r2)-e<,(r2))]„ + i~j 
^ ^g-7l[0p('-3)-ep('-3)W(0<,(r3)-e<,(r3))]^g;)5[0p(r3)-ep(r3)W(0„(r3)-e„(r3))]_^^ 

_ ^^gi%/80p.„(ri) |^g-;^[30p('^2)-ep(r-2)+<T(-0„(r2)-e„(r2))]„g^[-*p(r2)-9p(r2)+f7(-0<,(r2)-e<,(r-2))],_^ + i 

^g^[-<#'p('-3)-ep(r3)W(-0<,(r3)-e<,(r3))]_^g-;^[-0p(r3)-ep(r3)+a'(-0„(r3)-e„(r3)^^ 
_ ^^g-»%/80p.„(ri) |^g^[30p(r2)-0p(r2)+<T(-0„(r2)-e„(r2))]„g-^[-0p(r2)-0p(r2)+a(-0^(r2)-e,.(r2))]„ + i 
^ ^g-;7f [-0p('-3)-ep('-3)+<T'(-0<,(r3)-e„(r3))]_^g;^h<^p(r3)-ep(r3)+a'(-0„(r3)-e„(r3))]^^^^^ 
_^ ^^g-jy80p.c(ri) |^g-75h30p(r2)-ep(r2)+<T(0.(r2)-e„(r2))]„g^[</.p(r2)-ep(r2)+<T(0<,(r2)-e<,(r2))]„ + i~j 
^ |^g75K('-3)-ep('-3)W(0<,(r3)-e<,(r3))]^g-^[0p(r3)-ep(r3)+<T'(0<,(r3)-e„(r3))]^^^j ^^^^^^^^^gX 



Summing over the spins a and ct', and using result (5.17 1 we write the above expressions in 
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terms of functions Fi and F2 



-H{x2 + x^)a 



4e 



-i[6/fpFi(ri-r2)+2KpFi(ri-r3)-(/fp-/f-i)_Fi(r2-r3)] 



^g-i[-2F2(ri-r-2)+2F2(ri-r3)]g-[(iA:„ + i/f-^)Fi(r2-r3) + F2(r2-r3)] 

_4e~M6^P-^i(''i-''2)+2^P-F'i(''i-''3)-(i^p--ffp"')^'l('-2-r3)]g-i[2F2(ri-r2)-2F2(ri-r3)] 

a 



,-[(iA'„ + iK-i)Fi(r2-r3)-F2(r2-r3)] 



(A.4) 



The functions Fi{r) and F2{r) are given in Eq. (5.18 1. The factor (a/|r4|)'^ is the behavior 
of V(/)p at large distances (a;,ur) 3> a |16|. Because the model studied in Sec. 5.2.2 is 
spin rotational invariant, we have = 1. The function F2{r) gives the angular part in the 
above expression and thus, does not play any role in the scaling analysis, where only power- 
law dependencies are involved. Keeping only the terms with Fi functions and making the 
following change of variables: ri ~ r,r2 — 0, r^, — ri and r4 — 7-2; we obtain result ( |5.20| l. 
We take only the magnitude for 2:3 in the prefactor ^B{x3)ay because, as we said before, 
angular parts can be neglected in the scahng analysis. 



A.3 Full expression for {K^; Ky) with (/2 = 



For the analytical obtention of correlator {K^^Ky) with 92 = 0, we will need in several 



instances the average value (■0](Ti?,(^' ''')^j+i,(yR{y^ 0))o to first order in t^. We note that in 
the absence of magnetic field this quantity equals the free Green's function of the 2D lattice: 



(V']<TK(a;, T)V'j+i,ai?,(y> 0)>o = G^R{y - X, ay, -r) 



i E E ^-«(^' A:±, ic^„)e^'=(^-")e^'=^'^«e*""^ 



kk\ 



(3S 



i^n kk 



iujn — S,R{k) + 2tj_ cos(fc_Laj,) 



kk± . 



1 



1 



— 2t± cos{kj^ay) 

ik{y-x) iuj„T 



1 



J dxidTi GjaR{xi,Ti)gj„R{y ~ x - Xi,-t - Ti). 



oitl) 



^ik{y-x) ^ik^ 



With a finite magnetic field we can use another method: in general, in the presence of a 
time-independent perturbation, the Green's function obeys Dyson's equation: G(r,r' ,uj) = 
Go(r,r', w) + / dridr2Goir, ri, a;)S(ri, r2, cj)G(r2, r', tj). In our case Go(r,r',iu) is the 
Green's function in the absence of interaction and for t± = 0, i.e. Go{r, r' , uj) = 6jj'Q{x — 
x',oj) with j the chain index. On the other hand, the self-energy is simply S(ri,r2,CLi) = 
%'(a:i,a;2,w) = -tj_5{xi ~ a;2)('^j' j+ie"*'"='^i + j-ie*-^*^^!)- The magnetic field B 
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enters in the shift of momentum 6k = eB/c. At first order in t^, we have: 



Gjj'{x,x' ,u!) — Sjj'Q{x — x' + ''^ / dxidx2SjkG{x — xi,uj) X 

ki 

{-t^)5{xi - X2) {Si^k+ie-^'"''' + SLk-ie'""^') %'^(a;2 - x',lj) 
= Sjj'Gix - x' ^uj) - tj_Sj>^j+i J dxiGix - Xi,uj)e~'-^'''''^g{xi - x' ,uj) 

-i_L5j'j-i J dxig{x-xi,oj)e'^'''='g{xi~x',uj). (A.5) 

Thus, we obtain an expression for the Green's function to first order in t± and in presence of a 
magnetic field B, 

(^jVi?(2^'^)V'i+i,<Tfl(y,0))o = Gj+i,j{y,x,-T) = 

-t±^Yl I dxiGjaBiy ~ xi,iLJn)e'^''''^gjaR{xi - x,iujn)e~^'^"^~^'> 

= -t± J dxidTigj„R{y-Xi,Ti)e'-^''''^gj„R{xi-x,-T-Ti). (A.6) 



This last expression reduces to our previous result if i? = 0. The result (A.6 1 will be used in 
the following computation of the correlator {K^; Ky). Starting from Eq. (5.25 i, we write the 
complete expression for correlator {K^, Ky) 



(A.7) 



i{K^-Ky)o = 2e\vt^glay j dx^, J dy^, J 

ja- j'a' b=L,R ^ 



In the first parenthesis the operators are evaluated at position x and time r, while in the second 
parenthesis they are at position y and time 0. We consider the first term: 



"^^^{'^]aR'^],-aR'^3,-<TLll^]''L V'L'feV'l' -a'b^J^-cT',-&V'i' + l,<T^-6>0■ 



(A.8) 



{x,t) 



We can decouple this term in only two different ways. The first is 

e-''''Hi^]aR^r+i..'.-Mi^l-.Ei^f.-.',-b)^{^,,-a^^^^^ (A.9) 

and it clearly implies that <j' ~ a and b = L. The chain indices j and / must be such that the 
fermions are on the same or nearest-neighbor chains, since we work at first order in t±. The 
only solution is j = j', since the other possibility, / ~ j — 1, gives a term of order t^. Hence 
the term reads 

Saa'SbLSjj'e^'^''y{il;^j^j^{x,T)^jjj+i_aR{y,0))ogj,^„R{y-x, -T)gj^^„L{x-y,T)gj„L{x-y,T) 

(A.IO) 
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with G the free propagator Using Eq. (A. 6 1, we find that the first decoupling yields the follow- 
ing contribution to i{Kx; Ky)Q [we temporarily omit the prefactor in Eq. (A. 7 1]: 



—tj_ / dxcLydxidrdri e' 

GjcrRiy - a;i,ri)t/j>fl(xi - x,-t~ ti)^/^ - x,-T)Qj-aL{x - y,T)gjaL{x - y,T) 

= -^^4 E j dxdydxiY.^-''''^y-^"^ X 

iU\ ,11^2 ,i!-'3 

GjaRiy - Xi,Wi)Qj„r{xi - X, iiyi)gj-aR{y " X, ^2) 

y^Qj-crL^x - y,ivi)QjaL{.x - y, Wl + W2 " + i^) 

^ kik2q iu-i -ifz 

x^iL(fc2 - q,iv?,)Q]L{ki + (J, + W2 - it's + i^) 
= -2t^A{ifl,B). (A.ll) 

This corresponds to the diagram shown in Fig. |5.8| The 2nd decoupling of Eq. ( |A.8[ ) is 

e"^''^(V'UV',^-o^-6)o(4_,fi^,'+l,.^-fc)o(^,^-.L^lv'6)o(^ (A.12) 
and implies a' = —a, b = L, and j = j': 

Sa.-a'SbLSjj'e^'^^ygjaRiy-x, -r)('0] ^^^ix, T)i/'j+i.-o-fi(2/, 0))oGj-aL{x-y, T)gj„L{x~y, t). 

(A.13) 

This is the same expression as for the first decoupling provided we change a into —a. It will 
therefore give the same contribution —2t±A{iil, B). 
We now move to the second term: 

-e-'^''y{4^]_^^i:]„^^,,R^,,.,R4'],,,k4'],^_,,^->Pr,--',~b^^^^ (A.14) 



Again there are only two decoupling: 

= -6„^^„:6bR5jfe^'^^y{i]:]^^^^{x, t)V'j+i -CTL(y, Q))oQjaL{y - x, -T)gj„Ft{x ~ y, T)g.j^^„R{ 

= -S„^,dbR6jj>e~'^''ygj-aL{y - x,-T){i(;^j^^{x,T)tljj+i^„L{y,0))oGjaR{x ~ y,T)g.j^^„R{x ~ 

The calculation is very similar to the previous one: just exchange R and L and add a minus 
sign. Hence these two terms contribute 

= 4t_L ^ ^ g]L{ki,ivi)gi^L{ki + 6k,ivi)giL{k2:iv2) 

kik2q ivi,iv2,iv3 

GiR{k2 - q, wz)g]R{ki + q, ivi + W2 - ivs. + i^) 
= 4i_L X! X! G]R{ki,ivi)g^R{ki- 5k,ivi)giR(k2,iv2) 

kik2q ivi,iu2,iv3 

= 4:tj_A{ifl,-B), (A.15) 
where we have used the fact that g^Lik, iui) — gcrR{—k, iut). The third term is 

- e-"'=^(V'],^V'l,_,fl,V'j.-aLl^j<TL V']v'f,'/']' + i,_,'6V'j' + l,-a',-b^j' + l,a',-6)o (A.16) 

^ V V ' 

ix,r) (y,0) 
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and can be decoupled as follows 

-5<T.-<T'<56L<5/J-le~''^'=^(V']<,^V'J'+l,-a^-b)o(V'],_^/^^'J'+l.<T^-b)o(^^^ 

= -Sa-a'SbLSj' j-ier'^''ygj^R{y - X, -T)gj^^„R{y - X, -r)(?/>j t)'0]_i 0))otJj>L(a; - r) 

Both terms are equivalent except for the sign of a. From Eq. \A.5\ we read that 

(V'i<TL(a;,T)'(/'J_i.£,i(?;, 0))o = Gj^j-i{x,y,T) = -t± / dxidTigjaL{x-xi,Ti)e'^^'"''^gjaL{xi-y,T-Ti) 



J- 

(A.17) 

and we then find that these two terms contribute At±A{—iil., B). We continue the game with 
the fourth term: 

e''^^^ {i^],^^L^],L^3aR.i^3,-aRi^],,.b'^]'+i.-a'b^'y+i,-'y' ■-bi^j'+i,'y\-b)o (A.18) 
^ '"^ V ' 

(2:,t) (yfi) 

which has the decoupling 

= Sa-a'SbBSj'j-ie~'^''ygj^^aL{y ^ X, -T)g.j„L{y ~ X, -T)gj^R{x - y, T){'il)j^^„R{x, T)tl;j_^ _^j^{y, 0))o 

= Saa'SbRSj\j^ie^'^''^gj-aL{y " X, ~T)gj^L{y - X, ~t){iPj„r{x, T)-!/;]_-^^^^(y, 0))ogj-aR{x - y, t) 

and contributes —4t±A{~iil., —B). There are four additional terms to consider, for which we 
look only at the first of the two equivalent decoupling: 

■> V V ' 

(y,0) 

= 26a,^a'SbRSjj'e'^''ygjaR{y - X, -T)gj^^„R{y - X, -T){iljj^^^L{x, -aL(2/' 0))ogjaL{x - y, t) 

^ -UjLA{-in,-B), (A.19) 

V V ■■ 

(2;,r) (a,o) 
^ 4<_LA(-jf^, S) (A.20) 

-e'^^''{i}]^Ri}]^-^R^j,-aL^3<yLi}], + i^-^:-bi'], + l^a',-bi'r<y'b^f + l:~'y^ 

" V "-^ V ' 

{^:r) (j/,0) 

= -2<5^,_^'(5bfl5jv j_ie**'=^(V']^^V'i'+i,-<T'b)o(?A],_^flV'jV'b)o(V'j^ 

= -2(5CT-ff'(5bfl<5jVj_ie"'''^t/j^fl(?; - x, -t) _^i^(x, t)'0]_i 0))otJj -^1.(2^ - y,T)gj^L{x - y.r) 

^ Ui_A{in,-B), (A.21) 
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(x,t) iy,0) 

= 2S^^,5bLSfj-ie'^'''"gj-aL{y - a;,-r)(V'jaL(a;,r)V']_i_^i,(y,0))o^'j>fl(a; - y,T)gj_„R{x - y.r) 
^ -4t^A{in,B). (A.22) 

At this point we have calculate all the terms and we have found that 

i{K^; Ky)o = 2e^VT!t±glay{-8t±) [A{ifl, B) - A(iQ., -B) - A{-iQ., B) + A{-iVl, -B)] . 

(A.23) 

It is clear enough that i{Kx; Ky)Q vanishes at B = 0. The first-order term is 

. {K^;Ky}o ^ _^2e''vftlglayB [A'{in) - A'{-in)] + 0{B^). (A.24) 



withA'(if2) = ^dA{iQ,B)/dB\B=o = ^dA{iQ,B)/d5k\sk=o- For i?H this impUes 
The last step is to evaluate A'{u)): 



X 

5fe=0 



eayNyLx f dkidk2dqd^l 1 ^ ii ■ \ 

iui .iu2 

GlR{k2,ii^2)G\,L{k2 - q,iv^Q^L(k\ + q-,Wx + iv2 - W3 + i^)- 
The Matsubara sums are elementary (once third-order poles are under control, and yield 

Fi = In^t), F2 = n^i), ^3 = /(4+)-/(4".+,+4v,-4t) 

where / is the Fermi function /(^) = (e'^^ — At this stage we shift ki and k2 by fcp, and 
q by — 2fcF- Then the energy of the right movers ^^(^1,2) = Vp{ki^2 — fcp) + a(fci,2 — ^f)^> 
which depend only on one momentum, becomes ^ — vpki_2 + a^i 2- the same time 
^ki+q = -VF{ki+q+kp) + a{ki+q+kp)'^hecomes E~{ki + q) = -VF{ki + q)+a{ki + q)'^, 
and finally ^^^-g becomes Cfc2-9+4fcF ~ F~{k2 — q). This last equality holds because we are 
at half-filUng so that 4:kp = 0. We also move to adimensional variables k = (iv^k: 



el/"-- - - 

A'{i^)= -^—^ ^^^^3 j dkidk2dqw{l + 2aki)Fo 



Fi F2 ^ F3 



iCt — d {iCt — a)2 {iCt — d)^ ( 

(A.26) 

where a = and the F's and d are the F's and a with momenta k replaced by k/{Pvp) and 
^ replaced by E. For example, must be replaced by E^^^^^ ^ = + cJfc^-j ^ ^jjg^j. 



the temperature disappears completely from the momentum integrals. We then take away the 
tildes and rewrite: 



A'{in) = 



e vf(3^ 
Fi 



J dkidk2dq (1 + 2(5fci)i^o ^ 



F. 



F. 



5(42) 9i^k,+q) + 9{£k2-q) - 1 - 9{£k,+g)9{£k2-q) 



(A.27) 



Fi 
9{x) 



We treat this expression to first order in a. For this we introduce a new quantity: A"{iil) = 
(iA'{in)/d5i\a^o, in such a way that A'{in) = 2 7^A"(z») + ©(a^). Taking 

into account the various factors we arrive at 



{u + iO+) - A"{-Lj - iO+)] \ . (A.28) 



Taking the imaginary part in Eqs ( |A.28| l and ( |A.27| l leads to delta functions which allows one 
to do the ki integration directly using the identity: 



Im 



1 



1 1 

(/3w-a + iO+)" ~ {-Puj-a-iO+) 

{-iy'S'-'"-^^PLo -a)- <5("-i)(/3c^ + a) 



{n-iy. 

The remaining k2 and q integrals can be done exactly yielding 



i?H(o) = K\i 



-f^Y ['^ (/?c^/4)^-siiV(/?c^/4) 
IQ\ttvv) J Lo tanh(/3a;/4)sinh^(/3tj/4) 



(A.29) 



(A.30) 



We see that the integrand becomes — at high frequency, and therefore the oj integral di- 
verges as — log Wmax wherc uj^^ is the cutoff. We can extract the divergent term by rewriting 
the integral as 



duj (/3w/4)2 -sinh2(/3cj/4) 
tanh(/3w/4)sinh2(/3cj/4) 



dx 



I X — sinh X 
X tanh X sinh^ x 



1 



3 + .T 



^"^'^ dx 



2, + x 



-2.0622 - 2 log ( 1 + ^\ 



-2 log 



Our final result is thus 



i?H(0) ~ i?S 



log 



(A.31) 
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APPENDIX B 



Appendix for the study of the Hall effect on the triangular 

lattice 



B.l Some numerical details 



To obtain the full temperature and density dependence of Ru using Eq (6.6 1 we computed the 
sum over momentum numerically. For this we took a reduced zone of the reciprocal space 
of the 2D triangular lattice, defined by vectors bi = "Tf) ^"'^ ^2 = ^(0, ^) (see 



Sec. 2.3 1. We then performed an x iV grid discretization of this zone (discrete mesh), with 
the total number of sites. Due to the fc-periodicity of coefficients A^, Bk and Ck, and that 



of dispersion relation s^, which are the fc-dependent quantities in Eq (6.61, it can be proved 
that for the description of the whole reciprocal space, is enough to make the sum over a trian- 
gle defined by points Pi = (0, 0); P2 = ^(1, and P3 = ^(1, 0). Each fc-point in the 
triangle must be weighted according to the number of times it appears when the fuU reduced 
zone is recovered, by making reflections of the triangle. In the following, some short routines 
are presented (written in fortran 95). They allow to obtain the diamagnetic terms and therefore 
i?H (in unities of i?H x |e|c with c=l) at ?7 = 0. In this routines, "Nt" denotes the number of 
sites (the total number given by Nt^) and "dist" is the Fermi distribution function (where the 
temperature and density enter). 



function p(Nt) (gives the proper 


weigh to each point in 


the triangle ) 


Implicit None 






integer, intent (in) :: Nt 






integer, dimension ( (Nt/2+2) 


*Nt/2) : : p 




integer : : i, m 






p=4; p(l:Nt-l)=2 






do i=l,Nt/2; m=i* (Nt- (i-2) ) 


; p(m)=2; p (m-1 


)=2 


enddo 






p(Nt)=l ; p ( (2+Nt/2) *Nt/2) = 


1 




end function p 







This funtion p{Nt) is used in the two following routines. 
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function X(ciist) { returns the diamagnetic terms Xx{0) and Xyi^)) 
Implicit None 

real, dimension (:), intent (in) :: dist 
real, dimension ( 2 ): : X 

integer, dimension ( (Nt/2+2) *Nt/2) :: wei 
integer :: k,ki,k:j 
real : : kx, ky 
wei=p(Nt) ; X = zero 
do k=l, size (dist) 

k j=f loor ( (Nt+two) /2-sqrt ( (Nt+two) **2/4-k)+l.d-10) 

ki=k-k j* (Nt-k j ) 

kx=mod (ki , Nt ) *two*pi/ (a*Nt) 

ky= (-mod{ki,Nt) +2*kj) *two*pi/ (a*Nt*sqrt (3 .dO) ) 
X (1) =X (1) +wei (k) * (two* cos (kx*a) + 

(tper/ tpar ) *cos (kx*a/two) *cos (ky*sqrt (3.d0) *a/two) ) *dist (k) 
X (2 ) =X ( 2 ) +wei (k) *cos (kx*a/two) *cos (ky*sqrt (3.d0) *a/two) *dist (k) 
enddo 

X(l)= X(l)/Nt**2 ; X(2)= X(2)/Nt**2 
end function X 



Once we have the diamagnetic terms, we are ready to compute the Hall coefficient i?H by 
means of the following routine. 



function RH (dist) {returns the Hall constant Rh) 




Implicit None 




real (dp) , dimension ( : ) , intent (in) : : dist 




real (dp) , dimension (2) : : diamag 




real (dp) :: RH, terml , term2 , kx, ky 




integer, dimension ( (Nt/2+2 ) *Nt/2 ) :: wei 




integer:: k,ki,kj 




diamag=X (dist ) ; RH= zero; terml=zero; term2=zero; 


wei=p (Nt ) 


do k=l, size (dist) 




kj=floor ( (Nt+two) /2-sqrt ( (Nt+two) **2/4-k)+l.d-10) 




ki=k-k j* (Nt-k j ) 




kx=mod (ki , Nt ) *two*pi/ (a*Nt) 




ky= (-mod (ki,Nt)+2*kj) *two*pi/ (a*Nt*sqrt (3 . dO ) ) 




terml=terml+ 




wei (k) *cos (kx*a) *cos (kx*a/two) *cos (ky*sqrt ( 3 . dO ) * 


a/two) *dist (k) 


term2=term2+ 




wei (k) * (cos (kx*a) +cos (ky *sqrt (3.dO) *a) ) *dist (k) /4 


.dO 


enddo 




RH = - ( (a**2) *sqrt (3.dO) /2.d0) * (terml+ 




(tper/ tpar) *term2 ) / (diamag ( 1 ) * diamag (2 ) *Nt * *2 ) 




end function RH 





In the [/ 7^ case, the above routines must be adapted to obtain the self-energy and 
consequently (rifc), with the different methods explained in Sec. 6.3 Due to the weak momen- 
tum of S(fc, «cli„), we can compute it in a 16 x 16 mesh and then, using the function shown 
below, interpolate it into an x mesh as large as possible (in our work 1024 x 1024): 



92 



function spline (vl,v2,A) {returns a linear interpolation of the array A to the 
vector v) (vl must be in the range [0,size(A,l)] and v2 in the range [0,size(A,2)]. ) 

Implicit None 

real (dp) , intent (in) :: vl, v2 

complex (dp) , dimension ( :, :) , intent (in) :: A 
complex (dp) :: spline 

integer :: N, M, i, j; real (dp) x, y 
complex(dp) :: Al, A2 , A3, A4 
N=size(A,l); M=size(A,2) 
if (vl<zero.or.vl>real (N,dp) ) then 

stop 'Stopped => Out-of-range vector vl in lib: spline' 
else if (v2<zero . or . v2>real (M, dp) ) then 

stop 'Stopped => Out-of-range vector v2 in lib: spline' 
endif 

i=int(vl); j=int(v2); x=vl-i; y=v2-j; i=mod(i,N); j=mod(j,M) 
Al=A(i,j); A2=A (mod (i+1, N) , j) 

A3=A (mod (i+1, N) ,mod( j+l,M) ) ; A4=A (i,mod ( j+l,M) ) 
spline=Al+ (A2-A1) *x+ (A4-A1) *y+ (A1-A2+A3-A4 ) *x*y 
end function spline 



Then, for each element (i, j) of the iVo x A^o original complex matrix icun) we find 
by interpolation the element {iNo /N, jNq /N) with N x N the size of the final matrix we want 
to obtain. In the above routine, Vi, V2 and A correspond to iNo/N, jNq/N and S respectively. 



B.2 Calculation of the DMFT self-energy 

We evaluate the local self-energy in the DMFT framework using the Hirsh-Fye algorithm ||69]| 
as described in Ref |61. In this method the imaginary-time axis [0, /3[ is cut into L slices, and 
the Trotter formula is used in each time slice in order to single out the Hubbard interaction. 
In a second step the interaction is decoupled via the introduction of an Ising variable in every 
time slice. The Green's function Q{t) is finally calculated by averaging over the ensemble 
of configurations of the Ising variables using a Monte-Carlo sampling and local updates. In 
our calculations at /? = 1 and U ^ 20 we take L = 128 and we keep 10^ out of the ^ 10* 
configurations visited. The numerical accuracy of the calculated Q{t) is estimated to be ~ 
10^'^ at the highest U values, and closer to 10^^ at [/ < 8. We believe that the accuracy at 
J7 > 8 is not sufficient to get a reliable self-energy; hence we did not evaluate i?H at ?7 > 8 in 



Fig. 6.8 



In order to calculate the self-energy and solve the DMFT self-consistency condition, Eq. ( 6. 19 1 
we need to Fourier transform Q(t) from imaginary time to imaginary Matsubara frequencies 
iojn- In traditional implementations of the algorithm this step is performed through a cubic 
spline interpolation of G{t). Because cubic splines are non-analytic, however, the resulting 
Fourier series are unreliable at frequencies above ^ L/f3. Instead of an interpolation, we have 
performed a fit of G{t). Our fitting function is a discrete form of the spectral representation, 
G{t) = — J ds A{e)e~'^'^ f{—s), which we express as 

M 

g{r) = -J2A,e-'^^f{~e,) (B.l) 
with Aj ^ and jii = 1- The number M of poles Sj and their weight Aj are determined 



by adding more and more terms in Eq. ( |B.1[ ), until the fitted function matches all QMC data 
points within a numerical tolerance, which we take as the estimated accuracy of G{t). The 
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Figure B.l: DMFT imaginary-time propagators calculated at /3 = 1 and n = 1.54 for 
various interaction strengths. The symbols show the QMC results on the discrete-time 
mesh. The solid lines are the fit, Eq. ( |B.1[ ), used to evaluate the self-energy shown in 
Fig.lSTTl 



Fourier transform is then simply 



M 

E 



lUJn — £i 



(B.2) 



The calculated self-consistent propagators G{t) and GoiT) displayed in Fig. B.l together 



with the fits to Eq. (B.l 



Solving Eq. ( 6.19| l at fixed electron density requires to determine the chemical potential /i 
self-consistently. In our calculations we perform the search for both the self-consistent Go and 
/i in one shot using a global minimization procedure. As a result the self-consistent solution 
can be reached in typically less that 20 iterations. 



B.3 Self-energy and distribution function in the atomic limit 



In the case U ^ t one can treat the Hamiltonian Eq. (6.1 1 using a perturbative expansion in 
tij/U. The atomic limit is the zeroth-order term of this development, and it corresponds to a 
collection of disconnected sites with four possible states on each site. This limit is not very 
useful since there is no hopping and thus no possible transition below the Hubbard energy U. In 
order to retain the low-energy dynamics of the problem we adopt an hybrid approach, where the 
free dispersion is used in the lattice Green's function together with the self-energy evaluated in 



the atomic limit. The atomic self-energy is obtained by diagonalizing the Hamiltonian Eq. ( 6. 1 1 
with tij = 0, which leads to the atomic Green's function 



Gat(ja;„) 



l-n/2 

iuJn + Mat 



n/2 



Mat - ' 



(B.3) 



while the non-interacting Go, at =1/ {i^n + Mat) results by putting J7 = 0. Fr om D yson's equa 



tion, E = (S'o - we deduce the atomic self-energy displayed in Eq. (^6. 20b. Here //at is 



the chemical potential in the true atomic limit — i.e. the limit where the lattice Green's func- 
tion takes the form (B.3 i, and therefore the electron density is given by n (2 ~ n) f {— ji-^t) + 
nf{U — Mat) with / the Fermi function. We can invert this relation and express /iat explicitly 
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in terms of the electron density as 



Mat 



^-1 + ^(^-1)' + ^-/^^ (1-1) 



Using the atomic self-energy Eq. ( |6.20[ ) as an approximation to the exact self-energy in Eq. ( |6.14p , 
we evaluate analytically the lattice distribution function (n^). Let's first remark that 



with 



lUJn 



K = 

Zk = 



S.k - Sat(«t^,i) iUJn - El 



Afc + l7-/iat)/2 
^fc + Afc + [/ - /iat 

4A, 



(Cfc + Afc - t/ + /iat)(A^at - t/ + n?7/2) + nii^JJ 
^k{t-^,t-U + nU/2) + nfi,,U/2 

Afe = v/(Cfc + + Mat)2 + 2(n - 2)(efc + Mat)C/. 
As a result the Matsubara sum in Eq. ( |6.14[ ) is easily performed to yield 

{nk).,, = Zkf{E+) + {l-Zk)f{E^). 



(B.4) 



Within this approximation it is also straightforward to perform the infinite U limit. Taking 
into account that both /i and /^at are either of order i (if n < 1) or of order [/ (if n > 1) we find 
that Zk approaches n/2 as [/ increases toward +oo. Likewise, if n < 1 we have E^ ^ U and 
Er ~ t while if n > 1 we have Et ^ t and E, 



k 



-U. Hence we find 



U=oo 



(1- t)/[(l-f)a-fMat] n<l 



n f\ n ? 



1-t M 



^ 2 



n > 1 



where we have introduced flat = Hat — U = — log[(l — n/2) /{n — 1)] and S,k = £k — fl with 
fl ~ fi — U. For the purpose of evaluating the high-temperature behavior of the Hall coefficient 
atU — oo, we finally expand the distribution function in powers of f3 following the procedure 
described in Sec. 16.3.11 



\U=co 



f -n(l-n)efcf +0(/32) n<l 
f ~(n-2)(l-n)£fcf +0(/32) n>l 



Comparing with Eq. (6. 11 1, which is valid at ?7 = 0, we see that the only difference between 
the high-temperature behaviors of Ru at [/ = and [/ = oo is the n-dependent slope, and we 
easily deduce that 



T/t 



U=oa 
H 



(T> t) 



e n(l — nj 
T/t 1 



4 2+(t'/t)2 
a^\/3 3 



e (n-2)(l-ri) 



2+(t'/t) = 



n < 1 
n > 1 



By introducing 5 = |n — 1| which measures the doping with respect to half-filling, these two 
cases can be recast in one single expression shown in Eq. ( |6.21| i. 
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